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1 Introduction 



Supersymmetric extensions of the Standard Model (SM) are motivated by a solution of 
the hierarchy problem [1-5], an automatic unification of the running gauge couplings at a 
Grand Unified (GUT) scale Mqut [6-9], and the possibility to explain the dark matter relic 
density in terms of a stable neutral particle [10, 11]. 

It is well known that a supersymmetric extension of the Higgs sector of the SM [12, 13] 
requires the introduction of two Higgs SU (2)-doublets and Hd, where vacuum expecta- 
tion values (vevs) of Hu and Hd generate masses for up-type quarks and down-type quarks 
and charged leptons, respectively. The model with this minimal field content in the Higgs 
sector is denoted as the Minimal Supersymmetric Standard Model (MSSM) (for reviews see, 
e.g., [14-16]). The Lagrangian of the MSSM must contain a supersymmetric (SUSY) mass 
term /i for and H^, which has to be of the order of the SUSY breaking scale Msusy for 
phenomenological reasons (see below). This spoils a potentially attractive property of su- 
persymmetric extensions of the SM: the electroweak scale generated by the Higgs vevs could 
depend only on Msusy, which would be the only scale asking for an explanation to why it 
is far below Mqut ot the Planck scale Mpianck- The question how a supersymmetric mass 
parameter fi can assume a value of the order of Msusy is denoted as the "yU-problem" [17] 
of the MSSM. 

A simple and elegant way to solve this problem consists in generating an effective (su- 
persymmetric) mass term /i in a way similar to the generation of quark and lepton masses 
in the SM: the mass term fi is replaced by a Yukawa coupling of Hu and Hd to a scalar 
field, and the vev of the scalar field - induced by the soft SUSY breaking terms - is of the 
desired order. Since the fi parameter carries no SU{3) x SU{2) x U{1)y quantum numbers, 
the field to be introduced has to be a singlet 5* (the complex scalar component of a chiral 
superfield S), and the resulting model is the Next-to-Minimal Supersymmetric Standard 
Model (NMSSM), sometimes also denoted as the (M+1)SSM. 

In fact, already the first attempts to construct supersymmetric extensions of the SM 
employed such a singlet field [12,13,18]. A singlet was also present in most of the first 
globally supersymmetric GUT models [5,19-22]. Then one realised that spontaneous su- 
persymmetry breaking in the framework of supergravity (SUGRA) leads in a simple way 
to the desired soft SUSY breaking terms in the Lagrangian; see [14] for an early review. 
Within SUGRA, a /i term of the order of Msusy can actually be generated if one assumes 
the presence of a particular Higgs-dependent structure in the Kahler potential [23]. Still, 
the first locally supersymmetric extensions of the SM [24-26] as well as most GUT models 
within SUGRA [27-34] used a singlet field in the Higgs sector leading to variants of the 
NMSSM at the weak or SUSY breaking scale < 1 TeV. (See also SUGRA models motivated 
by string theory [35-44].) 

Expanding around the vacuum with non- vanishing vevs of the neutral CP-even compo- 
nents of Hu, Hd and S, one finds that the scalar components of S mix with the neutral scalar 
components of Hu and Hd leading, in the absence of complex parameters (corresponding 
to the absence of explicit CP violation), to three CP-even and two CP-odd neutral scalars 
(see [45-47] for some reviews). Likewise, the fermionic superpartner of S mixes with the 
neutral fermionic superpartners of Hu, Hd (and the neutral electroweak gauginos) leading 
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to five neutralinos. As a consequence, both the Higgs and the neutrahno sectors of the 
NMSSM can get considerably modified compared to the MSSM. 

In the Higgs sector, important alterations with respect to the MSSM are a possibly larger 
mass of the Higgs scalar with SM-like couplings to gauge bosons, and additional possibly 
light states with reduced couplings to gauge bosons. Notably a light CP-odd scalar with 
vanishing couplings to two gauge bosons like all CP-odd scalars (but with possibly even 
enhanced couplings to quarks and leptons) can appear in the Higgs spectrum, allowing for 
new Higgs-to-Higgs decays. Under these circumstances, the detection of Higgs bosons at 
colliders can become considerably more complicated. At present it is not even guaranteed 
that a single Higgs scalar can be observed at the LHC within the NMSSM, see Section 5. In 
addition, a light CP-odd scalar can affect "low energy" observables in B physics, T physics 
and the anomalous magnetic moment of the muon. 

The modifications within the neutralino sector are particularly relevant if the additional 
singlet-like neutralino is the lightest one and, simultaneously, the lightest supersymmetric 
particle (LSP). This would have an important impact on all decay chains of supersymmetric 
particles (sparticles), and hence on their signatures at colliders. For instance, the next-to- 
lightest supersymmetric particle (NLSP) can have a long life time leading to displaced 
vertices. Also, the LSP relic density has to be reconsidered in this case. 

Hence, apart from the theoretical motivations for the NMSSM, its phenomenological 
consequences must be worked out in order not to miss (or misinterpret) both Higgs and 
sparticles signals - or the absence thereof - at past, present and future colliders. 

In the present paper we review theoretical and phenomenological aspects of the NMSSM: 
Higgs masses and couplings including radiative corrections, issues related to the NMSSM 
Higgs potential as the nature of the electroweak phase transition, Higgs and sparticle spectra 
within specific assumptions on the origin of supersymmetry breaking as minimal SUGRA 
and gauge mediated supersymmetry breaking (GMSB), NMSSM specific Higgs and sparti- 
cle signatures at colliders, possible impact on low energy observables, and the dark matter 
relic density and its detection. Possible variants of the simplest NMSSM as explicit CP vi- 
olation in the Higgs sector, R-parity violation (in connection with neutrino masses) and 
extra U{iy gauge symmetries are sketched briefly. Relevant formulae like Feynman rules, 
renormalisation group (RG) equations and details of the radiative corrections to the Higgs 
masses are given in the Appendices. 

Clearly, we cannot present all details of all results that have been obtained within the 
NMSSM up to now. However, in all Sections we attempt to reference to the complete 
available literature, where the various subjects are discussed. 

Let us conclude the Introduction by recalling the arguments for a fi parameter of the 
order of MgusY, whose necessity constitutes the main motivation for the NMSSM: both 
complex Higgs scalars and of the MSSM have to be components of chiral superfields 
which contain, in addition, fermionic 5'[/(2)-doublets tpu and tpd- The Lagrangian of the 
MSSM can contain supersymmetric mass terms for these fields, i.e. identical positive masses 
squared /i^ for l-ff^P and and a Dirac mass fi for ipu and ipd- In the presence of a 

SUSY mass term ~ in the Lagrangian, a soft SUSY breaking mass term BfiHuHd can 
also appear, where the soft SUSY breaking parameter B has the dimension of a mass. 

For various reasons the mass parameter fi cannot vanish. First, a Dirac mass /i for ipu 
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and ipd is required for phenomenological reasons: both fermionic 5'[/(2)-doublets ijju and ijjd 
contain electrically charged components. Together with with the fermionic superpartners of 
the W"^ bosons, they constitute the so-called chargino sector (two charged Dirac fermions) 
of SUSY extensions of the SM. Due to the fruitless searches for a chargino at LEP, the 
lighter chargino has to have a mass above ~ 103 GeV [48]. Analysing the chargino mass 
matrix, one finds that this lower limit implies that the Dirac mass fi for tpu and ipd ^ ior 
arbitrary values of the other parameters - has to satisfy the constraint > 100 GeV. 

Second, an analysis of the Higgs potential shows that a non-vanishing term BfiHuHd 
is a necessary condition for that both neutral components of and if^ are non-vanishing 
at the minimum. This, in turn, is required in order to generate masses for up-type quarks, 
down-type quarks and leptons by the Higgs mecanism. Moreover, the numerical value of 
the product Bfi should be roughly of the order of the electroweak scale (M|). 

Third, fi = would generate a Peccei-Quinn symmetry in the Higgs sector, and hence 
an unacceptable massless axion [17]. 

However, |/i| must not be too large: the Higgs potential must be unstable at its origin 
Hu = Hd = in order to generate the electroweak symmetry breaking. Whereas the 
soft SUSY breaking mass terms for and Hd of the order of the SUSY breaking scale 
MsusY can generate such a desired instability, the //-induced masses squared for and Hd 
are always positive, and must not dominate the negative soft SUSY breaking mass terms. 
Consequently the n parameter must obey < Msusy- Hence, both "natural" values /i = 
and very large /i (~ Mqut or ~ Mpianck) are ruled out, and the need for an explanation of 
ji ^ Msusy is the //-problem. 

Within the NMSSM, where /t is generated by the vev {S) of a singlet 5", {S) has to be 
of the order of Msusy; this is easy to obtain with the help of soft SUSY breaking negative 
masses squared (or trilinear couplings) of the order of Msusy for S. Then, Msusy is the only 
scale in the theory. In this sense, the NMSSM is the simplest supersymmetric extension of 
the SM in which the weak scale is generated by the supersymmetry breaking scale Msusy 
only. 
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2 Lagrangian of the general NMSSM 



In this Section, we present the scalar potential (at tree level) and the mass matrices in the 
general NMSSM. Their diagonalization and conventions for the mixing matrices are deferred 
to the Appendix A, where we also list the Feynman rules for the physical eigenstates. 

2.1 Tree level potential and mass matrices 

As in any softly broken supersymmetric theory, the Lagrangian of the NMSSM is specified 
by the supersymmetric gauge interactions, the superpotential and the soft supersymmetry 
breaking gaugino masses, scalar masses and trilinear couplings. 

To begin with, we consider the general NMSSM defined as the MSSM with an additional 
gauge singlet chiral superfield S, including the most general renormalisable couplings in the 
superpotential and the corresponding soft SUSY breaking terms in £soft- (Here we limit 
ourselves, however, to the R-parity and CP conserving case. Chiral superfields are denoted 
by hatted capital letters; unhatted capital letters indicate their complex scalar components.) 
In the general NMSSM, the terms in the superpotential H^Higgs depending exclusively on the 
Higgs superfields Hu, Hd and S are (here we follow the SLHA2 conventions in [49] where, 
however, Hu is denoted as H2, and as Hi): 

Wffiggs = (/u + \S) Hu-Hd + ^fS + + (2.1) 

where the terms ~ A, k are dimensionless Yukawa couplings, the terms ~ /i, ^' are su- 
persymmetric mass terms, and of dimension mass^ parametrizes a (supersymmetric) 
tadpole term. To (2.1) we have to add the Yukawa couplings of the quark and lepton 
superfields: 

W^^ukawa = huQ ■ HuU% + hdH(i- Q -D^ + heHd ■ L E'}^ (2.2) 

where the Yukawa couplings hu, hd, he and the superfields Q, f/^, D^, L and E'^ should be 
understood as matrices and vectors in family space, respectively. 
In (2.1) and (2.2), the SU{2) doublets are 

and the products of two SU{2) doublets are defined as, e.g., 

Hu-Hd = H^HJ - HlH'd . (2.4) 

The corresponding soft SUSY breaking masses and couplings are, again in the SLHA2 
conventions [49], 

-^soft = m]jJ^Hu\^ + m]jJ^Hd\^ + ml\S\'' + ml\Q^\+ml\Ul\ 
+ml\Dl\+ml\L^\+ml\El\ 

+{huAu Q-HuU'r- hdAd Q ■ Hd D% - KA, L ■ Hd E% 

+XA^ Hu-HdS+ IkA^ S'' + mlHu-Hd+ ^m'| + S + h.c.) . (2.5) 
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(Sometimes the definitions m| = S/i, m'J = B' ^' are used.) 

Clearly, the dimensionful supersymmetric parameters /i, /i' and in the superpoten- 
tial (2.1) (and the associated soft SUSY breaking parameters mg, m'J and ^5 in (2.5)) have 
to be of the order of the weak or SUSY breaking scale, in contradiction to one of the theo- 
retical motivations for the NMSSM mentioned in the Introduction. Although some of these 
terms are non-vanishing in various scenarios, one considers mostly the simpler NMSSM 
with a scale invariant superpotential where fi = fi' = = 0, 



Wscinv. = A5i/,-i7, + ^53 (2.6) 

and vanishing parameters ml, m'J and in (2.5). Then, a vev s of S" of the order of the 
weak or SUSY breaking scale generates an effective //-term with 

fies = As , (2.7) 

which solves the /i-problem of the MSSM. 

As any supersymmetric theory with a scale invariant (cubic) superpotential, the com- 
plete Lagrangian - including the soft SUSY breaking terms - specified by (2.6) possesses 
an accidental Za-symmetry corresponding to a multiplication of all components of all chiral 
superfields by a phase e^'^*/^. In the following we denote the version with the scale invariant 
superpotential (2.6) as the "Zs-invariant NMSSM". Any of the dimensionful terms in the 
general superpotential (2.1) breaks the Zs-symmetry explicitly. Subsequently, the version 
corresponding to the general superpotential (2.1) will be denoted as the "general NMSSM". 
In the literature, "NMSSM" stands mostly for the Zs-invariant NMSSM. In the following 
we will retain this convention. 

The Higgs sector of the Zs-invariant NMSSM is specified by the seven parameters 

A, K,, m^^, mjj^, m|. Ax and . (2.8) 

Expressions for the Higgs mass matrices in the Za-invariant NMSSM can be found, e.g., 
in [45,46,50-52]; in the following we discuss, for completeness, the general NMSSM from 
which the Za-invariant NMSSM can always be obtained by setting m| = m'| = = fJ' = 

li' = iF = 0. 

From the SUSY gauge interactions, the F- and the soft SUSY breaking terms one obtains 
the Higgs potential: 

V^Higgs = \\{HtH, -HlH',)+KS^ + ^'S + iF\^ 

+ [ml^ + l/i + \S\') + \Ht\') + + l/i + A^l^) + \H^ 



|2 



d 



+ml\S\' + {\A, {H^HJ - HM) S + ^nA^ + mj {HtH~, - HlHf) 
+\m'IS^ + isS + \i.c) (2.9) 
where gi and g2 denote the U{1)y and SU{2) gauge couplings, respectively. 
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The neutral physical Higgs fields (with index R for the CP-even, index I for the CP-odd 
states) are obtained by expanding the full scalar potential (2.9) around the real neutral 
vevs Vu, Vd and s as 

Hu = Vu + ^ , H^ = Vd + -j= , S = s + — -j= — ; (2.10) 

where the vevs have to be obtained from the minima of 

^Higgs = {->^VuVd + KS^ + fi's + ^fY + g {vl - vf) ^ 

+ ("^k + (/" + As)') vl + {ml^ + (/i + Xsf) 



-m 



^ - 2XAx VuVdS + -kA^ - 2ml ^uVd + m'^ + 2^5 s , (2.11) 

The signs of some parameters in a Lagrangian have no physical meaning, since they can 
be changed by field redefinitions — ?■ — </>. Analysing all possible field redefinitions in the 
Lagrangian above, one finds that one can choose positive Yukawa couplings A, ht-, ht, hj. 
(the latter corresponding to the (3,3) components in family space of h^-, hd and in (2.2)) 
and positive vevs Vu and Vd, whereas k, s and all dimensionful parameters can have both 
signs. In addition, in the general NMSSM one of the dimensionful S'-dependent parameters 
can be removed by a constant shift of the real component of S. In particular, one can put 
the MSSM-like /i-term in the superpotential (2.1) to zero by a redefinition s s — fi/X 
(assuming A 7^ 0) and corresponding redefinitions of the other dimensionful parameters; in 
the following we will assume, for simplicity, that this convention is used (which does not 
imply, in general, that m| = !). 

We define, as usual, 



tan/3 = — (2.12) 

Vd 



(positive by assumption), and we have 



M|=^?V where = ^1±M , = vl + vj GeVf . (2.13) 

Furthermore it is convenient to define, together with fics as in (2.7), 

B^s = Ax + Ks, ml = ml + \{fi's + ^F) (2.14) 

where we have used the convention fi = 0. Bcs plays the role of the MSSM-like 5-parameter, 
and m| vanishes in the Zs-invariant NMSSM. 

It is possible to use the three minimisation equations of the potential (2.11) with respect 
to Vu, Vd and s in order to replace the three parameters m"^^, and m| by Vu, Vd and s. 
The minimisation equations are given by 



Mml^ + /ie\ + >?vl + '-^{vl -vl))- Vd (/iefri^cff + mi) = 



„2 I „2 



+ + >^'< + '-^K - <) j - Vu (/.cffBcfT + ml) = , 
s(ml + + Ai'^ + 2<F + ^iA^s + 2k^s'^ + Snsfi' + X'^{vl + vj) - 2\KVuVd 

+is + iF^' - Xv^Vd{Ax + /^') = . (2.15) 
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From the first two of these equations one can derive 



hence, in order to have both f„ and different from zero (tan/3 7^ 0, 00), we need 
ficsBcs + ml 7^ in the general NMSSM, and /iefr-Ses 7^ in the Zs-invariant NMSSM. 

Given Mz, one can choose as six independent parameters in the Higgs sector of the 
Za-invariant NMSSM 

A, K, Ax, A^, tan/3, fi^s, (2.17) 
to which one has to add in the general NMSSM (in the convention n = 0) the five parameters 



m 



3' 



/w', m's, and ■ (2.18) 



Let us consider the conditions arising from a phenomenologically acceptable minimum 
of the potential: both and must not vanish, and - in the absence of a /i-term - s 
must be large enough to generate a sufficiently large effective /i-term /ieg = As > 100 GeV 
(see the Introduction). In the Zs-invariant NMSSM, the dominant terms for large s in the 
potential (2.11) are 

2 

^Higgs(s) ~ m| + -kA^ + K^s^ . (2.19) 
One easily finds that Al > 8m| is a condition for s 7^ [53,54], and 

Al>9ml (2.20) 
a condition for an absolute minimum with 

s ^ i- (^-A^ - ^Al - 8m|) . (2.21) 

We note that in the case of a scale invariant superpotential (2.6) and without any soft 
terms for the singlet - as it can happen in GMSB models [55] - a vev s 7^ is still triggered 
by the s-dependent terms in the potential (2.11) neglected in (2.19). However, the resulting 
value of s is too small in order to give /icfr > 100 GeV [55]. 

Depending on the parameters, the Higgs potential of the Zs-invariant NMSSM can pos- 
sess several local minima (see, e.g., [56]). Notably, one should verify whether vacua where 
one of the vevs f^, Vd or s vanishes (typically preferred during the cosmological evolution) 
are not deeper. A general analysis (taking the radiative corrections to the effective potential 
into account) is quite involved, but typically one obtains upper bounds on n at least of the 
type < A^ [54]. 

The tree level Higgs mass matrices are obtained by expanding the full scalar poten- 
tial (2.9) around the real neutral vevs Vd and s as in (2.10). Then, the elements of the 
3x3 CP-even mass matrix read in the basis {HdR, Hur, Sr) after the elimination of 
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rrijj^, Tn^Hu ^^"^ "^1 (still in the general NMSSM, but assuming /i = 0) 

■^1,11 = a'^^d + (AicflF 5cfr + ml) tan (3 , 

-^5,22 = 9^^l + (/"eff 5cff + ml)/ tan /3 , 

•M|,33 = KA^ + ^l')'^ + ^^s{A^ + AKs + ?>^^')-{is + ^F^^')/s, 

-^1,12 = (2A^ - g^)vuVd - /icff 5cff - rhl , 

-^1,23 = A(2/icff Vu - (5eff + + /iO^^d) . (2.22) 

Rotating the upper left 2x2 submatrix by an angle /3, one finds that one of its diagonal 
elements reads ^ 

M| (^cos^ + ^ ^^^^ ^'^) (2-2^) 

which constitutes an upper bound on the lightest eigenvalue of A^|. The additional positive 
contribution ~ sin^ 2/3 (as compared to the MSSM) in the NMSSM is highly welcome 
in view of the present lower LEP bound of ~ 114 GeV on the mass of a Higgs scalar with 
SM-like couplings to gauge bosons [57]. However, this additional contribution is relevant 
only for not too large tan/3; in fact, the expression inside the parenthesis in (2.23) is larger 
than one only for > g"^, in which case it is maximal for small tan /3. Moreover, the 
actual lightest eigenvalue of is smaller than the value given in (2.23) in general, see 
the discussion in Section 3.2. 

In the general NMSSM, the elements of the 3x3 CP-odd mass matrix A^'p read in the 
basis {Hdi, H^i, Si) 

p,ii = (/"eff B^s + ml) tan/3 , 
p,22 = (Aicff ^cff + ml)/ tan (3 , 

A4p33 = A(5cff + 3Ks + /i')^-3KA^s-2m'|-K/i's-^F(4K + — ) - — , 

s \ s J s 

■M'p,12 = l^cS ^eflf + ml , 

-Mp23 = \vd{Ax - 2ks - fi') . (2.24) 



A^'p contains always a massless Goldstone mode G. Next we rotate this mass matrix into 
the basis {A, G, Si), where A = cos (3 Hui + sin /3 H^i'- 







nul 


H 


\ s, j 






(2.25) 

Dropping the Goldstone mode, the remaining 2x2 mass matrix A^p in the basis [A, Si) 
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has the elements 



Alpii — — 



sin 2(3 

A4p22 = A(i?cfr + 3ks + fi') ^"^'^ — SkAkS — 2me — k/x's — f 4k + — ) — — , 

s \ s J s 

Ml^^2 = XiAx - 2ks - fi') V . (2.26) 
In the Zs-invariant NMSSM, Mp simphfies to 

2 /icff Bcs 



Mi,, 



sin 2/3 



VuVd 



Mp22 = A(i?cff + 3/ts) 3kAk S 



s 



Mp^,2 = HAx-2ks)v (2.27) 

with i?ofj as given in (2.14). The matrix element A^pn corresponds to the mass squared 
Ml of the (only physical) CP-odd scalar A of the MSSM. 

Finally the charged Higgs mass matrix in the basis {H^, H^* = H^) is given by 

M'l = (/.e. B,, + ml + vM§ - A^)) ( ^ ,2 /3 ) • (^.28) 
It contains one massless Goldstone mode, and one eigenstate with mass 

Ml = 2(^effi?eff + mi) gj _ 

^ sin 2/3 ^2 ' ^ ^ 

Due to the term ~ A^, the charged Higgs mass in the NMSSM can be somewhat smaller 
than in the MSSM (for a given value of M| = Mj,,,). In contrast to the MSSM it is not 
even guaranteed within the NMSSM that f/(l)cm remains unbroken: considering again the 
Za-invariant NMSSM where rrig = 0, the expression for the charged Higgs mass squared 
becomes negative for s = fics = 0, > (7|/2, indicating a possible minimum in field space 
where the charged Higgs has a vev. Although radiative corrections have to be added and 
the depth of this minumum has to be compared to the physical one with s 7^ 0, A is bounded 
from above by the absence of a charged Higgs vev. 

The diagonalization of all scalar mass matrices is carried out in Appendix A (together 
with the Feynman rules); next we consider the fermionic sector. 

First, we have to consider the soft SUSY breaking gaugino mass terms, which do not 
differ from the MSSM. Denoting the U{1)y gaugino by Ai, the SU{2) gauginos by Ag 
{i = 1, 2, 3) and the SU{3) gauginos by A3 (a = 1, ... , 8), the soft SUSY breaking gaugino 
mass terms in the Lagrangian read 

C = ^MiAiAi + ^MsA^A'^ + ^MgA^A^ . (2.30) 
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In the neutralino sector, Ai and A| mix with the neutral higgsinos ip^^ ip^, ips and generate 
a symmetric 5x5 mass matrix JHq. In the basis = (— iAi, —iX^, ip^ ip^, ips), the resulting 
mass terms in the Lagrangian read 



■(^AY-Mo(V^°)+h.c. 



(2.31) 



where 
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and the term fi' in the (5,5) element appears in the general NMSSM only. 
As in the MSSM, the charged SU{2) gauginos are A" = ^ (A2 + iA2) 



(2.32) 



and A"*" 



V2 



{XI 



^Xl) 



which mix with the charged higgsinos ip^ and ip^ . Defining 



%p^ 



-iX+ \ ( -iX 

the corresponding mass terms in the Lagrangian can be written as 

1 



C 





X 



+ h.c. 



(2.33) 



(2.34) 



with 



X 



(2.35) 



M2 g2Vu 

92Vd /icff 

Again, the diagonalization of the neutralino and chargino mass matrices will be described 
in Appendix A. 

Finally we give the top and bottom squark and r slepton mass-squared matrices, t^, tf(, 
hi and denote the scalar components of the third generation quark superfields L^^g, t/^g, 

and D/jg ; and tl the scalar components of the third generation lepton superfields 
Z/L3 and -Eig, and tr the scalar component of the third generation lepton superfield £^^3. 
(The superfields were given in (2.2) and (2.3).) m\. = ^lizi = ""^Ds' ""^Qs' ^'s-i 

are the soft SUSY breaking masses squared for the third generation (see £soft in (2.5)), 
assumed to be diagonal in family space. 

The top squark mass matrix reads in the basis [tn, ti)'- 



ht{AtVu 



^d) 3 



ht{AtVu - ^icsVd) 



171' 



Q3 



+ hjvt + ivi 



12 



9l 



(2.36) 



The bottom squark mass matrix reads in the basis (6/j, hiP] 



ml + hM + K 

h{AbVd - fJ^eSVu) 



m, 



9? , 92 
12 4 



(2.37) 
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The tau slepton mass matrix reads in the basis (r/j, tl): 

2 



(2.38) 



The tau sneutrino (i'tl) mass squared is: 

ml-{vl-vl)('-l±^] . (2.39) 



Herewith we conclude the presentation of the tree level Lagrangian; clearly one has to 
add radiative corrections to the Lagrangian and all the resulting mass matrices. For the 
couplings/Feynman rules we refer to Appendix A. 

2.2 Limiting cases: the effective MSSM and approximate global 
symmetries 

2.2.1 The effective MSSM 

As it becomes clear from the superpotential (2.1) of the general NMSSM, all couplings 
between the components of the singlet superfield S and the MSSM Higgs superfields Hu 
and Hd vanish for A — ?■ 0. In order to generate a reasonably large value for fi^s (in the 
absence of an MSSM-like //-term in the superpotential), one should keep As > 100 GeV in 
this limit. From (2.21) one finds that the vev s scales as 1/k, hence a reasonable decoupling 
limit is 

A~K^O, s~1/k^oo, (2.40) 

while keeping all dimensionful parameters fixed. 

The corresponding parameters of the effective MSSM can be easily deduced from (2.7), 
(2.14) as well as the 1-2 components of the CP-even and CP-odd Higgs mass matrices (2.22) 
and (2.24) (and the charged Higgs and neutralino mass matrices): 

/ieff = As , 

"^Lff = ^1 + f^cs{Ax + Ks) . (2.41) 
In the Za-invariant NMSSM frig vanishes, and one has 

B,s = Ax + ^s. (2.42) 

In the limit (2.40) all couplings between the MSSM sector (including quarks, leptons 
etc.) and the CP-even, CP-odd and fermionic singlet states vanish; it seems a priori 
impossible to distinguish the NMSSM from the MSSM, since the singlet-like states would 
never be produced. However, in this limit the singlino-like neutralino X5 with its mass given 
by (A^o)55 = 2fi;s (in the Zs-invariant NMSSM) can easily be the LSP. Then, assuming R- 
parity conservation and a small non- vanishing value of A, all sparticle decay cascades will 
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first proceed as in the MSSM into the MSSM-hke NLSP (which could be a charged slepton!) 
which, at the end, will decay into X5 + SM-particles. The final decay of the MSSM-like 
NLSP involves necessarily a coupling ~ A, implying a possibly very long lifetime of the 
NLSP leading to displaced vertices [58]. Then, the difference between the NMSSM and the 
MSSM can be spectacular even in the decoupling limit (2.40). 

2.2.2 The Peccei-Quinn symmetry limit 

If the term A S Hu ■ Hd would be the only S dependent term in the superpotential of the Z3- 
invariant NMSSM, the Lagrangian would be invariant under a Peccei-Quinn-like symmetry 

Hu^Hue'^^^Q, Hd^ Hdc''^^'^, S ^ S e-^'"^"-^ (2.43) 

which allows to solve the strong CP-problem [59,60]. Since this global symmetry is spon- 
taneously broken by the vevs v^, Vd and s, a massless Nambu-Goldstone boson (the Peccei- 
Quinn axion) would appear in the CP-odd scalar sector, as can be verified by computing 
the determinant of A^p in the basis {A, Sj) in (2.27) for k = 0. 

The decomposition of the Peccei-Quinn axion ApQ in terms of the weak eigenstates H^j, 
Hdi and Sj is given by 

ApQ = — (w sin 2/3 A — 2 s Sj) , where 

N = Y^t;2sin2 2/3 + 4s2 , A = cos P Hui + sin P Hdi . (2.44) 

Hence, in most of the parameter space where s ^ v sin 2/3, ApQ is dominantly (but never 
purely) singlet-like. 

Apart from the strong cosmological constraints on a Peccei-Quinn axion, it is not 
straightforward to stabilise the potential for the vev s of the NMSSM in the Peccei-Quinn 
limit; obviously the approximate expressions in (2.19) and (2.21) are no longer appropriate 
for K = 0. In [61], couplings to additional singlets have been introduced in order to stabilise 
the potential for s. In [62-64], the limit A -C 1 has been considered, the potential including 
dominant radiative corrections in [65], and the Peccei-Quinn limit for large A in [66]. In 
any case, the stability of the scalar potential imposes strong constraints on the parameters 
of the NMSSM in the strict Peccei-Quinn limit k = [51,62-66]. 

If one prefers to avoid the cosmological constraints on a very light Peccei-Quinn axion, 
one can consider the situation where the Peccei-Quinn symmetry is explicitly broken by 
small additional terms in the superpotential of the NMSSM such as a small non-vanishing 
value for k in the Zs-invariant NMSSM [51] as it can be obtained in constructions of 
the NMSSM from the heterotic string [67]. Then the axion acquires a mass at tree level 
(without QCD contributions); if it is still light, it can lead to distinctive signatures at 
colliders due to the possible decay of the SM-like CP-even Higgs scalar into two pseudo- 
Goldstone bosons [68]. (For a scenario with an approximate Peccei-Quinn symmetry and 
large A, see [66].) 

In any case, the couplings of ApQ are important for its phenomenological signatures: 
the couplings of ApQ to gauge bosons, quarks and leptons are induced by its doublet 
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components indicated in (2.44), which are small for s ^ f sin 2/3. As a result, its couplings 
9 Add to down-type quarks and leptons are not enhanced for large tan/3, but given by 

gAdd-^gZd (2.45) 

where g^^i is the coupling of the SM-like CP-even Higgs scalar to down-type quarks and 
leptons. The couplings of A-pq to all CP-even scalars Hi vanish in the strict Peccei-Quinn 
limit, but become non-zero as soon as the Peccei-Quinn symmetry is explicitly broken. 
Then, while small, these couplings can still induce a branching fraction for Hi ^ 2 A-pq 
which is larger than the branching fraction for Hi — )■ hh [68]. 



2.2.3 The R-symmetry limit 

For A\, Af, — )■ 0, the Higgs sector of the Zs-invariant NMSSM (specified by the superpo- 
tential (2.6) and the soft terms in (2.8)) is invariant under an R-symmetry under which the 
scalar fields transform as 

H^-^H^e"^^ H^^H^e'^^, S ^ S e'^^ . (2.46) 

Again, this global symmetry is spontaneously broken by the vevs Vd and s, and a 
massless Nambu-Goldstone boson (now an R-axion) would appear in the CP-odd scalar 
sector, as can be verified by computing the determinant of A^p in (2.25) for i?efr = i^s, 
Ax = Ak = 0. The decomposition of the R-axion A^i in terms of the weak eigenstates H^i, 
Hdi and Si is given by 



Aji = ^{v sin 2l3A+sSi), where N = ^ sin^ 2/3 + (2.47) 

and A as in (2.44). Again, in most of the parameter space where s ^ v sin 2/3, Aji is 
dominantly (but never purely) singlet-like. However, in the full Lagrangian this R-symmetry 
is explicitly broken by the gaugino mass terms, which induce non-vanishing values for Ax 
(and subsequently for Af^) through radiative corrections, as can be seen from the /3-functions 
for Ax and in Appendix B.3. Even if the statement Ax = A^ = is scale dependent 
(and hence unnatural), one can still consider the situation where both trilinear couplings are 
relatively small [69-71] implying a light "pseudo" R-Goldstone boson given approximately 
by Ar as in (2.47). 

The couplings of Ap to down-type quarks and leptons (in the R-symmetry limit and for 
large tan/3) differ by a factor 2 from (2.45) and are given by 

9Add-^-^9Tdd, (2.48) 

while the couplings to CP-even scalars Hi depend on the values of Ax and A^. Again, 
this scenario can allow for possible (even dominant) iJj — )■ 2 Aji decays with important 
consequences for Higgs searches at colliders [69-71]. 

Hence, light pseudoscalars can easily appear in the NMSSM in the form of (pseudo-) 
Nambu-Goldstone bosons; however, light pseudoscalars in the NMSSM can also result from 
accidential relations among the parameters in which case they would not suffer from sup- 
pressed couphngs to down-type quarks and leptons. Phenomenological consequences of 
light pseudoscalars will be discussed in Sections 5 and 6. 
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3 Radiative corrections 



3.1 Renormalisation group equations 

Supersymmetry allows to formulate models at a very high scale, such as the GUT or Planck 
scale, while avoiding quadratically divergent quantum corrections involving an ultraviolet 
cutoff A ~ Mgut or a ~ Mpianck- However, assuming Msusy ~ ^weak, the quantum 
corrections still generate large logarithms ~ In (A/Msusy)- 

Fortunately, these can be summed up by the introduction of scale dependent parameters 
in the Lagrangian, where the scale dependence is described by the renormalisation group 
equations (RGEs) or /3-functions. If the parameters are assumed to be given at a large scale 
A ~ Mgut or A ~ Mpianck? the large logarithms are accounted for by the integration of the 
RGEs down to a low scale Msusy ~ M^cak- 

In the case of the three gauge couplings of the SM, this procedure allows for a successful 
unification at Mqut ~ 2 x 10^^ GeV in the MSSM [6-9]. It should be underlined that 
gauge coupling unification remains valid in the NMSSM: the additional gauge singlet field 
has no effect on the one-loop /3-functions of the gauge couplings, and the additional Yukawa 
couplings A and k only appear in the two-loop terms of the gauge /3- functions (see (B.l) 
in Appendix B). The resulting effect on the numerical values of the gauge couplings at the 
GUT scale is negligibly small (of the order of the unknown threshold effects at the GUT 
scale) . 

However, this is true only if the running Yukawa couplings A, k, hf and hi, remain suffi- 
ciently small so that perturbation theory remains valid; typically one requires that Yukawa 
couplings remain below 1 or a/Itt. Once the running Yukawa couplings have assumed values 
of 0{1) at a given scale, naive extrapolations of the perturbative RGEs result in singulari- 
ties (so-called Landau singularities) for the running couplings very close to that scale; as a 
result, upper bounds of 1 or -y/ivr have nearly the same consequences in practice. Since the 
running Yukawa couplings increase towards large scales (unless the one-loop contributions 
from gauge couplings are dominant, see the /3-functions in Appendix B), avoiding Landau 
singularities implies upper bounds on the couphngs at the scale Msusy ~ M„eak- 

The /3-functions for the Yukawa couplings have been derived for the NMSSM to one- 
loop order in [53], and to two-loop order in [72,73], see Appendix B.l. (Studies of analytic 
solutions, quasi-fixed points and RG invariants have been performed in [74-78] ; the impact 
of the NMSSM specific Yukawa couplings on hb-hr unification was studied in [79].) Upper 
bounds on the Yukawa couplings at the scale M„eak are given in [46, 80-87] to one-loop 
order, and in [72] to two-loop order, with the following results: The top quark Yukawa 
coupling hf at the scale rrit depends not only on the top quark mass mj, but increases 
with decreasing tan /3 as given in (A. 11). Hence, for small tan/3, hf can become too large 
( > a/Itt) at the GUT scale. Depending on the precise value of rrit and on the threshold 
corrections between the scales rrit and Msusy, one obtains a lower bound on tan/3 close to 
unity (typically between 1.5 and 2) both in the MSSM and the NMSSM. A similar reasoning 
based on hf, leads to an upper bound on tan/3, tan/3 < 80. 

The NMSSM specific Yukawa coupling A plays an important role for the upper bound 
on the mass of the lightest CP-even Higgs scalar, see (2.23). The validity of perturbation 
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theory up to the GUT scale imphes A < 0.7 — 0.8 at the weak scale, depending on the precise 
value of TTif, on the threshold corrections, on k and notably on tan/3: a large value of 
(i. e. small values of tan /3) amplifies the increase of A towards large scales, thus implying a 
decreasing upper limit on A at low scales. Larger values of k, have the same effect on A. 

These features are visible in Fig. 1, where we plot Amax as function of tan (3 and n (using 
NMSSMTools [52,88], see Appendix D). The black (dark) and red (grey) bands correspond 
to different choices of the Higgs and sparticle spectrum: the black bands correspond to 
a light spectrum with sfermion masses of 200 GeV, /igff = 100 GeV, a bino mass Mi 
of 50 GeV (with M2, M3 from (3.1) below) and Higgs masses below ~ 200 GeV. The 
red bands correspond to a heavy spectrum with sfermion masses, fiea and heavy Higgs 
bosons of 1 TeV, and Mi = 200 GeV. (The squark/slepton trilinear couplings are chosen 
to vanish in both cases for simplicity.) Inside the black and red bands, the top quark 
mass is 171.2 ± 2.1 GeV [89], a larger top mass leading to a slightly lower upper bound 
Amax- However, the variation of Amax with the threshold corrections of the running gauge 
and Yukawa couplings induced by the unknown sparticle and (heavy) Higgs spectrum is 
numerically more important. 




Figure 1: Left panel: upper bound on A (Amax) as a function of tan /3 for fixed n = 0.01. 
Right panel: Amax as a function of k for fixed tan/3 = 10. Black (lower) bands: light 
spectrum, red (upper) bands: heavy spectrum. Inside the bands the top quark mass is 
171.2 ±2.1 GeV. 

On the left-hand side in Fig. 1, k is fixed to 0.01; one sees both the decrease of the 
upper bound on A for tan/3 ~ 1 and very large tan/3 induced by very large values of ht, 
hf,, respectively. On the right-hand side, tan /3 = 10 and one sees the decrease of the upper 
bound on A with increasing k. 

Motivated by simple models for supersymmetry breaking, the soft SUSY breaking pa- 
rameters are often assumed to satisfy particularly simple relations at the GUT or the Planck 
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scale such as universal gaugino masses M1/2, universal trilinear couplings Aq and universal 
scalar masses mg (where, sometimes, the singlet mass and trilinear coupling are allowed to 
play a special role, see Section 7.1.3). Then, the RG evolution of these parameters from 
the high to the weak scale plays an important role. 

The one-loop /3-functions for all soft SUSY breaking parameters of the MSSM were 
computed in [90,91], for the NMSSM in [53,72], and for a general supersymmetric model 
in [92,93]. The general two-loop /3-functions for the soft terms can be found in [94,95], 
from which those of the (general) NMSSM can be deduced. They are given in Appendix B. 

Universal gaugino masses Mi = M2 = M3 = M1/2 at the GUT scale imply specific 
relations among these parameters (defined in (2.30)) at a low scale. At 1 TeV one finds (for 
squark masses < 1 TeV) 

M3 : M2 : Ml ~ 5.5 : 1.9 : 1 (3.1) 

with somewhat larger ratios (depending on threshold effects) below 1 TeV. These values 
hold for the MSSM as well as for the NMSSM. 

In the case of universal scalar masses mo at the GUT scale (with > 0), the RG 
evolution down to the weak scale is essential for the Higgs masses: at the weak scale, at 
least one Higgs mass squared has to be negative in order to trigger electroweak symmetry 
breaking. Fortunately rnj^^ < at the weak scale is induced nearly automatically by the 
large top Yukawa coupling [19], given the large top quark mass. 

In models with GMSB, the soft SUSY breaking terms are induced radiatively at a 
messenger scale Mmcss; which is typically far above the weak scale. In principle, all possible 
terms in the Lagrangian of the general NMSSM can be generated (see [96] and Section 7.2). 
Then, all these parameters have to be evolved by the RGEs from the messenger scale down 
to the weak scale. For completeness, we give the 2-loop /3-functions for all parameters of 
the general NMSSM in Appendix B. 

3.2 Radiative corrections to the Higgs masses, and the upper 
bound on the hghtest CP-even Higgs mass 

The Higgs mass matrices in Section 2, and notably the upper bound (2.23) on the mass 
of the lightest CP-even Higgs scalar, have been derived from the tree level potential (2.9), 
where all quartic terms are determined by supersymmetry through the superpotential and 
(electroweak) supersymmetric gauge interactions. This hypothesis would be justified only 
if the scale of supersymmetry breaking Msusy would be smaller than the Higgs vevs (or 
Mz), which is obviously not the case. 

However, for Msusy > Mz the deviation of the Higgs potential from (2.9) is calculable 
- as a function of the mass splittings among the superpartners - in the form of quantum 
corrections involving scales (momenta) with M| < < ^susy- "^^^ dominant con- 
tributions to ^Higgs originate from top quark/squark loops, since these particles have the 
largest couplings to Higgs fields (to Hu), and lead to an increase of the upper bound (2.23) 
on the mass of the lightest CP-even Higgs scalar. This phenomenon had first been discussed 
in the MSSM to one-loop and subsequently to two-loop order in [97-103] where, at tree 
level, (2.23) with A = would imply a mass of the lightest CP-even Higgs scalar below Mz- 
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An approximate formula for the mass M^m of the SM-hke Higgs scalar in the NMSSM 
in the limit ks ^ |^k|, |^a| (corresponding to a heavy singlet-like scalar), including the 
dominant top/stop radiative corrections, is given by 

MIm - M|cos2 2/3 + AVsin2 2/3- — t;2(A-Ksin2/3)2 

where v is defined in (2.12), the soft SUSY breaking stop masses squared in (2.36) are 
assumed to satisfy m"^ ~ mg^ ^ mf, At is the stop trilinear coupling assumed to satisfy 
\At\ ^ rrit, /iefr; the terms ~ are specific to the NMSSM, and the last term in the first 
line originates from the mixing with the singlet-like scalar. In the MSSM, where A = 0, 
the LEP bound on M^m implies that tan /3 has to be large such that cos 2/3 ~ 1, above 
~ 300 GeV for maximal mixing {Af ~ Qm^, maximising the second line in (3.2)), or 
> 1 TeV otherwise. 

In order to maximise Msm in the NMSSM, A should be as large as possible, and tan /3 
should be small in order to avoid a suppression from sin^ 2/3. (As discussed before, A is 
bounded from above by A < 0.7 — 0.8 if one requires the absence of a Landau singularity 
below the GUT scale.) However, the negative contribution from the mixing with the singlet- 
like scalar should vanish; without neglecting Ax, the relevant mixing term is proportional 
to (A — sin2/3(K + Ax/{2s)))'^ [104]. If this expression is not small, a larger value of A can 
even generate a decrease of the mass of the Higgs scalar with SM-like couplings to the Z 
boson in the NMSSM. 

The resulting upper bound on the lightest CP-even Higgs mass in the NMSSM has been 
studied in the leading log approximation in [82-84,103,105-109]. Full one- loop calculations 
of the corresponding upper bound involving top/bottom quark/squark loops have been 
carried out in [72,85,86,110-117]. (Analyses at large values of tan /3 have been performed 
in [118-120], and upper bounds for more general supersymmetric Higgs sectors have been 
considered in [121-123].) 

At present, additional known radiative corrections to the Higgs mass matrices in the 
NMSSM include MSSM-like electroweak together with the NMSSM-specific Higgs one-loop 
contributions [124,125] and dominant two-loop terms [73,87,125-127]. In order to discuss 
these in detail, it is convenient to separate the quantum corrections involving scales 
with Q'^ > M|ugY from those with scales < M^^^y. 

The result of the quantum corrections with > M^^^y is still a supersymmetric 
effective Lagrangian (including soft SUSY breaking terms), where all running parameters 
(couplings and masses) are defined, within a given subtraction scheme, at the scale ~ 
Mlugy. (If desired, the parameters at the scale can be obtained in terms of parameters 
at a higher scale with the help of the RGEs.) Subsequently, the quantum corrections with 
Q2 < M|ugY (i- e. with an ultraviolet cutoff M^^^y) have to be evaluated, generating a non- 
supersymmetric effective action including an effective Higgs potential, effective couplings 
of fermions and wave function normalisation constants. From the effective potential and 
couplings one can derive the so-called running masses, which still differ somewhat from the 
physical pole masses (the poles of the propagators). 
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The terms in the effective action can be computed in a systematic expansion in powers 
of couphng constants and large logarithms where, to start with, the couplings are defined 
at the scale ~ M^^^y. Hence, some additional effort is necessary in order to relate the 
gauge couplings to those measured at the scale M|, the Higgs vevs to Mz, and the Yukawa 
couplings to the quark and lepton pole masses: the gauge and Yukawa couplings at the 
scale ~ ^susy have to be obtained from the measured couplings at lower scales by the 
integration of RGEs where, however, all possible threshold effects have to be taken into 
account. 

The electroweak and NMSSM specific Higgs one-loop contributions, and the two-loop 
contributions oc hf ag, have recently been computed in [125] without an expansion in large 
logarithms. However, subsequently we confine ourselves to electroweak and NMSSM-specific 
Higgs one-loop contributions involving large logarithms, and to two-loop contributions pro- 
portional to two powers of large logarithms. In Appendix C we summarise the corresponding 
formulae, which allow to determine the lightest CP-even Higgs mass in the NMSSM with 
the following accuracy (considering the tree level mass squared following from the first line 
in (3.2) to be of order g"^ or A^, where g stands for gi or (72): 

(i) One-loop contributions from top/bottom quarks/squarks are fully included, which gen- 
erate contributions of order h^, and g"^ /i^ etc., possibly multiplied by a large logarithm 
In (Mlugy/m^). (Here and below we identify Mgugy with an average value of the soft SUSY 
breaking squark mass terms, Mgygy = f^Qs friT- It is possible to allow for different defini- 
tions for M|ugY [52] in cases where other soft terms are larger than mq^ , m^, which we do 
not consider here for simplicity.) 

(ii) Additional one-loop corrections considered are of the orders (7'^, g^X^., g^n"^, and k"^, 
if multiplied by a large logarithm. Here, loops from squarks and sleptons, charginos and 
neutralinos, Higgs and gauge bosons provide relevant contributions. 

(iii) At two-loop order, only dominant logarithms ~ In^ {MQ^jgy/mf) multiplied by or 
a<j (and hf replaced by h^, possibly relevant at large tan /3) are considered, including 
some subdominant effects due to the evaluation of ht and at the scale MgusY- 

Taking these loop corrections into account, and requiring perturbative running Yukawa 
couplings ht, A and k below the GUT scale, the upper bound on the lightest CP-even 
Higgs mass has been studied in [104] as a function of tan /3 and for different values of 
rrit in the NMSSM, and compared to the MSSM with the result shown in Fig. 2. (In 
Fig. 2, the upper bound is denoted as mmax-) The squark mass terms (and hence Mgusy) 
have been chosen as 1 TeV; the upper bound would still increase slowly (logarithmically) 
with MsusY- In order to maximise the one- loop top/bottom (s) quark contributions to 
the lightest CP-even Higgs mass for these squark masses, the trilinear soft couplings are 
chosen as At = Af, = 2.5 TeV. The threshold effects depend somewhat on the gaugino 
masses, which are Mi = 150 GeV, M2 = 300 GeV and M3 = 1 TeV (roughly in agreement 
with (3.1)); the remaining parameters A, k, Ax, A,^ and fi^s of the Zs-invariant NMSSM 
have been chosen such that the upper bound is maximised, without encountering Landau 
singularities (which requires k, as small as possible) nor violating other constraints such as 
an unstable potential, which forbids k — > 0. 

The lower dashed lines in Fig. 2 refer to the MSSM, where the mass of the CP-odd 
scalar - which can be chosen as the other independent parameter in the Higgs sector 
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Figure 2: Upper bound on the lightest Higgs mass in the NMSSM as a function of tan /5 
for rrit = 178 GeV {Ma arbitrary: thick full line, Ma = 1 TeV: thick dotted line) and 
rrit = 171.4 GeV (thin full line: Ma arbitrary, thick dotted line: Ma = 1 TeV) and in the 
MSSM (with Ma = 1 TeV) for mt = 178 GeV (thick dashed line) and mt = 171 A GeV (thin 
dashed line). Squark and gluino masses are 1 TeV and At = Ah = 2.5 TeV. (From [104].) 

apart from tan /3 - is set to Ma = 1 TeV. The other parameters (and approximations) are 
the same as described above. In the MSSM, the increase of the upper bound with tan^ 
originates from the tree level term (the first term ~ cos^2/3 in (3.2)), according to which 
minax is maximised for large tan/3. Due to the one-loop top (s)quark contributions, the 
upper bound rrtmax increases with mt. Numerically, a variation Arrit of rrit implies nearly 
the same variation Ammax for large tan/3. 

In the NMSSM, the second term ~ sin^ 2/3 in the tree level expression (3.2) dominates 
the first one for sufficiently large A, and accordingly mmax is maximal for low values of 
tan/3. On the other hand, the absence of a Landau singularity for A below the GUT scale 
implies a decrease of the maximally allowed value of A at MgusY with increasing ht, i.e. 
with increasing rrit and decreasing tan/3, see Fig. 1. (At large tan/3, arbitrary variations 
of the NMSSM parameters A, k, A\, and /Xeff can imply a mass Ma of the MSSM-like 
CP-odd scalar far above 1 TeV. For comparison with the MSSM, mmax in the NMSSM with 
Ma < 1 TeV is depicted as dotted lines in Fig. 2.) 

Numerically, for rut = 171.4 GeV, mmax in the NMSSM (thin full line in Fig. 2) is 
140 GeV, and assumed for tan/3 = 2, A = 0.70 and k = 0.05, whereas for rrit = 178 GeV 
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(thick full line in Fig. 2), mmax increases just to 141.5 GeV for tan/3 = 2.2, A = 0.68 and 
K = 0.07. 

It must be emphasised, however, that the non-observation of a Higgs boson in standard 
search channels with a mass below ~ 140 — 145 GeV (allowing for larger Msusy and 
remaining theoretical uncertainties) would not rule out the NMSSM: here the lightest CP- 
even Higgs scalar could have small couplings to all quarks, leptons and gauge bosons if it is 
dominantly singlet-like, in which case the mass of the Next-to-lightest Higgs scalar with SM- 
like couplings to gauge bosons is bounded from above, and can be ~ 20 GeV larger [127], see 
also Section 5.1. Also, the lightest CP-even Higgs scalar could have unconventional decay 
channels which can make its discovery at colliders quite difficult. In such a case one could 
possibly detect one of the heavier CP-even Higgs scalars of the NMSSM, whose masses and 
couplings can vary over a wide range. (Radiative corrections to their masses are also given 
in the Appendix C.) 

Finally, mmax in the NMSSM obviously increases if one allows for larger values of A 
[128,129]. Since this would imply a Landau singularity below Mgut for the particle content 
of the NMSSM, one can consider modifications of the theory at higher scales. For instance, 
one can allow for new strong gauge interactions at intermediate scales [130-132], large extra 
dimensions implying a lower effective ultraviolet cutoff [133], or replace the singlet S with 
its coupling A by an unparticle operator [134] (in which case the Higgs content would not 
be NMSSM-like). Modifications of the /3-functions due to additional SU{5) multiplets with 
masses in the TeV range have been considered in [66]: additional SU{5) multiplets lead to 
larger gauge couplings at larger scales which, in turn, affect the /3-function for A such that 
larger values of A at the weak scale are allowed. In [135] an extra SU (2) gauge symmetry has 
been introduced, acting on the third generation and the Higgs doublets only and allowing 
for A ~ 1 at low energy, still consistent with perturbativity and GUT-scale gauge coupling 
unification. This, together with small values for tan/3 (possibly below 1) leads to an upper 
bound on the lightest Higgs mass of mmax ^ 250 GeV. The phenomenology of large-A 
scenarios with A up to 1.5 has been studied in [136], again with the result that tan /3 must 
be small for A — > 1.5. 

3.3 Radiative corrections to coupling constants 

Radiative corrections to coupling constants can be ultraviolet divergent or (ultraviolet) 
finite. In SUSY extensions of the Standard Model, all ultraviolet divergent radiative cor- 
rections can be absorbed into redefinitions of the SUSY preserving or soft SUSY breaking 
masses and couplings, which are effectively described by the RG equations for these param- 
eters (valid between Msusy and, e.g., Mgut)- 

However, ultraviolet finite quantum corrections involving scales between M^eak and 
^susY can induce new couplings (or modify SUSY relations between existent couplings) 
which cannot be described in terms of an effective supersymmetric theory plus soft terms. 
Amongst others, such ultraviolet finite quantum corrections to the Higgs mass terms are 
responsible for the increase of the upper limit on the lightest CP-even Higgs mass with 
the logarithm of the top squark masses, as discussed in the previous subsection and in 
Appendix C. 

Subsequently we briefly discuss quantum corrections to hb and comment on the coupling 
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Zbb. Radiative corrections to the Higgs self-couplings are given at the end of Appendix A. 
We leave aside the threshold effects for the running gauge and Yukawa couplings at the 
various sparticle and Higgs mass scales; for a recent NMSSM specific analysis see [137]. 

Large tan 13 corrections to 

The superpotentials of the MSSM and the NMSSM contain couplings of to up-type 
quarks, and of to down-type quarks and leptons. It is well known [138-140] that squark- 
gluino- and squark-chargino loops can induce a coupling of to down-type quarks and 
leptons giving rise to an effective Lagrangian [141] 

£efr = hH^bb + AhHlbb+ ... . (3.3) 

As a result, the physical 6-quark mass mb is given by 

mf, = hvd + AhbVu = hhVd{l + Anib) (3.4) 

with 

Amfe = — — ^ tan /3 . (3.5) 

h 

Notably for large tan /3, large values of /i (or ^^s in the NMSSM), large gluino masses and/or 
large values of \At\, the quantity Ami, can be sizable; explicit expressions for Am^ in the 
MSSM can be found in [138-141]. These hold also in the NMSSM, provided fi is replaced 
by fies- As a first consequence, the formula for the Yukawa coupling hf, as function of 
in (A. 11) has to be corrected and should read 

' t;cos/3(l + Am,) ■ ^ ^ 

The corresponding value for hf, has to be used as boundary condition at the weak scale for 
the integration of the RGEs up to the GUT scale. Moreover, the couplings of the Higgs 
bosons to 6-quarks in (A. 12) have to be corrected as follows: 

^corr 

SibLb% : {Sii + Ambtanf3 Si2) 

v2 

ijCorr 

AibLb% : ^^Pii 

H-tLb% : -hrsm(3, (3.7) 

where the couplings of the CP-odd and charged Higgs bosons follow from the SU{2)- 
invariant completion of the effective Lagrangian (3.3). 

In the NMSSM, squark-gluino- and squark-chargino-loops also induce direct couplings 
of the singlet 5" to quarks and leptons [142]. They give rise to additional terms ~ A^hh Sbb 
in the effective Lagrangian (3.3), which have been worked out in the limit of large tan/3 
in [142], together with the additional corrections to the Higgs 6-quark couplings (3.7). The 
terms A^hb Sbb are proportional to the singlet-components of the Higgs mass eigenstates 
and to Vu,d/s; hence they are small for s ^ Vu,d- 
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The coupling Zbb 

In view of the discrepancy between the Standard Model and the forward-backward 
asymmetry A^p^ in 6-quark production at SLC and LEP [143], radiative corrections to 
the couphng Zbb are of particular interest. However, given that the measurements of 
Rb = T{Z — bb)/T{Z hadrons) agree quite well with the Standard Model, explanations 
of the discrepancies only in the asymmetries are not easy. A recent analysis of radiative 
corrections to the Zbb vertex in supersymmetric extensions of the SM including the NMSSM 
is given in [144], including possible NMSSM specific contributions to other electroweak 
precision observables. The net result of this study is, however, that the NMSSM cannot 
improve the agreement with the measurements, since NMSSM specific contributions have 
the opposite sign of the one required. 
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4 The tadpole and domain wall problems, baryoge- 
nesis, hybrid inflation 



4.1 The tadpole problem 

In order to avoid a fine tuning problem in the NMSSM, none of the dimensionful terms in 
the superpotential W^Higgs (2.1) or in Csoft in (2.5) should be much larger than the weak or 
SUSY breaking scale. 

However, if a singlet superfield 5" couples in the most general way to heavy fields (as 
it is possible in GUTs or GMSB), radiative corrections can induce very large terms in the 
effective action which are linear in S in the superpotential (denoted by in (2.1)) or 
linear in S in £soft (denoted by C,s in (2.5)). These terms are called tadpole terms, and if 
they are too large, the corresponding model has a "tadpole problem". This problem has 
been discussed in the context of globally supersymmetric theories with soft supersymmetry 
breaking terms - as induced by a hidden sector in supergravity - in [29,145-149], using 
general power counting rules and/or by explicitly evaluating Feynman diagrams. 

Once non-renormalisable supergravity interactions suppressed by powers of Mpianck are 
taken into account, power counting rules signal again a potential tadpole problem [150], as 
confirmed by explicit calculations in [151-153]. 

In the case of soft supersymmetry breaking terms ~ Msusy induced by a hidden sector 
in supergravity, the orders of magnitude of the radiatively induced tadpole terms are 

~ A Msusy , 

~ AMsV, (4.1) 

where A ~ Mgut in the presence of general couplings of S to GUT fields, or A ~ Mpianck 
in the presence of the most general non-renormalisable supergravity interactions of S. In 
both cases, the induced orders of magnitude of and are too large. 

In the case of GMSB (see Section 7.2), the phenomenologically required soft super- 
symmetry breaking mass terms of sparticles are induced radiatively; the source of super- 
symmetry breaking are so-called messenger fields (f with SM gauge quantum numbers and 
supersymmetric mass terms ~ Mmess, but whose real and imaginary scalar components have 
different masses, split by a scale m. This kind of supersymmetry breaking is also denoted as 
-F-type splitting, since it can be represented in terms of a non-vanishing F component of a 
spurion superfield which couples to the messenger fields (p. Then, Mgusy in (4.1) is given by 
Msusy ~ fn'^/M^acss, and A by Mmess- Hence, if Mmess as well as the F-type splitting m are 
not much larger than the weak scale, the singlet tadpole problem is absent [33,34,154-156]; 
if these scales are larger than the weak scale, the tadpole diagrams can be suppressed by 
small Yukawa couplings, see [96] and Section 7.2. 

In order to circumvent the tadpole problem in the NMSSM in the case of soft super- 
symmetry breaking terms induced by a hidden sector in supergravity, we note that the 
dangerous terms (4.1) with A ~ Mgut, Mpianck can be generated only if is a singlet 
with respect to all continuous and discrete symmetries of the full theory. As soon as the 
heavy/hidden sector is invariant under a discrete symmetry under which S transforms, the 



25 



terms (4.1) are absent if the discrete symmetry is unbroken, or A is of the order of the scale 
of the breakdown of the discrete symmetry (possibly multiplied by high powers of coupling 
constants and loop factors l/167r^). Hence, if the NMSSM is embedded in a GUT or su- 
pergravity theory with an exact or approximate discrete symmetry such that A < Msusy, 
the tadpole problem is avoided. 

4.2 The domain wall problem 

Discrete symmetries can generate another problem, however [157]: once they are spon- 
taneously broken after a symmetric phase in the hot early universe, domain walls are 
generated which can dominate the energy density of the universe, creating unacceptably 
large anisotropies of the cosmic microwave background radiation and spoiling successful 
nucleosynthesis [157]. 

In particular, the NMSSM with a scale invariant superpotential is invariant under a Z3- 
symmetry (see the discussion following (2.7)). Then, after electroweak symmetry breaking, 
the universe would consist of different "bubbles" of the same vacuum energy, but in each of 
which the phases of {H^) and {S) would differ by a Za-transformation; these bubbles 

are separated by domain walls [158-162]. 

However, since the Zs-symmetry is just an accidental symmetry of any scale invariant 
superpotential, it is not expected that it will be preserved by Planck scale suppressed (non- 
renormalisable) gravitational interactions. Whereas such a violation of the Zs-symmetry 
is typically sufficient to avoid the domain wall problem (since the vacuum energy within 
different bubbles would be slightly different, leading to a collapse of the ones with higher 
energy) [158-161], it may again lead to dangerously large tadpole diagrams [150-153,162]. 
It has been believed [162], that the conflict between the domain wall and tadpole problems 
cannot be solved. 

Subsequently solutions of this problem have been proposed in [163-165]: it is possible to 
impose constraints on Za-symmetry breaking non-renormalisable interactions or (renormal- 
isable) hidden sectors in the form of various additional symmetries, such that Zs-symmetry 
breaking renormalisable terms - as the tadpole terms above - are generated radiatively, but 
with very small coefficients. These Za-symmetry breaking terms can still solve the domain 
wall problem of the otherwise Za-invariant NMSSM, without having a visible impact on its 
phenomenology. 

In [166] a Zs-R-symmetry of the non-renormalisable interactions has been considered, 
as a consequence of which a tadpole term ~ is generated only at the six-loop level. 
If the order of magnitude of is ^ ^susy^ i^ constitutes one of the possible terms 
that can solve the domain wall problem of the otherwise Zs-invariant NMSSM. However, 
one can also assume that C,s ~ ^susy [166]; in this case the trilinear term ~ ^S^ in the 
superpotential of the NMSSM is not even phenomenologically required and can be omitted. 
The resulting model has been denoted as the New Minimal MSSM (or nMSSM [167], if a 
term ~ Mgugy is added to the superpotential as well), which we will discuss in more 
detail in Section 7.3. 
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4.3 Electroweak baryogenesis 

A strongly first order electroweak phase transition is a necessary (albeit not sufficient) 
condition for baryogenesis relying on anomalous baryon number violating Standard Model 
processes [168]: during a first order phase transition bubbles are formed, within which the 
Higgs fields assume non-vanishing values corresponding to the absolute minimum of the 
zero temperature effective potential. These bubbles are separated by expanding "walls" 
from the phase with vanishing Higgs vevs corresponding to the minimum of the effective 
potential at high temperature. If the dynamics of the processes across the expanding walls 
sufficiently violates CP and baryon number, the baryon asymmetry of the present universe 
could be explained [168]. 

A first order electroweak phase transition is difficult to achieve in the SM and in simple 
supersymmetric extensions as the MSSM, given the present lower bounds on Higgs and 
stop masses (which play an important role for the one-loop corrections to the effective 
potential); see, e.g., [169]. On the other hand, the more complicated tree level Higgs 
potential in various versions of the NMSSM can easily give rise to a sufficiently strong 
electroweak phase transition [170-179]. However, the trajectory in the space of fields S, 
and must be studied carefully in order to avoid a phase transition in two steps where 
first S alone assumes a vev (implying an insufficient violation of CP) and subsequently 

and assume their vevs as in the MSSM with fixed fi^s- To this end, the vevs of 
5", and should be of the same order along the trajectory in field space [172,178], 
without violating present bounds on the Higgs sector. This seems easier to achieve in the 
presence of additional terms in the superpotential and/or £soft beyond those of the Z3- 
invariant NMSSM [171,173,174]. Simultaneously, such additional terms explicitly break 
the Zs-symmetry avoiding the domain wall problem. The simplest additional terms are 
actually the tadpole terms ~ and/or ~ ^p, which allow to omit the terms ~ k of the 
Z3-invariant NMSSM; the resulting model is the nMSSM mentioned above and discussed 
in Section 7.3. Electroweak baryogenesis in the nMSSM has been studied in [176,177,179]. 

In addition to a first order phase transition responsible for the formation of bubble walls, 
the processes inside the expanding bubble walls have to violate CP in order to generate 
a baryon asymmetry. CP violation can originate from the chargino sector (through a 
non-vanishing phase of M2 x fi^s in the NMSSM) or the Higgs sector (the prospects for 
CP violation in the Higgs sector of the NMSSM are reviewed in Section 8.1). In both cases, 
strong constraints from the upper bound on the (neutron) electric dipole moments have to 
be respected, which is possible in the nMSSM [176,179]. 

The possibility of spontaneous CP violation inside the bubble walls without CP viola- 
tion at zero temperature has been advocated in [180,181]; then, bounds on the neutron 
electric dipole moment are trivially satisfied. A similar mechanism inside the domain walls 
separating phases related by the Zs-symmetry of the Zs-invariant NMSSM has been con- 
sidered in [161]. Electroweak baryogenesis in a f/(l)' model with a secluded ?7(l)'-breaking 
sector (see Section 8.2.1), where one can have a sufficiently strong first order electroweak 
phase transition, has been studied in [182,183]. Once the criteria for successful baryogenesis 
are fulfilled in the NMSSM, the background density of gravitational waves produced during 
the electroweak phase transition can be within reach of the gravitational wave experiment 
LISA [184,185]. 
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4.4 Hybrid inflation 

A variant of the NMSSM without domain wall problem, solving the strong CP problem 
via an invisible axion, and allowing for hybrid inflation has been proposed and discussed in 
[186-191]. Two singlets (f) and are introduced, and the relevant part of the superpotential 
is given by 

W = \<pHuHd + K(PN^ . (4.2) 

The soft SUSY breaking mass squared of (p is assumed to be very small and negative, and 
both Yukawa couplings A and k are assumed to be A ~ /t ~ 10"^*^ which can be motivated by 
brane constructions and/or type I string theory. The model has a Peccei-Quinn symmetry 
(see Section 2.2.2) with Peccei-Quinn charges -2 for (p and +1 for A^. During the inflationary 
epoch the vevs of A^, Hu and vanish, and the infiaton rolls slowly along an almost 
flat direction of the scalar potential. Only near the end of inflation A^, Hu and Hd develop 
vevs. Domain walls, created by the spontaneous breaking of the Peccei-Quinn symmetry, 
are diluted by the inflation. The vev of after inflation is automatically of the correct 
order such that fics = A0 is of the order of the SUSY breaking scale, and the constraints on 
the parameters associated with inflation (number of e-folds, curvature perturbations and 
the spectral index) can be satisfied. 
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5 NMSSM phenomenology at colliders 



Due to the additional singlet superfield S, the phenomenology of the NMSSM can differ 
strongly from the MSSM in the Higgs and neutralino sectors. Subsequently we discuss these 
subjects separately. 



5.1 The Higgs sector 

The mass matrices of the three CP-even and the two (physical) CP-odd neutral Higgs 
sectors are given in Section 2.1 at tree level, and the dominant radiative corrections in 
Appendix C. In general, after the diagonalization of the mass matrices as in Appendix A, 
the SU{2) doublets mix with the singlet states in both sectors. As a consequence, the 
reduced couplings 

= sin /3 Si2 + cos /3 Sn (5.1) 

of the 3 CP-even mass eigenstates Hi to the electroweak gauge bosons (normalised with 
respect to the couplings of the SM Higgs scalar) can be very small; however, they always 
satisfy the sum rule 



E^' = l- (5-2) 



1=1 

Another useful sum rule involving the eigenvalues of the CP-even mass matrix 
has been given in [167]: 

3 

= cos' ^ {Ml)^^ + 2cosf3sm(3 {Ml)^^ + sm' ^ {Ml)^^ 

1=1 

= Ml l^cos^ ^'^ + ^ ^^^^ ^'^) + ' '^^•^^ 

where the last line is the upper bound rn^^^ on the lightest CP-even Higgs mass in the 
NMSSM (corresponding to (2.23) at tree level) including the radiative corrections which 
lift it up to ~ 140 GeV, see Fig. 2. Using Mh^ > (by definition) and eliminating 

^2 + ^3 with the help of (5.2), one can derive from (5.3) [192] 

< (mL. - eiMl) , (5.4) 

which will be useful below. 

The extended CP-even and CP-odd Higgs sectors allow for the possibility of additional 
Higgs-to-Higgs decays compared to the MSSM. Subsequently we will discuss separately the 
regions in the parameter space of the NMSSM where Higgs-to-Higgs are allowed/disallowed. 



5.1.1 LEP and e"*" e colliders 

To start with, the constraints from LEPl and LEP2 must be translated carefully into 
constraints on the parameter space of the NMSSM. Early studies [113,127,193-196] confined 
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themselves to the standard Higgs search channels 

e+e^ ^HZ, H ^ bb, t+t~ , (5.5) 

neglecting possible unconventional Higgs decay modes. 

Since the lightest CP-even scalar Hi can have a dominant singlet component, its reduced 
coupling ,^1 to Z bosons can be smaller than in the MSSM with, possibly, .^i — )■ 0. The 
LEP Working group for Higgs Boson Searches [57] has published upper bounds on C,"^ as 
function of Mh (valid for any Hi) combining results from the ALEPH, DELPHI, L3 and 
OPAL collaborations, as shown in Fig. 3. Clearly, even very light CP-even scalars are not 
ruled out if is sufficiently small, and such scenarios can indeed be realised within the 
NMSSM [113,114,127,175,192-199]. 




mjj(GeV/c^) 

Figure 3: Upper bounds on C,"^ from LEP [57], where SM branching ratios H ^ bb and 
H ^ T~ are assumed. Full line: observed limit; dashed line: expected limit; dark 
(green) shaded band: within 68% probability; light (yellow) band: within 95% probability. 

In the limit .^i — )■ corresponding to a light singlet-like CP-even Higgs scalar, (5.4) 
shows that the upper bound on the Higgs mass discussed in Section 3.2 applies now to 
the Next-to-lightest CP-even Higgs scalar iJ2, which helps to establish no- lose theorems at 
colliders (stating that at least one Higgs boson is detectable). In the case 7^ 0, but ^1 
small, the lightest CP-even scalar Hi would still be very difficult to detect at the Tevatron 
and at the LHC, and Higgs searches will only be sensitive to H2 (and possibly H^) in 
the CP-even sector. For ^1 7^ 0, but Q < the upper bound on the mass of the 
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observable H2 is somewhat alleviated [127,192]. Using again (5.4) and the upper bound 
on ^1 as a function of Mh^ from Fig. 3 (assuming SM decay branching ratios for Hi), one 
finds that the maximal possible value of Mh^ is assumed for ~ ^| ~ 1/2, in which case 
Fig. 3 gives Mh^ > 110 GeV implying Mh, < 162 GeV. 

The relevance of Higgs-to-Higgs decays in the NMSSM for Higgs searches was first 
mentioned in [195], and underlined in [68,69] in the framework of scenarios with an ap- 
proximate Peccei-Quinn or R-symmetry (see Section 2.2) where a CP-odd scalar Ai can 
be very light allowing for Hi — )■ Ai Ai decays. Then the branching ratios for H^ — )■ bb 
and Hi — )■ r"*" r~ would be suppressed, and the upper bounds on in Fig. 3 should be 
re- interpreted [52,200-202]. In particular, a CP-even scalar H with a SM coupling to the 
Z boson [C, = 1) but with Mh < 114 GeV can be compatible with the constraints from 
LEP [70,71,200-207]. 

On the other hand, once Higgs-to-Higgs decays as Hi — )■ Ai Ai (possibly also Hi — )■ 
Hi Hi) are kinematically allowed, many additional Higgs search channels studied at LEP 
can be relevant. The constraints on the corresponding cross sections, after combining results 
from the ALEPH, DELPHI, L3 and OPAL collaborations, have also been worked out by 
the LEP Working group for Higgs Boson Searches [57], and those particularly relevant for 
the NMSSM have been summarised recently in [208]: 

1. Searches for e~ H Z independent of the H decay mode, looking for a peak of the 
Mx recoil mass distribution in e+ e" X Z, by OPAL [209]; these give Mh > 82 GeV 
if H couples to Z with SM strength (,^ ~ 1). 

2. Searches for if -> $ $ and $ 66 ($ = CP-even or CP-odd Higgs) by OPAL [210] 
and DELPHI [211]; once combined, these give Mh > 110 GeV for ^ ~ 1 [57]. 

3. Searches for if $ $ and $ ^f^f, $ cc, $ r+ r~ by OPAL [212] exclude 
45 GeV < < 86 GeV for ^ ~ 1. A new analysis of ALEPH data [213] implies 
stronger limits on if -> 4ras Mh > 109 GeV for M$ = 10 GeV and 
^2 X BR{H ^ $ $) X BR{^ rrf ~ 1. 

Also possibly relevant could be [57] e" H Z with H jj (two jets), if — )■ 77 
and H decaying invisibly (as, e.g., into neutralinos) , and e" — > HA together with 
H A 3A 6b, all of which are considered in the pubhc code NMHDECAY [52,88], 
which is part of NMSSMTools (see Appendix D). 

In Fig. 3 one can note a slight excess of events for Mh ~ 95 — 100 GeV of 2.3 cr statistical 
significance [214]. In the NMSSM this excess could be explained by 

• a SM-like CP-even scalar of corresponding mass and ^ ~ 1, but with reduced BR{H — )■ 
bb) ^ 0.1 [70,202,205] due to a dominant decay H ^ AA with a branching ratio 
~ 0.8. This scenario could alleviate the little fine tuning problem related to the non- 
observation of a CP-even scalar at LEP (see Section 7.4), but requires Ma < 10.5 GeV 
due to the constraints (2) on if — )■ $ $ and ^ bb; 

• a CP-even scalar ifi of corresponding mass with ^ < 0.4 [215] as it can appear, e.g., 
in the constrained NMSSM (see Section 7.1) [216,217] or for large A ^ 0.7-0.8 [218]. 
Then, another CP-even scalar if2 with C, > 0.9 and a mass not far above 115 GeV is 
present in the spectrum. 
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In the first scenario, difficulties for Higgs searclies at tlie Tevatron and at tlie LHC are 
foreseeable (see below), whereas probably only H2 would be detectable at hadron colliders 
within the second scenario. 

If Hi is singlet-like {^l < and decays unconventionally, the upper bound on 
deduced above {Mh2 ^ 162 GeV from Mh^ > 110 GeV) is no longer valid. Using the upper 
bound on Mhj^ from constraints on Hi Ai Ai ^ 46, (5.4) gives ^ 167 GeV, whereas 
the weaker upper bound on M^^ independent of its decay mode (using only the Z recoil 
mass) allows for ^ 192 GeV. 

The prospects for NMSSM Higgs searches at e"*" e~ colliders of higher energy, as an ILC, 
have been investigated in [177,192,194,219-222]; for a review see [223]. Due to the sum 
rule (5.2), one can derive a theoretical lower limit on the production cross section cTj for 
at least one CP-even Higgs boson, obtained as (Xj > 0.04 pb for y/s = 300 GeV in [192]. 
(Production cross sections at a/s = 500, 1000 and 2000 GeV were computed in [194].) 

As at LEP, one can search for Higgs bosons at an ILC independently from their decay 
modes using the recoil spectrum (the invariant mass Mx) in e~ Z X , Z ^ e~ 
and Z — )■ /i"*" This allows to establish a no-lose theorem for arbitrary singlet extensions 
of the MSSM (assuming perturbativity of all Yukawa couplings up to 10^^ GeV, implying 
upper bounds on the Higgs masses) for ^/s = 500 GeV and an integrated luminosity C = 
100 fb~^ [219]. Also the ZZH couplings can be determined by this method [223]. The clean 
environment of e~ colliders would allow for the study of NMSSM specific Higgs-to-Higgs 
decays [220], and for the discovery of light Higgs states that would be difficult to observe 
at the LHC [221,222]. 

5.1.2 Tevatron and LHC 

Next we turn to hadron colliders, starting with regions in the NMSSM parameter space 
where Higgs-to-Higgs decays are kinematically forbidden or occur only with small branching 
ratios. Then, the search for Higgs bosons in the NMSSM can proceed as in the MSSM 
(see [224] and references therein); of course, Higgs boson production cross sections and 
decay branching fractions in the NMSSM have to be rescaled (for tree level couplings) 
or re-evaluated (for loop corrected or loop induced couplings), employing the tree level 
couplings given in Appendix A. 2. Important search channels are 

1. gg ^ H (gluon-gluon fusion) with H W'^*'> l^Tuu, H ^ Z Z'-*^ 4 leptons, 
H ^ '-ff (the latter also for A instead of H); 

2. Vector boson fusion corresponding to qq^qqVV with VV ^ H at the Tevatron, 
and qq -> qqVV with VV H at the LHC (where V = W^, Z), and -> 77, 
H ^ T~ , H ^ WW^*\ H decaying invisibly 01 H ^ bb (the latter decay suffers 
from a too large background at the LHC); 

3. qq ^ W* ^ W H and qq Z* ^ Z H (Higgs-strahlung); 

4. Associated production with heavy quark pairs qq/gg — )■ QQ H, where Q = t, b. 

In some regions of the parameter space of the NMSSM, more Higgs bosons than in 
the MSSM could be observed (e.g. in the two photon channel [225]), which would allow 
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to distinguish the two models. On the other hand, since all can be simultaneously less 
than 1 (respecting the sum rule (5.2)), it was not clear at first instance whether any of 
the Hi could always be detected after combining the results for Higgs searches at LEP2 
with the prospects for the LHC [195,198], even if Higgs-to-Higgs decays are assumed to be 
kinematically forbidden. Later on, the analysis of Higgs searches at the LHC was refined and 
additional Higgs production channels (notably vector boson fusion) have been considered. 
On the theoretical side, radiative corrections to the Higgs masses have been improved by 
dominant two-loop corrections; these lower the CP-even Higgs masses, increasing the region 
in the NMSSM parameter space excluded by LEP2. As a result, a no-lose theorem for the 
detection of at least one CP-even Higgs scalar in the NMSSM at the LHC with an integrated 
luminosity of 600 fb^^ (two detectors with 300 fb~^ each) could be established [201,220,226], 
however, under the assumption that Higgs-to-Higgs decays are kinematically forbidden. 

Once Higgs-to-Higgs decays are possible (or even dominant), Higgs searches at the 
Tevatron and the LHC can become considerably more complicated; see [66,201,203,208,227- 
229] for discussions of possible scenarios and proposals for search channels. In principle, CP- 
even scalars H2 can decay into a pair of CP-even Hi in LEP-allowed regions of the NMSSM 
parameter space [201,217,228] (see [129] for very large A ^ 2 where M^i ^ 200-300 GeV, 
~ 350 — 700 GeV); but most of the recent studies concentrated on decays into light 
CP-odd scalars A motivated by an approximate Peccei-Quinn or R-symmetry, and/or the 
H — > AA explanation of the light excess of events at LEP. 

For Ma > 2mb (in which case LEP constraints require Mh ^110 GeV for SM-like 
couplings of H), the cascade H^AA would end mostly in a 4 6 final state. This final state 
has a very large background at hadron colliders (which seemed to make it invisible at least 
at the Tevatron, unless the H production rate is enhanced relative to the SM [230]). Hence, 
a first study proposed to consider the subdominant H ^ A A ^ hhT^T~ final state, with H 
produced via vector boson fusion at the LHC [220,231]. However, here the visibility of the 
signal depends strongly on the poorly known background from 1 1 production. Subsequently 
the Higgs-strahlung process has been added, where the lepton(s) from W , Z decays can 
help to trigger on the relevant events [232-234]. Then, both 46 and 26 2r final states can 
be relevant for both the Tevatron and the LHC [234] (the 2h 2r final state still being more 
promising at the LHC). 

For Ma < 2 m^, the cascade H^AA would end mostly in a 4r final state; clearly, the 
visible r-decay products would not generate narrow peaks in the invariant masses of Ma 
or Mh- At the Tevatron, the prospects to detect the Higgs scalar in this case look pretty 
dim [204] unless the H production cross section is enhanced; the 47 decay mode would 
require enhanced branching ratios for A — )■ 77 [235]. At the LHC particular efforts would 
also be required. Proposals for signals and cuts appropriate for the A A — )■ 4 r — 2 /i-|-2 jets 
final state have been made in [236]; with 100 fb^^ of integrated luminosity, the expected 
rates after cuts are ~ 8 x 10^ from H production via vector boson fusion, and ~ 10'^ from 
H production via Higgs-strahlung {W^ * ^ H + W^) where one can trigger on a lepton 
from decays. In [237] it has been proposed to consider diffractive Higgs production 
{pp —> pp + H) in order to be sensitive to H 4r, which requires to install additional 
forward detectors. Using a track-based analysis in which all events with more than 6 
tracks in the central region are discarded, a viable signal seems possible after accumulating 
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300 fb~^ of integrated luminosity. In [238,239], the subdominant — )■ — t- 2r2/i final 
state was discussed: in spite of the small branching fraction (with 2/i from direct A decays) 
it was argued that, for Mh ~ 102 GeV, the Tevatron can see a signal over background for 
an integrated luminosity £ ~ 10 fb~^, and the LHC already for £ ~ 1 fb~^, with H being 
produced via gluon-gluon fusion. 

Light A production in association with charginos [240] (which requires A > 1 to be 
observable), and from neutralino decays [241] has also been considered. The LHC discovery 
potential for A = 2 - where M^-^^ can be as large as 200 to 300 GeV - has been analysed in 
[129] . Further details of current ATLAS and CMS studies of benchmark scenarios including 
the H ^ A A —i- At final state can be found in [228]. 

The charged Higgs decays — )■ + H or — j- + A play only a marginal role in 
the MSSM. In the NMSSM, larger regions in the parameter space exist (at low tan /3) where 
^ T^t, but where these processes are kinematically allowed and the corresponding 
branching fractions are important [242]. In the case of a light A with a non-negligible 
SU{2) doublet component, the second process (and also pp + A ^ + AA) can 

be observable [243]. 

Higgs searches at the Tevatron have started to test regions of the NMSSM parameter 
space. Recall that, in the NMSSM, a CP-even Higgs H2 with nearly SM-like couplings can 
have a mass up to 160 GeV (if a light singlet-like Hi exists, escaping LEP constraints), 
or even up to ~ 190 GeV, if the singlet-like Hi decays unconventionally (see above). The 
present results of searches for a SM Higgs at the Tevatron (CDF and D0 ) from various 
production and decay channels have recently been combined [244], and exclude a Higgs 
scalar with SM couplings in the mass range 160 < mu < 170 GeV. While such a Higgs scalar 
would be impossible in the MSSM, the corresponding mass range touches the parameter 
space of the NMSSM. 

In [245] the D0 collaboration has looked for CP-even or CP-odd scalars $ produced 
in association with h quarks and decaying into r'^r", assuming a branching ratio of a 
SM-like Higgs scalar. In the mass range 90 GeV< Af$ < 150 GeV, very large values 
of tan/3 > 40 — 80 (which would imply an enhanced production cross section) could be 
excluded. This result is relevant for both the MSSM and the NMSSM, although one usually 
expects tan/3 < 60. 

The D0 collaboration has also searched for H ^ AiAi decays in the range Ma-^ < 
2 rriT- where the branching ratio Ai — )■ fi~^ fi~ is not too small. The 4 /i signature has 
been searched, and an upper limit of about 10 fb on a{pp — )■ H X) x BR{H AA) x 
{BR{A lJi^lJi~)f has been set [246]. Assuming Mh = 120 GeV and BR{H ^ AA) = 
100%, the negative result implies an upper bound on the BR{A — )■ ^~^fj.~) < 10%. The 
final state /i"'"/U.~ t~^t~ (relevant for M^^ > 2 rrir) has also been investigated, but the limits 
are still a factor ~ 4 above the theoretical expectations. Further results from the Tevatron 
can be expected in the near future. 

5.2 The neutralino sector 

The neutralino mass matrix in the general NMSSM is given in (2.32) as function of the 
vevs Vu, Vd and s. In the Za-invariant NMSSM, the last term 2/i' in the (5, 5) element is 
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absent. In [58,247-249] approximate formulae for the mass eigenvalues and mixing matrices 
have been derived. From M^o < M'^ + M| sin^ 9w, and assuming the relations (3.1) for the 
gaugino masses, an upper bound on the lightest neutralino mass as function of the gluino 
mass M3 can be deduced [250]. In the nMSSM (where k = 0, see Section 7.3), the bound 
M^o < \v (with V as in (2.13)) holds [251]. 

The full set of neutralino couplings in the Zs-invariant NMSSM, necessary to study 
neutralino production processes at e"*" e~ colliders, in sfermion (squark/slepton) decays, 
chargino/neutralino cascade decays and Higgs decays, is given in [50]. 

Clearly, in the decoupling limit A^/t— ^-O, s~1/k— s-ooof the Za-invariant NMSSM 
(see Section 2.2.1), it becomes practically impossible to distinguish the neutralino sector of 
the NMSSM from the one of the MSSM: in this limit, the singlino-like neutralino x^s ceases 
to mix with the MSSM-like neutralinos, and all production cross sections as well as partial 
widths of decays into x's tend to zero. An important exception to this rule is the scenario 
where x^s (^^^ '^^ small enough), and R-symmetry is conserved: then all sparticle 

decay chains will first proceed as in the MSSM into the NLSP (which can now be charged!), 
and only at the end of the decay chain the NLSP finally decays into x% This additional 
decay (leading typically to extra leptons in the final state) should allow to distinguish the 
NMSSM from the MSSM in the singlino-LSP regime. Due to a possibly tiny value of the 
NLSP— coupling, the NLSP lifetime can be quite long leading to displaced vertices at 
colliders, see below. Subsequently we briefly discuss neutralino production processes at 
e"*" e~ colliders, in sfermion and cascade decays. 

5.2.1 Neutralino pair production at e"*" e~ colliders 

Neutralino pair production at e"*" e~ colliders can proceed via a Z boson in the s channel 
and via the exchange of a selectron in the t/u channels. Formulae for the corresponding 
cross sections are given in [252] where, however, the "photino/zino" -basis in the gaugino 
sector is employed, differing from the one used here. In [249], the formulae for the cross 
sections are given in the "bino-wino" -basis employed here. 

Present constraints on the neutralino sector originate from its contribution to the invis- 
ible Z width and from bounds on — )> x^ Xj (^ 7^ j) from LEP. Note that bounds on 
the mass of the lightest neutralino Xi relying on MSSM-like relations between neutralino 
and chargino masses cannot be applied to the NMSSM. For earlier discussions of LEP 
constraints on the neutralino sector of the NMSSM see [253,254]. 

Assuming 3 massless neutrinos, additional contributions to the invisible Z width P™^ 
should not exceed ~ 2 MeV (see the Section on z/ mass limits in [89]). This limit has to be 
compared to the partial Z width into neutralinos (if m^o < M^/2), given by 

where the mixing matrix elements Nij are defined in (A. 7). Clearly, the mixing matrix 
elements Ni^ and (proportional to the higgsino components of the neutralinos) are tiny 
for singlino and/or gaugino like neutralinos, which are thus not much constrained. 
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Signals from e+ e — >■ Xj with subsequent x'i decays (i > 1) have been searched for by 
DELPHI [255] and OPAL [256] at LEP, with ^ up to 209 GeV. The non-observation of such 
signals imposes upper bounds on the production cross sections (j(e+e~ — )■ XiX^) ^ 10~^ pb 
and a{e+e- ^ 10"^ pb {i,j>l). 

Phenomenological analyses of pair production of neutralinos in the NMSSM at e+ e~ 
colliders at higher energies have been performed in [249,252,257-260]. Since the information 
on the neutralino sector from the LHC will be quite limited, an ILC-like e+ e~ collider can 
be crucial to distinguish the NMSSM neutralino sector from the one of the MSSM [259], 
although it cannot be guaranteed that the difference is visible if one is close to the decoupling 
limit mentioned above. This question has also been addressed in the radiative production 
of the lightest neutralino pair, e+ e~ XiXil: at an ILC with i/i = 500 GeV in [260]. 

5.2.2 Decays into and of neutralinos 

As in the MSSM, in the NMSSM sfermions fi (squarks and sleptons) can decay into neu- 
tralinos or charginos via fi fi x° [249] or fi fj x"^ , where the latter decay does not 
exhibit any NMSSM specific features. However, due to the modified neutralino sector, 
the processes fi — )■ fi x^ can be quite different in the NMSSM: even a mostly singlino-like 
neutralino can be produced in decays of squarks and sleptons of the third generation, if a 
non-negligible higgsino component induced by a not too small value of A [261] is present. In 
addition, gaugino components of the neutralinos also allow for their production in sfermion 
decays of all generations [249]. 

Higgs bosons could also decay invisibly into neutralinos; these processes can still be 
interesting for Higgs searches via the recoil spectrum of Z decay products in Z* — )■ Z H, 
see above. An exception are GMSB models (see Section 7.2), where the neutralinos can 
subsequently decay into a gravitino and a photon; the present constraints on and the future 
prospects of this Higgs decay mode have been studied in [262]. 

Once produced, neutralinos can decay in many different ways (except for the LSP): 

• into a fermion + sfermion (if kinematically allowed) ; 

• via two-body decays (if kinematically allowed) x^ ~^ x° — > x^ W^, x^ Xj ^) 
where $ is a CP-even or CP-odd Higgs, and via the radiative decay Xi ~^ Xj 7! 

• via three body decays Xi Xj ^"^ ^~ 5 X? ~^ Xj 9 ^ and x'i ~^ X^ ^ ^ (via Z, Higgs and 
sfermion exchange). 

Formulae for the relevant partial widths can be found in [58,249,252]. Note that the decay 
X° — Xj ^1 can be relevant for the search for a light pseudoscalar Higgs [241]. If the 
m^o — m^o mass difference is small (e. g. in the case x^ = X2 ~ bino, Xj = Xi ~ singlino), 
soft leptons from X2 ~^ Xi ^~ can be an important signal for the NMSSM, but their 
detection at the LHC may require low cuts on lepton transverse momenta [263]. 

5.2.3 Displaced vertices 

As stated above, for small A the couplings between the NLSP and a singlino-like LSP in the 
NMSSM can be so small that the large NLSP lifetime leads to macroscopically displaced 



36 



vertices. This phenomenon has been studied first in the case of a bino NLSP in the scale 
invariant NMSSM in [58], and applied to searches at LEP in [264]. Such additional cascades 
in sparticle decays (as compared to the MSSM), with possibly displaced vertices, can also 
occur in more general extensions of the MSSM by singlets [265-267] as, for instance, in 
models designed to explain the CDF multi-muon events [268] where the cascades originate 
from Higgs decays. 

Considering the case of a bino- like NLSP X2 (with a mass ~ Mi) and a singlino-like 
LSP Xi) approximations for the relevant neutralino mixing parameters can be derived for 
small A, and an approximate expression for the (typically dominant) bino to singlino decay 
width via right-handed slepton {Er) exchange can be obtained [58]: 

^ xTn ^ 2 X IQ^'XHh Hv, cj) (5.7) 

where rj = m^o/m^o, u = m^o/m^^ and the phase space integral I{ri^uj) is of 0{1Q~^) for 
masses compatible with universal soft terms at the GUT scale [58]. The important point 
is that the bino decay width is proportional to A^, implying a long lifetime for A — 0. 

In the fully constrained NMSSM with its nearly unique sparticle spectrum (once a dark 
matter relic density complying with WMAP constraints has been imposed, see Section 7.1), 
the LSP is always singlino-like [216,217], and the NLSP is the lighter stau. The stau mass 
must only be a few GeV above the LSP mass in order to allow for a sufficient reduction of 
the relic density via coannihilation, see section 9. Then, again using approximations for the 
relevant neutralino mixing parameters as in [58] , the expression for the stau decay width is 
approximately given as [217] 



r(ri ^ xlr) ^ A^ (aAm - ^m^ , (5.8) 

where the coefficients a and /3 still depend on trif^ and m^o (of the order 0.01 > a ~ 
13 > 0.0001, decreasing with m^J, and Am = mf^ — m^o. The stau decay width can be 
small due to a small value of A and a possible phase space suppression for small Am — m-r . 
For Am ~ a few GeV and A ~ 10""^, stau lengths of flight of 0{mm) are possible [217]. 

In the near future it will be important to study the possible impact of NMSSM speciflc 
displaced vertices for sparticle searches at hadron colliders. 
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6 b physics and the anomalous magnetic moment of 
the muon 



In this Section we briefly discuss NMSSM specific effects in B physics, T physics and 
for the anomalous magnetic moment of the muon. In most cases an extensive literature 
exists on contributions from MSSM-like supersymmetric extensions of the Standard Model, 
which are beyond the scope of the present review. In general, the contributions in the 
NMSSM can differ from those of the MSSM due to the different Higgs and neutralino 
sectors. Correspondingly, most of the results within the MSSM - like constraints on the 
parameters or a possible improved agreement between theory and experiment - are valid 
in the NMSSM as well. 

A notable exception is the region in the NMSSM parameter space corresponding to a 
light CP-odd Higgs scalar Ai. Here, important new contributions "beyond the MSSM" 
can arise, which lead to constraints on this part of parameter space and/or to possible new 
phenomena (in the case of T /rif, physics). 

6.1 B physics 

The following B physics processes have been studied in the NMSSM: BR{B — Xg'-f) 
in [269,270], BR{Bs _and AM^^d) in [269-271], B+ t+u^ in [270], BR{B 

X,in-) in [269,272], and BR{B -fl+l') in [272]. 

In most cases, the NMSSM specific contributions to the B physics observables origi- 
nate essentially from a possibly light CP-odd Higgs scalar Ai, and the fact that a flavour 
violating vertex h-q-Ai (with q = s,d) is generated, amongst others, by squark-chargino 
loops (see, e.g., [269,273]). Here and below, minimal flavour violation is assumed, i.e. 
the only flavour violation originates from the Yukawa sector and is parametrized by the 
CKM matrix. For large tan/3 and \At\, the loop-generated flavour violating vertex b-q-Ai 
is roughly proportional to tan^^S \ At\] as a consequence, the NMSSM speciflc contributions 
induced by a light CP-odd Higgs scalar Ai increase strongly with tan /3 (similar to most of 
the other MSSM-like contributions). Subsequently we briefly discuss the various B physics 
processes. 

BR{B X,7) 

The branching ratio BR{B — )■ ^^7) is one of the most intensively studied quantity in 
B physics. In the past, constraints from b s'-f have been particularly severe, since the 
experimental world average for BR{B — )■ ^^7) was somewhat below the Next-to-leading 
order (NLO) SM prediction, while the beyond the SM (BSM) contribution involving a 
charged Higgs boson in the relevant diagram is positive. 

This situation has changed during the last years: the present world average estimated 
by the Heavy Flavour Averaging Group [274] reads (for a lower cut on the photon energy 
E^> Eo = 1.6 GeV) 

BR{B ^ Xs-f)\ = (3.52 ±0.23 ±0.09) X 10"^ . (6.1) 
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The SM Next-to-Next-to-Leading Order 0{a1) corrections to the total BR{B — )■ Xs'^) 
branching fraction have recently been combined [275,276], giving 

BR{B X,7) Igj^ = (3.15 ± 0.23) x 10"^ . (6.2) 

In [277] the cut E^y > 1.6 GeV on the photon energy has been treated differently, leading 
to an even lower SM prediction: 

BR{B X,7) Igj^ = (2.98 ± 0.26) x 10"^ (6.3) 

This result can be interpreted as a (weak) hint for positive BSM contributions to 6 — 57; 
in any case constraints on the parameter space of supersymmetric models have become less 
stringent. 

In supersymmetric extensions of the SM, positive contributions to BR{B — > Xg'j) arise 
from loops with charged Higgs bosons. Additional contributions involving stop quarks and 
charginos are roughly proportional to tan/3 and to the trilinear coupling At, and can have 
either sign. For an analysis of the corresponding constraints on the parameter space of the 
MSSM see, e.g., [278]; in the meantime, the full two- loop SUSY QCD corrections of the 
MSSM contributions have been computed in [279]. 

NMSSM specific contributions to BR{B — )■ Xg'y) from the extended Higgs and neu- 
tralino sectors arise only at the two-loop level. These were considered in [269] and, in 
more detail, in [270] with the result that the effects are negligibly small, at least for a fixed 
charged Higgs mass (which can be somewhat lower in the NMSSM due to the negative term 
~ in (2.29)). Even contributions from a light CP-odd Higgs Ai are not very important 
here, since Ai appears only in loops and not in a possibly resonant s-channel. 

Also in the case of B^ — r+z/,- the NMSSM specific effects are merely indirect: this 
process is dominated by exchange at tree level in the SM; in the MSSM and NMSSM, 
additional contributions come from charged Higgs exchange, which can be somewhat lighter 
in the NMSSM. 

BR{Bs /i+/i^) 

The most up-to-date experimental result on this process is an upper bound (at 95% 
confidence level) from CDF [280]: 

5i?(5, ^ /i+/i^)|^^p < 5.8 X 10"^ . (6.4) 
This is still one order of magnitude beyond the SM estimate [281] 

BR{Bs /iV") Ism = (3-8 ± 0.1) x 10"^ , (6.5) 

so that there is some room for potentially large new physics contributions. 

The SM contributions originate from box and gauge-penguin diagrams, which are small 
with respect to the experimental bound. The corresponding MSSM contributions can be 
found in [273] and have been generalized to the NMSSM in [269-271]. Here, the flavour 
changing b-s-Ai vertex ~ tan^ /3 can lead to significant contributions arising from penguin 
diagrams, notably if the CP-odd scalar Ai in the s-channel (decaying into /i+/x~) is close 
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to its mass shell, i.e. if its mass is close to Mg_^. (Since the Ai-fi~^-fi~ vertex is proportional 
to tan /3, the Ai contribution to BR{Bs — )■ fi~^fi~) is proportional to tan^/3.) In fact, CP- 
odd Higgs masses equal to Mg^ are always excluded by the experimental bound on this 
process, but the width of the forbidden region depends on tan/3 and \At\. An example for 
a forbidden region in the (M/i^, tan/3)-plane is shown in Fig. 4 below (taken from [270]); 
for a study at very large tan/3 ~ 50, see [271]. 

Information on the mass differences AM^ ^ = mg^ ^ — tub^ ^ originates from measure- 
ments of B meson oscillations. The present result for AMg obtained by the CDF collabo- 
ration [282] is 

AMf P = 17.77 ± 0.12 ps"^ . (6.6) 

The SM prediction depends on the hadronic factor J'b^ and CKM matrix elements; for a 
comparison with BSM contributions, the CKM matrix elements should be determined by 
tree level measurements such that their values are not "polluted" by BSM effects [283]. 

Using \V*V^^\ = (41.3 ± 0.7) x 10-^ and Jb^sJ^s = 0.281 ± 0.021 GeV given by the 
HPQCD collaboration [284], one obtains [270] 

AM™ = 20.5 ±3.1 ps"^ (6.7) 

Hence, despite the large uncertainty, a negative contribution from new physics seems 
favoured. 

The measurement of AM^ is also quite precise [285]: 

AMf P = 0.507 ± 0.004 ps~^ . (6.8) 

However, again the SM prediction suffers from large uncertainties: 

AM|^ = 0.59 ±0.19 ps-\ (6.9) 

obtained with a tree level determination of IV^^^VJ^I = (8.6 ± 1.4) x 10^'^ [283] and a hadronic 
factor fB,\l^, calculated from fs^^f^J {jB,\l^}j = 1.216 ± 0.041 [286]. 

Supersymmetric contributions arise from box diagrams at the one-loop level, but also 
from double-penguin diagrams involving two flavour changing vertices like b-s{d)-Ai [269- 
271], which are particularly important if Ai is exchanged in the s-channel. Now the Ai- 
contribution is proportional to tan^ /3, and leads again to an exclusion of Ai masses near 
Mb ~ 5 GeV, as shown in Fig. 4. In the MSSM, a relation between AMj,(rf) and BR{Bs — )■ 
/x"'"/i~) at large tan/3 can be deduced [287]; as pointed out in [269], this relation is spoiled 
by the Ai-exchange diagrams in the s-channel in the NMSSM. 

BR{B XJ+1-) 

This branching ratio (for / = e or /i) has been measured by Babar [288] and Belle [289] . 
In order to avoid contributions from c-quark resonances, the regions 1 GeV^ < M?^, < 
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6 GeV^ (low) and 14.4 GeV^ < M;+;_ < 25 GeV^ (high) are considered separately: 



( (1.493 ± 0.504+°:^^}) X 10-^ Belle 
BR{B -> X,/+/")io„ = < (1.8 ± 0.7 ± 0.5) x 10"^ BaBar 

I (1.60 ± 0.50) X 10-6 average 



(6.10) 




(6.11) 



The SM analysis has become quite refined (for a review, see [290]): 



jBR{B Xsfi+fi~)f^ = (1.59 ± 0.11) X lO-'^ 

\bR{B X,^+^-)^[^, = 2.40 X 10-^(11°;^^) 

jBR{B Xse+e-)f^ = (1.64 ± 0.11) x 10"*^ 

\BRiB ^ X.e+e")^[J, = 2.09 x 10-^(l^°:i) 



(6.12) 



(6.13) 



These values are well within 2 cr of the experimental measurements. 

In the NMSSM, the process B — t- Xgl'^l' is also sensitive to a light CP-odd scalar [269, 
272]: as in the case of the BR{Bs — j- fi~^fi~), the fiavour changing b-s-Ai vertex proportional 
to tan^ P (with Ai decaying into the pair) can lead to significant contributions, and 
again exclude regions of light Ai masses as shown in Fig. 4. The contributions to the 
total branching ratio {BR{B X,/i+/i-)'=^P = (4.3 ± 1.2) x 10"*^ according to [291]) were 
studied in [272] , but only the dependency on tan /3 was shown explicitly. Other processes 
such as BR{B — )■ XsT^t~) and BR{Bs — )■ l^l^'~f) are interesting for Higgs phenomenology 
as well [272], but here experimental data is not yet available. 

Combined constraints from B physics 

The combined constraints from BR{Bs — )■ fi~^fi~), AM^ and BR{B — )■ X^'-f) on the 
mass of a light CP-odd scalar in the NMSSM have been studied in [270], focusing on the 
dependency on parameters like tan/3, M^^ and At. (In addition, LEP constraints on the 
Higgs sector were taken into account.) In Fig. 4 we show the excluded domains for the 
mass of a light CP-odd Higgs scalar as a function of tan /3 for At = —2500 GeV where, in 
addition, constraints from BR{B — )■ Xs/i^/i") are indicated. 

As mentioned before, constraints from BR{Bs — fi^fi~) and AM^ exclude domains 
for a mass of a light CP-odd Higgs scalar around Mb ~ 5 GeV for all tan/3. Apart from 
the visible strong increase of the excluded region with tan /3, the excluded region decreases 
for smaller values of \At\ due to the smaller fiavour violating Ai-quark couplings induced 
by stop-chargino loops. 

Constraints from BR{B X^fi+fi') from 1 GeV < M^+^- < ^6 GeV and VUA GeV 
< ^- < nib exclude always masses of a light CP-odd Higgs scalar inside these domains 
(where Ai would be on-shell). For larger tan /3, larger regions of M^^ are excluded by this 
observable, although constraints from Bg — )■ are typically more significant. 
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Figure 4: Excluded regions in the tan/3-Myij plane for At = —2500 GeV: the gridded 
region is excluded by LEP, the green (vertically shaded) region is excluded by AM^^s, the 
brown (diagonally shaded) region by BR{Bs — yU+/i~), the region inside the dashed lines 
by BR{B X,/i+/i-) with 1 GeV < M^+^- < VQ GeV or VUA GeV < M^+^- < m,,, 
and the lower yellow corner at tan /3 > 12 and small Ma^ by BR{B — )■ ^^7) (from [270]). 



6.2 T and 775 physics 

Radiative T decays 

Light CP-odd scalars could also be produced in radiative T decays T{nS) — 7A1 
with Ai — I = fi,T [269,292-299]. The branching ratio depends essentially, apart 
from Maj^, on the coupling of Ai to 6-quarks. It is useful to introduce a reduced coupling 
Xfi, which denotes the coupling of Ai to 6-quarks relative to the corresponding coupling 
of the SM Higgs boson. Using the decomposition of the mass eigenstate Ai according to 
(A. 2) in Appendix A (where P/^ are the elements of an orthogonal 2x2 matrix) and the 
/3-dependence of hh as in (A. 11), one obtains 

Xd = tan f3P[^ (6.14) 

where can be of 0{1). (Sometimes is denoted as cos^a-) Hence, for tan/3 ^ 1, 
can also be much larger than 1. 

The BR (T(15) ^ 7A1) is given by the Wilczek formula [300] 

PP(T(lg)^7^i) _ GFmlXW m% 
where a denotes the fine structure constant, and P is a correction factor which includes 



j X P (6.15) 
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three kinds of corrections to the leading-order Wilczek formula (the relevant formulae are 
summarised in [301]): bound state, QCD and relativistic corrections. 

The CLEO collaboration recently published their negative results on the search for a 
CP-odd scalar (with a width below 10 MeV) [302]. Using the Wilczek formula (6.15), these 
results allow to put upper limits on Xd [293,298]. 

Two caveats must be mentioned, however: first, the various (e.g. relativistic) correc- 
tions parametrized by F in (6.15) become large and unreliable for m^^ > 8.8 GeV, where 
F vanishes if the presently known corrections are extrapolated naively. Second, the exper- 
imental search assumed a width of Ai below ~ 10 MeV which could be violated for m^^ 
close to the T mass and/or for very large Xd- 

Nevertheless, CLEO results impose strong constraints on CP-odd scalars with masses 
below ~ 8.8 GeV, which are shown as function of m^^ and Xd (together with other con- 
straints) in Fig. 6 in Section 6.4. 

Mixing between the light CP-odd Higgs and 7]^ states 

If the mass of the CP-odd Higgs is above ~ 9 GeV (as favoured by CLEO constraints), 
but below the BB threshold of ~ 10.5 GeV, a mixing between CP-odd hadronic resonances 
rib{nS) and the Ai can become relevant in case of a large coupling Xd [294,299,303,304]. 
Such a mixing could have a direct impact on the masses of the rjb system. 

The BABAR collaboration has recently determined the rih{lS) mass 11^1^^(13) with an 
error of only a few MeV from radiative decays T — )■ 777^ of excited T states and the 
observation of peaks in the photon energy spectrum. The result from T(3S') decays is 
m^,(i5) = 9388.9l|^ (stat) ± 2.7 (syst) MeV [305], and the resuh from T(25) decays is 
= 9392.9^^8 (stat) ± 1.9 (syst) MeV [306]. The average gives [306] m^^^is) = 
9390.9 ± 3.1 MeV, implying a hyperfine splitting E^fg^lS) = mx(i5) — rrir^^f^is) of 

£;^^P(15) = 69.9 ± 3.1 MeV . (6.16) 

This result can be compared to predictions from QCD. Recent results based on pertur- 
bative QCD are in good agreement with each other and give Egg°(l^) = 44±11 MeV [307] 

and E^'^^ilS) = 39± 14 MeV [308] (see [309] and refs. therein; in phenomenological mod- 
els as quark models Ehfsi^S) varies over a wider range [310]). The discrepancy between 
-E'^j^(lS') and E^^^{1S) could easily be cured in the presence of a mixing of the observed 
rjiy with a CP-odd scalar Ai with a mass somewhat above ~ 9.4 GeV [298]. 

On the one hand, a CP-odd Higgs scalar Ai with a mass very close to 9.389 GeV (before 
mixing) and a strong r/;, — Ai mixing is excluded, since then the mass of the physical 
eigenstate (after mixing, i.e. after the diagonalization of the 2x2 mass matrix) could not 
be given by 9.389 GeV. This implies an upper bound on Xd for Ma^ near 9.389 GeV [298], 
which is also shown in Fig. 6. 

On the other hand, a CP-odd Higgs scalar Ai with a mass below 10.5 GeV can mix 
with all rit,{nS) states (n = 1,2,3), possibly generating both a distorted spectrum as well 
as unusually large branching ratios into r"*" r~ in the rib{nS) — Ai system [299]. Assuming 
that the large observed hyperfine splitting _E^j^(lS') compared to E'^'^J^{1S) is induced by 
a mixing with Ai, Xd can be determined as a function of M^^ (leading to an increasing Xd 
with M^J, and predictions for the masses of the r]h{nS) — Ai system can be made (within 
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errors). Denoting the 4 eigenstates of the rif,{nS) — Ai system by 77^, i = 1 . . . 4, these 
masses are shown in Fig. 5 as function of M^-^ [299]. {M^^ is denoted by in Fig. 5.) For 
clarity we have indicated the masses m^o^^g^ before mixing as horizontal dashed lines. For 
^ 9.8 GeV, the state 772 has a large Ai component, but for Ma^ > 9.8 GeV its mass 
drops below Ma^ due to the strong mixing with 'qti^S). 
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Figure 5: Masses of all eigenstates rii of the rib{nS) — Ai system as function of Mai, once 
m^i is forced to coincide with the BABAR result for mr^^{is) (denoted by mobs! from [299]). 



6.3 The anomalous magnetic moment of the muon 

The measurement of the anomalous magnetic moment of the muon = {g^ — 2)/2 by the 
E821 experiment at BNL has reached a level of precision which is sensitive to supersym- 
metric contributions. The latest experimental value obtained is [311]: 

a;fP = 11659 208.0(5.4)(3.3) x 10"^° . (6.17) 

However, the evaluation of in the SM suffers from a relatively large uncertainty of the 
hadronic contribution (see, e.g., [312-314] and refs. therein). If it is determined using data 
from e^e~ — ?■ hadrons, one obtains a ~ 3 a deviation between the SM and the experimental 
value [312-314] as, for instance, in [315]: 

«r - = (27-7 ± 9.3) X 10"^° , (6.18) 

where the leading uncertainty originates both from the experimental measurement and the 
estimate of the hadronic contribution from e^e~ data. Hence, a positive contribution to 
seems desirable. 
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However, if the hadronic contribution to is estimated from hadronic r decays [312- 
314,316], the discrepancy is reduced. The disadvantage of this approach is that it rehes 
on assumptions on the pion form factor, isospin violating effects and vector meson mix- 
ings. On the other hand, recent (prehminary) measurements using the radiative return of 
e^e~ — 7- 7r+7r~ seem in better agreement with the estimate from hadronic r decays [317,318], 
according to which the discrepancy is reduced to ~ 1.9 cr. 

The various contributions to in MSSM-hke supersymmetric models are reviewed and 
summarised in [319]. In the MSSM, the (dominant) one-loop contributions to originate 
from chargino/sneutrino or neutralino/smuon loops, which increase linearly with tan/3 and 
are positive for a positive /i-term of the MSSM. This phenomenon persists in the NMSSM 
with 11 replaced by /ieg; the effects of the extended neutralino sector in the NMSSM are 
small [320]. 

Studies of in the NMSSM including a possibly light CP-odd Higgs scalar were per- 
formed in [315,321,322]; details of the necessary generalizations of all MSSM-like one- and 
two-loop contributions (see [319] and refs. therein) are given in [315]. The relevance of 
light CP-odd Higgs scalars for was actually first discussed in the context of general Two- 
Higgs-doublet extensions of the SM [323-325]: it is remarkable that positive contributions 
to from two-loop diagrams involving a closed fermion loop can dominate the negative 
contributions from one-loop Higgs-diagrams. The sum of both contributions from a light 
Ai to has a (positive) maximum for Ma^ ~ 6 GeV. Below Ma^ ~ 3 GeV, the negative 
one-loop contribution dominates, which increases the discrepancy with respect to the mea- 
surement. Since the contributions of Ai to are proportional to Xj (its reduced coupling 
to charged leptons), one obtains upper limits on Xd for Mai ^ 3 GeV from (6.18), even 
after allowing for 2 a deviations, which are shown in Fig. 6 below. 

6.4 Combined constraints on a light CP-odd Higgs 

It is interesting to compare the constraints on a light Ai in the NMSSM from CLEO, 
B physics, and from the measured //b(15') mass. Clearly, constraints from B physics and 

(notably at large tan/3) also depend on MSSM-like parameters as Mh±, chargino and 
slepton masses, which we will not review here and for which we refer to the corresponding 
literature. To some extent, the MSSM-like parameters of the NMSSM also affect the con- 
straints from B physics and on a light Ai in the NMSSM. It is possible, however, to scan 
over these parameters and obtain regions in the [X^, v4i)-plane, which are always excluded. 
The results are shown in Fig. 6 taken from [298], where more details can be found. 

We learn from Fig. 6 that, for Ma^ < 9 GeV, CLEO constraints are generally stronger 
than those from Bg — )■ fi'^fJ'~, AMq (g = d, s) and a^, with the exception of a small window 
for Ma^ near Mb ~ 5 GeV. In the range 9 GeV < Ma^ < 10 GeV, constraints due to 
the measured rih{lS) mass by BABAR play a role, but we recall that m^i slightly above 
9.4 GeV can even have a desirable effect in the form of a negative shift of the rij,{lS) mass, 
which would be interpreted as an unusually large T(IS') — ?7b(15') hyperfine splitting. In 
addition, constraints from LEP2 on a light CP-odd scalar below the B — B threshold are 
quite weak (cf. Section 5), which makes this mass range particularly interesting. 

Future searches for a light Ai in this mass range in radiative T decays will be challenging 
due to large backgrounds and the softness of the emitted photon. This problem can be 
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Figure 6: Upper bounds on versus the Ai mass for all parameters scanned over (see [298] 
for details). Indicated are constraints from Bg — )■ fi~^fi~ and AMg, q = d,s as a. green 
dashed line, constraints from blue dashed line, the latest bounds from CLEO on 

BR(T jT~^T~) as a black line and constraints due to the measured rif,{lS) mass by 
BABAR red line. 

overcome if one studies a breakdown of lepton universality in the inclusive leptonic decays 
of the T system as advocated in [298,326,327]: since Ai would decay dominantly into 
T~ , the branching fraction into 2r of an T{nS) state, decaying partially via Ai, would 
be enhanced. 

For a very light Ai {Ma^ < 1 GeV), the decay Ai -> fi"^ fi~ could be dominant. The 
CLEO bounds mentioned above are valid for Ai masses down to ~ 250 MeV, but the recent 
BABAR constraints from the absence of dimuon signals in T{3S), T{2S) — 7^4 can be 
more relevant. These read [328] 

BR{T{3S) ^ 7 A) X BR{A ^+ /i") < 5.2 x 10"^ for 0.212 < Ma < 9.3 GeV (6.19) 

and give [329] 

Xd < 0.4 if BR{A ^ /i+ ^^) ~ 1 . (6.20) 

Very light CP-odd scalars can also show up in Q, S, K, t] and vr decays, which have been 
studied in the NMSSM in [269,295,329-333]. For Ma, < 360 MeV, the decay K+ n+Ai 
is possible; for Mk+ — > Ma, > 2 M^, the result of the HyperCP collaboration [334] 

BR{K^ ttVV^) = 9.8 ± 1.0 ± 0.5 x 10~* (6.21) 

implies X^ (1 — tan~^ /3) < 0.06 [329] within this specific range for Ma,- 
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7 Specific scenarios 



7.1 The constrained NMSSM 

Even if one confines oneself to the Zs-invariant NMSSM with a scale invariant superpoten- 
tial, the number of independent parameters in the Higgs sector at tree level (six according 
to (2.17)) is considerably larger than the two parameters of the Higgs sector of the MSSM 
(typically chosen as tan/3 and M^)- However, one can make the assumption that the soft 
SUSY breaking terms are universal at a large scale like the GUT scale (not very different 
from the Planck scale), as it would be the case in mSUGRA with spontaneous super- 
symmetry breaking in a hidden sector and the mediation of supersymmetry breaking to 
the observable sector via flavour blind gravitational interactions [14]. The corresponding 
version of the Zs-invariant NMSSM is denoted as the constrained NMSSM or cNMSSM. 

Here, as in the constrained MSSM (cMSSM), all soft scalar masses squared, trilinear 
soft couplings and gaugino masses are assumed to be given by common values ml, Aq 
and Ml/2, respectively, at the GUT scale. In fact, the number of independent parameters 
in the cNMSSM is the same as in the cMSSM: the additional parameters fi and B of 
the cMSSM [16] are replaced by the Yukawa couplings A and k, in the superpotential (2.6). 
Apart from the gauge and quark/lepton Yukawa couplings, the Lagrangian of the cNMSSM 
depends on just five parameters ml, Aq, M1/2, A and k, and the correct value of Mz reduces 
the dimension of the parameter space from five to four. 

It is possible that the gravitational couplings (in the Kahler potential in supergravity) of 
the singlet 5* in the NMSSM differ from those of all other fields; then the singlet- dependent 
soft SUSY breaking terms could differ from universality (at the GUT or Planck scale). 
Therefore it is interesting to analyse models with relaxed universality conditions of the soft 
terms, which are slightly more general than the fully constrained NMSSM. 

Below we will discuss first the cNMSSM, where the singlet- dependent soft SUSY break- 
ing terms ms and A^ are assumed to coincide with mg and Aq, respectively, at the GUT 
scale. Subsequently we will turn to the cNMSSM with relaxed universality constraints, 
where the values of ms and/or A^ at the GUT scale are allowed to differ from mo and Aq. 

First analyses of the phenomenologically acceptable regions of the cNMSSM parameter 
space [45,46,53] were handicapped by the still unknown top quark mass, which plays an 
important role due to the impact of the top Yukawa coupling on the RG equations [53] 
and the radiative corrections to the Higgs masses [110]. Somewhat later - but still before 
all present LEP2 constraints were available - the phenomenologically acceptable regions in 
this parameter space were investigated in [54,72,113,197,335,336] (see [337] for an analysis 
of non-universal soft terms from orbifold string theory, and [338] for constraints on soft 
terms in dilaton-dominated SUSY breaking scenarios). 

LEP2 constraints on the parameters of the Higgs sector of the cNMSSM have been 
analysed in [104, 127]. Concerning these, it is important to recall that larger values of 
A do not necessarily imply an increase of the mass of the SM-like CP-even Higgs scalar: 
as indicated in and discussed below (3.2), the mixing of the SM-like Higgs with a heavy 
singlet-like Higgs (also proportional to A) leads to a decrease of the SM-like Higgs mass. In 
the constrained parameter space of the cNMSSM, the off-diagonal matrix elements cannot 
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be fine tuned to 0. As a consequence, LEP constraints on the Higgs sector lead to an upper 
bound on A within the cNMSSM: 

A<0.3. (7.1) 

Then the lower bound on \fj^cs\ = l-^-^l of \fj^cs\ ^ 100 GeV from the non-observation of a 
chargino lighter than ~ 100 GeV at LEP implies 

|s| > 300 GeV . (7.2) 

For |s| ^ Vu, Vdi the dominant s-dependent terms in the potential (2.11) are given by 
(2.19) leading to the condition (2.20) for an absolute minimum with s 7^ [53,54]: 

A^>9m|. (7.3) 

For small A and hence small n from vacuum stability [54], the parameters and 
are hardly renormalised between the GUT and the electroweak scales in the cNMSSM, and 
the above condition translates to (assuming > 0) 

mo < l\Ao\ . (7.4) 

Next, we consider the CP-odd Higgs boson mass matrix (2.27) whose diagonal matrix 
element A^p22 ~ —SkA^s (for large s as relevant here) must be positive. For positive s 
and K this implies negative trilinear couplings 

A. ~ < . (7.5) 

In the region of the parameter space of the cNMSSM where A, k <^ 1, its phenomenology 
would be very close to the one of the cMSSM with corresponding soft SUSY breaking terms, 
since the additional Higgs and neutralino states of the NMSSM decouple and would never be 
produced - unless the additional singlino-like neutralino is the LSP (see below). For larger 
A ~ 0.3 and low tan/3 ~ 2, the mass of the lightest CP-even Higgs scalar with substantial 
couplings to gauge bosons and fermions can be about 10 GeV larger than in the cMSSM. 
However, the lightest CP-even Higgs scalar can also be singlet-like [54,72,113,127,197,336], 
and escape LEP constraints due to its reduced coupling to the Z boson in spite of its small 
mass. The search for such a Higgs scalar would prove to be quite difficult at the LHC. 



7.1.1 Constraints from the absence of charge and colour breaking minima 

Additional constraints on the parameter space follow from the absence of charge and/or 
colour breaking (CCB) minima of the potential, where squark and/or slepton vevs do not 
vanish. A dangerous direction in field space is along the D-fiat direction {Ept-^ \ = \Li\ = \Hd\ 
(with sleptons Eji-^^ and Li of the first generation with the smallest Yukawa coupling he)- 
A deeper minimum is avoided if 

^Al<ml^+ml^+ml^ (7.6) 

at the corresponding scale A^/he [53,339,340]. (In the MSSM, an additional term /i^ has 
to be added to the right-hand side of (7.6) [26]; in the NMSSM, s and hence /^efr can vanish 
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in the dangerous minimum of the potentiaL) In principle, another dangerous direction in 
field space corresponds to |Tr| = IQ^l = \Hu\', however, the corresponding constraints are 
never more relevant than the ones following from (7.6) in the cNMSSM. 

In terms of the universal soft SUSY breaking parameters Aq, M1/2 and mo, (7.6) at the 
scale Ae/he becomes [340] (using the corresponding RG equations) 

{Ao - O.5M1/2)' <9ml + 2.67M2/2 (7.7) 

implying an upper bound on \Aq\. 

Equally delicate could be the so-called unbounded-from-below (UFB) directions in field 
space, which are both Z)-fiat and F-fiat. In [339,340], it has been clarified that such dan- 
gerous directions in the field space of the MSSM are still present in the NMSSM, although 
the singlet vev s gives an additional positive contribution to the potential. Analytic ap- 
proximations to the potential along such dangerous directions have been studied in [341], 
with the conclusion that the inequality 

mo > (0.3- 1.0) Ml/2 (7.8) 

(where 1.0 corresponds to low tan /3, 0.3 to large tan /3) is an approximate condition for the 
absence of deeper minima in these directions. However, since the decay rate of the standard 
vacuum is usually much larger than the age of the universe [341], (7.8) can be violated if 
we assume that the early cosmology (temperature-induced positive masses squared for the 
squarks and sleptons) places us into the local standard minimum of the scalar potential. 

7.1.2 The constrained NMSSM with dark matter constraints 

Next one can require that the LSP of the cNMSSM provides the correct dark matter relic 
density, see Section 9. Early studies of correspondingly allowed regions of the cNMSSM 
parameter space have been performed in [335,342]; in [339], the bounds (7.7) and (7.8) 
from the absence of CCB and UFB minima were taken into account. In the meantime, 
constraints on Higgs and sparticle masses as well on the dark matter relic density [343,344] 
have become tighter (and RGEs/radiative corrections are known to a higher accuracy), 
with the result that the regions in the parameter space of the cNMSSM considered in these 
early studies are no longer phenomenologically viable. 

As discussed above, the UFB constraint (7.8) can in principle be violated, and in the 
cNMSSM this is necessary in view of the up-to-date experimental constraints [216,217]. 
On the other hand, the constraint on mg (7.4) (from s 7^ 0) has to be respected in the 
cNMSSM, and it turns out that only regions where mo <^ M1/2 are phenomenologically 
viable. 

In the cMSSM [16], small values of mo < |^i/2 result in a (charged) stau LSP ri, 
which is excluded. In the cNMSSM, the additional singlino-like neutralino Xi can still be 
lighter than the lightest stau and be the true LSP - this allows for very small (or vanishing) 
values of mo. However, the Xi annihilation rate in the early universe must be large enough in 
order to avoid an excess of singlino-like dark matter. It turns out that, within the restricted 
parameter space of the cNMSSM, no s-channel resonances (Higgs or Z) with a mass twice 
the Xi mass are present, which could enhance the Xi ~ Xi annihilation rate (see Section 9). 
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Figure 7: The Higgs masses in the cNMSSM as a function of M1/2 (in GeV). = 
in the left panel, while mg ~ M1/2/IO in the right panel. From below, the displayed 
lines correspond to the states (blue/dotted), /ig (pink/dashed), a? (full/black) and 03 
(full/red) which is degenerate with the and states (from [217]). 

The only possibility to reduce the Xi abundance is then via coannihilation with ti (now the 
NLSP), whose mass has to be just somewhat above the Xi mass (see Section 9). Then the 
Xi abundance is reduced by the process Xi+^ ^ t^i +^', and the Ti — ti annihilation cross 
section is generally large enough. It should be mentionned that the singlino-like neutralino 
x\ would be practically invisible in direct or indirect dark matter detection experiments in 
this scenario. 

The condition that the mass of ri is close to (just somewhat above) the Xi mass leads 
to rriQ < and Aq ~ — |Mi/2 [216,217]; with these restrictions on the parameters, 

the LEP constraints on the Higgs sector are even more severe and require 

A < 0.02 . (7.9) 

(As discussed in [216,217] and in Section 9, the successful coannihilation of Xi with ri 
implies a lower bound on A of ~ 10~^.) 

With Aq being determined in terms of M1/2 and the strong upper bounds on rriQ and A, 
the spectrum of the cNMSSM is nearly completely determined by Mi/2- tan /3, which is no 
longer a free parameter, comes out relatively large with tan/3 > 25. 

The SM-like Higgs mass is 115-120 GeV for M1/2 > 400 GeV, increasing with M1/2. 
Within this range of M1/2, all sparticle masses satisfy lower bounds from direct searches 
and from precision observables as those from B physics. If one requires that the SUSY 
contribution to the anomalous magnetic moment of the muon explains the ~ 3 a deviation 
between the SM and the experimental value (see Section 6.3), the value of M1/2 should be 
below ~ 1 TeV [216,217,315]. 

In Fig. 7 we display the masses of the neutral CP-even, CP-odd, and charged Higgs 
bosons as a function of the parameter Mi/2- On the left-hand side, we take mo = 0, while 
on the right-hand side we take the maximally viable value (from the dark matter relic 
density) mo ~ ^Mi/2. 

The CP-even Higgs boson with the dominant singlet component is the only Higgs state 
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whose mass depends - apart from M1/2 - on niQ. For small M1/2 it is lighter than the SM- 
like Higgs boson, escaping LEP constraints due to the very small coupling to the Z boson. 
For increasing values of M1/2, its mass increases until it becomes nearly degenerate with the 
SM-like CP-even Higgs state. In this region of parameter space, the singlet-like and SM-like 
Higgs states strongly mix; for a mass of the singlet-like Higgs state of ~ 100 GeV, the excess 
of events at this mass observed at LEP (cf. Section 5.1.1) could be explained [217]. 

The neutralino and slepton mass spectra are shown in Fig. 8. As for the Higgs bosons, 
mo = in the left-hand panel, and mo ~ in the right-hand panel. The two nearly 

degenerate sets of lower lines in both panels correspond to the masses of the Xi singlino-like 
LSP (blue/dotted) and the lighter stau ri NLSP (red/full). The mass difference between 
these two states is smaller than ~ 8 GeV, as required in order to obtain a cosmological 
relic density for the singlino Xi compatible with WMAP. The pattern for the masses of the 
charginos and the heavier neutralinos (blue/dotted lines) follows the one of the MSSM, once 
the proper relabeling of the states is made. Since the effective higgsino mass parameter //gff 
is generally quite large, fies ^ the heavier neutralino states Xi and X5 are higgsino-like 
with masses ~ /ies- The states X2 and Xs are, respectively, bino and wino-like with masses 
m^o 2m^o ^ M2 (with M2 ~ 0.75 M1/2). The charginos Xi and X2 are nearly degenerate 
in mass with, respectively, the wino-like Xs and the higgsino-like X4 5 states. 

In any case, the "smoking gun" for the cNMSSM would be the ri NLSP: all sparticle 
branching ratios into the singlino-like LSP are tiny for A satisfying (7.9), hence all sparticles 
would decay at first into the ri NLSP. Only then the ri will decay into the Xi LSP and 
a r lepton, which will thus appear in every sparticle decay chain. For very small A (still 
larger than 10~^, however) or a very small Ti — Xi niass difference, the ri life time can 
be so large that its decay vertices are visibly displaced [216,217] (see Section 5.2.3). This 
signal would definitively be spectacular and serve to distinguish the cNMSSM from other 
supersymmetric extensions of the Standard Model. 
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Figure 8: Neutralino (blue/dotted lines), selectron (black/full lines) and stau (red/full lines) 
masses in the cNMSSM as a function of M1/2 (in GeV); on the left-hand side mo = 0, while 
on the right-hand side mo ~ M1/2/IO. In both panels the states are ordered in mass as 
m^o < < mg^ < m^o < m^^ < mg^ < m^o < m^o ^ . The charginos Xi" and X2' are 
degenerate in mass with, respectively, Xs and X45 (from [217]). 
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7.1.3 The constrained NMSSM with relaxed universality conditions 

The constrained NMSSM with relaxed universality conditions in the singlet sector, allowing 
for both m| 7^ rriQ and 7^ Aq at the GUT scale, has been first analysed - without 
requiring a dark matter relic density in agreement with WMAP constraints - in [345] . Then 
the bounds (7.4) (from the bound (7.3) on 777,5, A^) and (7.5) no longer apply; notably mo 
is no longer bounded from above. Also the upper limit (7.1) on A is relaxed to A < 0.55. 
The freedom in the choice of allows for a possibly quite light CP-odd Higgs scalar Ai. 
This opens new phenomenologically viable regions in parameter space, where the lightest 
CP-even Higgs mass Mj^^ is well below 114 GeV, but LEP constraints are satisfied since 
Hi decays dominantly into a pair of CP-odd Higgs scalars (see Section 5.1). On the other 
hand, the value of at the GUT scale must be chosen inside a narrow window in order 
to obtain < Ma, < \Mh, [228,345]. 

Since the detection of Higgs-to-Higgs decays will be particularly challenging at the LHC, 
some particular points in the parameter space of the semi-constrained NMSSM have been 
proposed as benchmark points in [228]. These include scenarios corresponding to a light 
CP-odd Higgs scalar Ai\ allowing, in addition, the Higgs masses and to differ 
from ?77o at the GUT scale, scenarios with light CP-even singlet-like Higgs scalars can be 
obtained as well. Both cases require particular search strategies at the LHC, which are 
briefly surveyed in [228] and discussed in Section 5.1. 

A complete analysis of the parameter space of the semi-constrained NMSSM, including 
dark matter constraints combined with those from LEP and other colliders on Higgs and 
sparticle searches, S-physics and the anomalous magnetic moment of the muon, has been 
performed in [346,347]. In contrast to the fully constrained NMSSM, where the small values 
for mo imply a Ti lighter than the bino-like neutralino, the LSP can now either be bino-, 
higgsino- or singlino-like. This, and the possible presence of light Higgs states in s-channels, 
lead to many additional regions in parameter space where the WMAP constraints on the 
dark matter relic density can be satisfied. The corresponding processes, as well as the 
prospects for the direct detection of the neutralino LSP [347], will be discussed in more 
detail in Section 9. 

The potential of the LHC to unambiguously identify the semi-constrained NMSSM 
by observing at least 3 neutral Higgs states is confined to the regime of large values of 
tan/3 ~ 50 and A ~ 0.1, where one can benefit both from the enhanced couplings of the 
heavy Higgs states to 6-quarks and from substantial Higgs doublet /singlet mixings [347]. 

The semi-constrained NMSSM with ms 7^ mo, but imposing = Aq at the GUT scale, 
has been investigated in [348], with the result that the phenomenologically viable regions 
in parameter space are still considerably larger than in the fully constrained NMSSM. 

7.2 The NMSSM and Gauge Mediated Supersymmetry Breaking 

The yU-problem of the MSSM is particularly acute in the framework of GMSB (see [349] for 
a review). The essential ingredients of GMSB models are a hidden or sequestered sector, 
where supersymmetry is spontaneously broken in such a way that the component Fx of 
at least one superfield X - often considered as a spurion without kinetic term - does not 
vanish. In addition, a messenger sector (pi with vector-like Standard Model gauge quantum 
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numbers exists, which couples to X in the superpotentiah 

W = X^i^, + ... (7.10) 

The messengers ipi have a supersymmetric mass Mmoss (the messenger scale), which can be 
parametrized by a vev {X) = M^^ss/'^ of the scalar component of X. The non- vanishing 
vev of Fx induces a SUSY breaking mass splitting rh between the scalar/pseudoscalar 
components of (pi. Typically, dynamical supersymmetry breaking is assumed as the origin 
of Fx 7^ [349], but Kahler potentials in supergravity of the No-Scale type, together with 
Giudice-Masiero-like terms for the messengers, can also lead to models with GMSB [156]. 
Since the messengers (pi carry Standard Model gauge quantum numbers - often in the form 
of a complete SU (5) multiplet - they induce SUSY breaking gaugino masses of the order 

Oil Fx ai /„ . , x 

Mi — (7.11) 

47r Mmess 47r Mmess 

at one loop, and SUSY breaking scalar masses squared of the order 
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at two loops. However, neither a //-term nor a 5/i-term are generated, which would lead 
to serious phenomenological problems as summarised in the Introduction. The simplest 
solution to this problem is the introduction of a singlet superfield S with the scale invariant 
superpotential (2.6). However, since S* is a gauge singlet, the radiative corrections men- 
tioned before will not generally generate a soft mass m| and/or a trilinear coupling of 
the order Msusy ~ [55,350], whereas at least one of these parameters (with m| < 0) is 
required in order to trigger a sufficiently large vev of s, cf. (2.19) - (2.21). 

Numerous proposals have been put forward in order to successfully solve the /i-problem 
in GMSB models within the NMSSM or slight modifications thereof [55,62,71,96,156,349- 
367]: 

(i) Several gauge singlets (or non-renormalisable self-interactions of S) could trigger a 
vev s which is large enough [350-353,357,359,363,364,366,367]. Since gauge symmetries 
do not constrain the possible couplings of singlets, the choice of the terms kept in the 
superpotential should preferably be justified by a discrete (R-) symmetry. 

(ii) If the messenger scale Mmess is large enough, the soft mass m| can become negative 
at low scales - as desired - due to a positive /3-function and the RG evolution; however, 
additional positive terms in the /3-function of m| beyond the ones of the NMSSM (see (B.6)) 
are generally necessary to this end. Such terms appear once S couples to additional light 
matter Qi with vector-like Standard Model gauge quantum numbers in the form SQiQi, 
provided the soft SUSY breaking scalar masses of Qi are large enough [71,350,351,355]. 
Alternatively, the positive terms in the /3-function of m| within the NMSSM (B.6) can 
aheady be large enough, if one allows for more general soft SUSY breaking Higgs (and 
squark) masses. This requires more than one spurion X and S't/(5)-breaking spurion- 
messenger couplings [365,367]. 



53 



(iii) The singlet 5* can couple directly to the messengers tpi. The simplest coupling of 
the form r] Sipiipi generates tadpole terms and [62,96,156,351,353] (defined in (2.1) 
and (2.5)) of the order 



6 




mess 



mess 



(7.13) 



(see (4.1) with A ~ Mmcss)- On the one hand, such terms would trigger a vev s even 
for negligibly small parameters m| and A^; on the other hand, the values of C,f and ^5 
can easily be too large for a messenger scale M^ess ^ ^susy ^ this is nothing but the 
singlet tadpole problem already discussed in Section 4. However, a small Yukawa coupling 
7] ~ 10~^ can render such a scenario phenomenologically acceptable [96]. 

(iv) The singlet S can couple directly to a generalized messenger sector invariant under 
a discrete symmetry such that tadpole diagrams are forbidden [356]. Then parameters 
and m| < of the desired order are radiatively generated to one- and two-loop order, 
respectively, leading to a consistent phenomenology [96,362,363]. 

The phenomenological consequences of the diverse scenarios can be quite different. First, 
the squark and slepton spectrum depends on the quantum numbers of the messengers and 
is generally easy to distinguish from, e.g., mSUGRA. Second, additional light matter fields 
Qi as in (ii) above would typically be observable at the LHC Finally, an approximate 
[/(1)r symmetry in the NMSSM Higgs sector [68-70, 205, 368] in the context of GMSB 
[71, 96] can imply light CP-odd Higgs scalars, into which the SM-like Higgs scalar can 
decay; as discussed in Section 5, this scenario would have an important impact on the 
search for Higgs bosons at the LHC. 

Finally we should mention that the /i-problem is also present in models with super- 
symmetry breaking from extra dimensions. Anomaly Mediated Supersymmetry Breaking 
and models where the MSSM sector couples to a superconformal field theory; again, the 
NMSSM and its variants as described above have been proposed as a way out [369-375]. 



As discussed in Section 4, a solution of the domain wall problem of the scale invariant 
NMSSM invariant under a discrete Z3 symmetry consists in assuming the presence of tad- 
pole terms ~ ^5 and/or ~ C,f in the soft SUSY breaking Higgs potential and/or the super- 
potential. In order to solve the domain wall problem, these parameters can be very small 
such that their impact on the phenomenology can be neglected. On the other hand it has 
been proposed in [166] that notably could be of the order ^5 ~ ^susy ^^^^ [1^7], 
~ Mlugy. Then, the tadpole terms in the potential trigger a vev s ~ MgusY even for 
ft — >■ 0, i.e. the terms ^S^ in the superpotential and ^hiA^S^ in £soft can be omitted. The 
resulting version of the NMSSM, where the singlet only appears in the term XSHu ■ Hd in 
the superpotential and in the tadpole terms ~ and/or ~ ^p, has been denoted as the 
Minimal Non-minimal or new Minimal Supersymmetric Standard Model. 



7.3 The nMSSM 
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The phenomenological consequences of this class of models consist in a possibly quite 
light charged Higgs boson (lighter than the SM-like Higgs boson) into which top quarks 
can decay [167,376], and a neutralino LSP which has always a large singlino compo- 
nent [167,377]. Even so, as shown in [176,378-380], the dark matter relic density can 
well be of the required magnitude: the neutralino LSP can annihilate either through 
the Z resonance (due to its small, but non-vanishing higgsino component) or through a 
CP-odd Higgs resonance in the s-channel (see Section 9). Bounds from Z — x^x? have 
been checked; interestingly, the SM-like Higgs scalar could dominantly decay invisibly into 
XiXi [176]. All additional phenomenological constraints on supersymmetric models can 
be satisfied in this quite "economic" class of models [379,380], and the prospects of neu- 
tralino, chargino and Higgs discoveries at the LHC and the ILC in the nMSSM have been 
worked out in [379]. Finally we recall that the nMSSM also allows for successful electroweak 
baryogenesis [176,177,179] (see Section 4). 



7.4 The fine tuning problem in the NMSSM 

One of the most important motivations of supersymmetric extensions of the SM is the 
solution of the hierarchy problem: if one assumes that the Lagrangian of the SM is valid 
up to a very high scale as Mqut, due to quadratically divergent radiative corrections the 
Higgs mass squared has to be tuned with a relative precision of (M^eak/^GUx)^ ~ 10~^* 
in order to obtain a weak scale (given by the Higgs vev) 14 orders of magnitude below the 
GUT scale. 

In supersymmetric extensions of the SM, the weak scale is naturally given by the SUSY 
breaking scale Msusy (if the /i problem is solved). The details of this relation depend, 
however, on the supersymmetric model considered: on its field content, and on the values 
of the unknown parameters. It is well known that the lower LEP bounds on a Higgs scalar 
with a SM-like coupling to the Z boson (and a large branching ratio into 66), as well as 
the present lower bounds on sparticle masses, induce a "little fine tuning problem" in the 
MSSM: its parameters have to be tuned with a relative precision of ~ 10~^ in order to 
satisfy these constraints. It has been suggested that the NMSSM (or variants thereof) can 
alleviate this problem, which we will discuss below (see [381] for a recent review). 

There are different ways to quantify the amount of fine tuning within a given model: an 
analytic approach consists in identifying the independent parameters of the fundamental 
Lagrangian, in terms of which all particle masses including Mz are determined. To each 
parameter one associates a fine tuning measure defined by 



A 



a 



dluMz 



dlnpa 



(7.14) 



(or Mz replaced by M|), and the "amount of fine tuning" is given by the maximum of all 
Afl. For the SM with a fundamental Lagrangian defined at the GUT scale one would obtain 
Amax ~ 10^"^, but all values much larger than 1 are considered as unnatural. In practice, 
after radiative corrections, the dependence of Mz on all Pa is so complicated that the 
various A^ can only be approximatively estimated which is, however, sufficient to identify 
a potential fine tuning problem. 
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Alternatively one can use codes which compute the particle spectrum numerically in 
terms of input parameters: then one can perform scans over a large number of input 
parameters, and study which proportion of the parameter space is still consistent with 
present experimental constraints. However, the result of such an analysis is related to the 
required amount of fine tuning only if the input parameters correspond to the independent 
parameters Pa of the fundamental Lagrangian. Then, although the result depends somewhat 
on the measure and on the initial range used to scan over the parameter space, its physical 
interpretation is clearer than the one of the formal fine tuning parameters introduced 
above. 

Let us now consider the origins of the "little fine tuning problem", starting with the 
generation of the weak scale in terms of the soft SUSY breaking Higgs masses , f^^Ha 
and Hcs (in the NMSSM): neglecting the radiative corrections (which play no important 
role here), the minimisation equations (2.15) of the Higgs potential imply 



M|^-2/i,\+ ^; :7 'V'""^ . (7.15) 



In the absence of fine tuning, all terms on the right-hand side of (7.15) should be of com- 
parable magnitude, and no large cancellations should occur; hence both /i^g and —rrij^^ 
should not be much larger than 0{M'^). The bound \fies\ ^ 100 GeV from the lower bound 
on chargino masses already starts to generate a slight tension, but the main problem stems 
typically from the value of fn'jj^ once the soft SUSY breaking terms are generated at a large 
scale as Mgut: 

Fortunately, given the large value of the top Yukawa coupling, the dominant positive 
terms proportional to the stop masses squared in the one-loop RGE for rnjj^ (cf. (B.6)) 
generate easily m'jj^ < at the weak scale as desired [19]. However, one typically obtains 
~ —m'^, which is often much larger (in absolute value) than : from the RGBs for 
m|. one finds that rriT is never much smaller than the gluino mass M3 (at Mgusv) which, 
in turn, is bounded from below by ~ 300 GeV from searches at the Tevatron (see [89] and 
refs. therein). 

Clearly these arguments are the same for the MSSM and the NMSSM, and could be 
alleviated in the case of the generation of the soft SUSY breaking terms at a much lower 
scale than Mgut, as possible in GMSB. (See [62,350,354,355,367,382] for discussions of the 
fine tuning problem in GMSB and proposals for solutions within the NMSSM or variants 
thereof.) 

At least within the MSSM, another strong argument for large values of rriT originates 
from the LEP bound on the SM-like Higgs scalar mass, whose theoretical expression was 
given approximately in (3.2) according to which it increases proportional to In ^"t^ In 

the MSSM, large values for rriT are unavoidable in order to satisfy the LEP bound. Albeit 
large stop masses are consistent with the non-observation of stops, they would generate a 
too large value for —Tn'jj^ as discussed above. 

In the NMSSM, two different strategies can be pursued in order to alleviate the fine 
tuning problem without assuming large values of mj-: first, one can try to use the first 
positive term in the second line of (3.2) to push the Higgs mass above the LEP bound. Then 
A should be as large as possible, and tan (3 must be small in order to avoid a suppression 
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by sin^2/3. As discussed in Section 3, within the NMSSM A is bounded from above by 
A < 0.7 — 0.8 if one requires the absence of a Landau singularity below the GUT scale; 
nevertheless, a reduction of the required amount of fine tuning by the LEP bound has 
been observed in [175]. Otherwise one can enlarge the particle content or reduce the range 
of validity of the NMSSM to scales far below Mqut in order to justify larger values of 
A [128-133,136,383,384]. 

Instead of trying to push the Higgs mass above the LEP bound with the help of the 
NMSSM specific terms, it is conceivable that the LEP bound of 114 GeV on the mass of 
a SM-like Higgs scalar H does not apply within the NMSSM (see Section 5): ii H decays 
dominantly into a pair of light pseudoscalars H^AA, the branching fraction for H ^ bb 
can be much smaller than 1, and the Higgs mass can be well below 114 GeV even for 
SM-like couplings of H to the Z boson. If A decays into bb, combined LEP constraints 
still require Mh > 110 GeV [57], hence A should be lighter than 10.5 GeV. Then, Mh 
can be as low as 86 GeV [212] which is easily compatible with (3.2) without large stop 
masses. A reduction of the required fine tuning in the region of the parameter space of the 
NMSSM corresponding to light pseudoscalars has been discussed first in [200,202] and was 
elaborated further in [70,205,215,293]. (However, the latest constraints from ALEPH [213] 
should be taken care of.) 

At first sight, light pseudoscalars appear naturally in the NMSSM in the case of an 
approximate global Peccei-Quinn- or R-symmetry (in the Higgs sector, see Section 2.2). 
It is not always straightforward, however, to realise a consistent scenario in these limits. 
In the Peccei-Quinn limit (k — in the scale invariant NMSSM), two possible minima 
exist generically for the vev of the singlet s: one is the "large s solution" where s ~ 1/k 
(see (2.21)) would be very large. Then, in order to avoid an excessively large fics = 
As, A must be of the order of k which corresponds to the decoupling limit (2.40). Since 
the light pseudoscalar would be nearly a pure singlet, the coupling H AA and hence the 
desired branching ratio would disappear. In the other possible minimum for s where s ~ 
AAa sin 2/3/A^ + m|/f ^ (from the third of eqs. (2.15)), it is difficult to have |/ieff| = | As | large 
enough (above 100 GeV) and to ensure that it is the global minimum of the potential [368]. 

In the R-symmetry limit both trilinear soft terms A\ and are assumed to be small; 
however, since the R-symmetry is broken by the gaugino masses Mj, Ax and receive 
radiative corrections proportional to Mj (cf. the RGEs in (B.3)). Hence it is unnatural to 
assume that Ax and are smaller than these radiative corrections, which requires some 
tuning in the (5 — 10)% range [70]. Hence, at least within the minimal NMSSM without 
enlarged particle content up to the GUT scale, the present experimental constraints also 
require some adjustment of parameters. 

Besides these qualitative arguments for the required amount of fine tuning, probability 
distributions or likelihood analyses in parameter space can be studied numerically using 
scans over input parameters and suitable codes, which compute spectra and couplings 
and allow to check experimental constraints. As stated above, the input parameters should 
preferably correspond to the parameters of the fundamental Lagrangian (at the GUT scale) 
up to an overall mass scale of the soft terms, which can be fixed by Mz- Then, quantities like 
tan P are obtained as output from the minimisation of the Higgs potential. This approach 
had been pursued in the fully constrained NMSSM in [54, 113, 127] with the result that 
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~ 90% of the points in the parameter space still allowed by LEPl are now eliminated by 
fruitless searches for the Higgs boson, charginos and sleptons at LEP2, implying a necessary 
fine tuning of at least ~ 10%. (Also, large values tan (3 > 10 require fine tuned parameters 
at the GUT scale in the cNMSSM.) 

Apart from the investigation of the required fine tuning, such analyses serve to study 
the reach of future experiments in the parameter space within specific scenarios. In the 
semi- constrained NMSSM (where the soft singlet mass squared m| is not required to unify 
with TTT-Q, see Section 7.1.3) such studies have recently been performed in [385,386], where 
also an LSP dark matter relic density compatible with WMAP constraints was imposed. In 
these approaches, the code NMSPEC [345] was used, where tan (3 is an input parameter and 
m| and k are determined by the minimisation equations of the Higgs potential, which can 
hide the required amount of fine tuning, notably for large soft SUSY breaking terms. Large 
values of tan/3 and/or quite heavy sparticle masses appear with relatively high probabilities 
in [385,386], which would imply that the LHC cannot test large parts of the corresponding 
parameter space (whereas the prospects for direct dark matter detection in the future look 
better). It can be expected that, once fine tuning criteria are taken into account, the 
prospects of sparticle detection at the LHC would look much brighter. 
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8 Variants of the NMSSM 



8.1 CP violation 

In the SM, CP is explicitly violated by charged current interactions: complex Yukawa 
couplings lead to one physical phase in the Cabibbo-Kobayashi-Maskawa (CKM) quark 
mixing matrix (5ckm); which can successfully explain the observed direct and indirect CP 
violation (CPV) in the neutral kaon sector [e! je and Ek^ respectively) and in the 5-meson 
sector [89]. On the other hand, the non-observation of electric dipole moments (EDMs) 
of the electron, neutron and atoms (like Hg) puts severe constraints on flavour conserving 
phases, such as the strong CP phase, forcing them to be extremely small. 

Supersymmetric extensions of the SM, like the MSSM and the NMSSM, introduce sev- 
eral potential new sources of explicit CPV: all couplings and fermion masses appearing in 
the superpotential or in the soft SUSY breaking Lagrangian can be complex, leading to non- 
trivial CPV phases. SUSY CPV can manifest itself in several low energy observables: in 
addition to the new mixings in the Higgs sector and the implications for the Higgs spectrum 
and decays, one can expect contributions to the EDMs, lepton polarisation asymmetries 
in semi-leptonic decays, CPV in 5-meson mixings and decays, among many others. The 
measured values of (or bounds on) these observables lead to severe constraints on the new 
phases, which is sometimes denoted as the "SUSY CP problem". 

In particular, flavour conserving SUSY phases can induce very large contributions to the 
EDMs, several orders above the current bounds, even if the SUSY spectrum is relatively 
heavy (~ 1 TeV). In the absence of some cancellation mechanism (as it would be natural in 
minimal flavour violation, such as in the case of scenarios based on mSUGRA or GMSB), 
complying with observation forces the new flavour conserving phases to be very small, 
typically in the range 10"^ - 10"^ This holds both in the MSSM and the NMSSM. 

In principle, generalizations of the Higgs sector of the SM allow for spontaneous CPV at 
the electroweak scale, an attractive scenario at first sight. However, in all supersymmetric 
models (including the NMSSM) where couples to up-type and H^, to down-type quarks 
only, possible phases of H^^d can be rotated away by redefinitions of quark fields, and the 
CKM matrix is real. (This remains true if explicit CPV occurs in the Higgs sector only.) As 
a consequence, SM contributions to the CPV observables are absent, and all the observed 
values must be generated from SUSY contributions, which is difficult (but not excluded a 
priori). Furthermore, the structure of the Higgs potential often forbids stable minima with 
spontaneous CPV as in the MSSM and the Zs-invariant NMSSM (see below). 

Independently from CPV in charged or neutral current interactions, spontaneous or 
explicit CPV in the Higgs sector can have important phenomenological implications: in the 
CP conserving case, scalar and pseudoscalar Higgs bosons do not mix, leading to physical 
states with well defined CP parity. CPV in the Higgs sector may lead to mixings between 
scalar and pseudoscalar states (as the physical states no longer are CP eigenstates): the 
full neutral Higgs mass matrix will no longer be block diagonal, but a 5x5 matrix which 
can be parametrized as 




(8.1) 
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where A^lp 7^ implies CPV. Actually, at tree level no CP violation exclusively within the 
Higgs doublet sector occurs, since it can be "rotated away" by field redefinitions; at tree 
level only NMSSM specific doublet-singlet couplings can violate CP. 

In general, the CPV couplings of the Higgs scalars to SM gauge bosons and fermions, 
to superpartners, and to Higgs bosons themselves can significantly differ from the CP 
conserving case, which requires to re-interpret present bounds on the Higgs spectrum, 
and has consequences for Higgs and sparticle detection at colliders. Finally, given the 
possibility of having an enhanced first order electroweak phase transition in the context of 
the NMSSM (see Section 4.3), additional sources of CPV, beyond those of the SM, could 
play an important role in obtaining successful electroweak baryogenesis, which is difficult 
within the MSSM. 

In what follows, we briefiy review some relevant aspects of NMSSM-like models with 
spontaneous or explicit CPV, focusing on the Higgs sector and the relevant literature. (For 
a review of CPV in connection with non-standard Higgs sectors see [387].) 

8.1.1 Spontaneous CP violation 

Assuming a CP conserving Lagrangian before electroweak symmetry breaking implies that 
all bilinear, trilinear and quartic terms can be chosen real (after possible field redefinitions). 
The Higgs fields can develop complex vevs, so that (2.10) generalizes to 

= e'^- , (i/,) = , (5) = s . (8.2) 

The two physical phases (up to a f/(l)y gauge transformation), 

e = Lfy^^ Lpd + Lf^, 5 = 3(^^ , (8.3) 

open the possibility of spontaneous CPV (SCPV). CP can be spontaneously broken at the 
minimum of the tree level Higgs potential, or only once radiative corrections to the effective 
potential are included - the latter case is sometimes denoted as "radiative SCPV". 

Tree level SCPV in the Z^-invariant NMSSM 

As stated in a no-go theorem in [388], spontaneous CP violation is not possible in the 
Za-symmetric NMSSM with a superpotential as in (2.6) at tree level. Minimising the scalar 
potential 

1/ fi'l "I" fi'l / 2 2\2 , / 2 I \2 2\ 2 I / 2 , \2 2\ 2 I 2 2 

s'^ + v'^v'^ + 2 K XvuVdS cos{9 — 5) 
2 

+2 XAxVuVdS cos 9 + -hiAf^s cos 5 (8.4) 

with respect to the absolute values of the vevs and their phases, one finds that CP violating 
extrema are always local maxima and never local minima of the scalar potential, implying 
negative masses squared for at least one Higgs boson. 
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Radiatively induced SCPV in the Z^-invariant NMSSM 

Radiatively induced SCPV is theoretically possible in the NMSSM [389-391], but these 
scenarios are severely constrained by the non-observation of Higgs pair production off Z 
bosons at LEP: analogous to the MSSM, a very light spin-0, non-Goldstone boson appears 
in the spectrum. The presence of such a light boson is a consequence of the Georgi-Pais 
theorem [392] , and current experimental bounds seem to lead to the exclusion of this class 
of models. 

SCPV m the general NMSSM 

In the general NMSSM, where dimensionful terms are present in the superpotential (2.1), 
the no-go theorem of [388] no longer holds and spontaneous CP violation is indeed viable 
already at tree level [393,394]. Combining the minimisation conditions for spontaneous CP 
violation with the constraints from Ek, it was noticed that the theoretical upper bound 
on the lightest Higgs boson mass becames stronger, leading to Mh^ ^100 GeV (still in 
agreement with LEP bounds due to reduced couplings to SM gauge bosons) [393,394]. 
However, these models re-introduce a //-term in the superpotential so that a solution to the 
/i-problem of the MSSM is not obtained. 

SCPV m the NMSSM with a tadpole term 

As discussed in Section 7.3, the nMSSM contains singlet tadpole terms ~ ^5, in 
the scalar potential and/or superpotential [165-167,376,377], but k = 0. Then, as in the 
MSSM, neither spontaneous nor explicit CPV in the Higgs sector is possible at tree level. 
Keeping k 7^ 0, but adding just a soft SUSY breaking tadpole term ~ to the otherwise 
Zs-invariant NMSSM, allows already to circumvent the no-go theorem in [388], and one 
can find true minima of the scalar potential associated with non-trivial spontaneous CP 
violating phases which are phenomenologically viable [395]. Again some Higgs states can 
be quite light, but not excluded by current bounds due to their reduced coupling to gauge 
bosons. The new Higgs mass matrices, together with the relevant radiative corrections, can 
be found in [395]. 

Problems with SCPV 

However, although elegant and apparently simple to realise, scenarios of SCPV are 
seldom viable: recalling that the CKM matrix is real, there are no SM contributions to 
either flavour conserving (EDMs) or flavour violating CP violating observables. SUSY 
contributions, with 5 and 9 as the only sources of CPV, would be the only means to 
saturate the observed values of e^, e'/e, the CP asymmetry of the Bd meson decay, etc.. 
In particular, the dominant contributions to Ek should arise from the chargino-mediated 
box diagrams [395,396]. In general, complying with observation requires the phases 5 and 
6 to be quite large, as well as maximal left-right squark mixing [396,397]. Large values for 
flavour conserving CP violating phases nearly inevitably lead to sizable contributions to the 
EDMs of the electron, neutron and atoms. Complying with the EDM bounds, either via 
a cancellation mechanism, or a heavy SUSY spectrum (not always possible as excessively 
heavy squarks would preclude saturating the values of Ek etc.) represents the most serious 
challenge to the survival of all models with SCPV in the Higgs sector. 
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Spontaneous CP violation at finite temperature 

As discussed in Section 4.3, additional sources of CP violation in the Higgs sector would 
be desirable for baryogenesis. An interesting hypothesis is to assume that CP is indeed 
conserved at T = 0, and to have spontaneous CP violation at finite temperature in the 
Za-invariant NMSSM. This would allow to evade constraints from EDMs, and to avoid 
a light spectrum in the Higgs sector (arising in the case of radiative CPV, as discussed 
before). Finite temperature effects can indeed trigger spontaneous CP violation inside the 
walls of the propagating bubbles present after the phase transition [180,181]. An extremely 
tiny explicit phase is nevertheless required in the effective potential, in order to lift the 
degeneracy between vacua (each generating identical baryon asymmetries, but of opposite 
signs). Such a phase, typically C(10~^ — 10"^), would give rise to negligible contributions 
to the EDMs. 

8.1.2 Explicit CP violation 

The NMSSM is the simplest SUSY extension of the SM where one can have explicit CP 
violation in the Higgs sector at tree level [45,398]. In the MSSM the explicit CP violating 
phases in the Higgs sector (in fi and B) can be rotated away by a redefinition of the Hu 
and Hii fields; this is not possible in the NMSSM due to the presence of the additional 
singlet couplings. Since now the standard CKM mechanism is assumed to be the dominant 
source of CPV in the quark sector, we will not discuss the SUSY contributions to flavour- 
dependent CP observables. MSSM-like phenomena (e.g. CP violation in the i?-mesons) 
will also not be discussed here, as we will focus on the implications of CP violation for the 
NMSSM Higgs sector. 

Explicit CPV at tree level in the Z-^-invariant NMSSM 

In the NMSSM with a scale invariant superpotential (2.6), the parameters A, k, Ax and 
Ak can be complex. The associated phases only appear in a subset of terms of Vniggs (2.9), 

Vm,7 = H'uH'.S*' - XA, H'^H'.S + ^kA^ S' + h.c. . (8.5) 

By a redefinition of Hu, Ha and S, only one physical phase remains [399,400] which can 
be taken to be the phase of Xk,*: 

(j) = arg(AK*) . (8.6) 

The stationary conditions for the phases induce necessarily complex vevs {6, 6^0) provided 
^ 0, TT [400]. 

Regarding the tree level Higgs mass matrix (8.1), the diagonal blocks and are 
those already given in (2.22, 2.24), once the appropriate redefinitions have been performed: 

XAx^XAx cos 9, kX ^ kX cos{(p + 9 — 6) , — > cos 5 . (8.7) 

The entries of the new off-diagonal block A^|p in (8.1) are all proportional to sin{(f) + 9 — 6) , 
and can be found in [399-403]. However, it is important to stress that, at tree level, explicit 
CP violation only induces scalar-pseudoscalar mixings between Hu,d and S; the scalar and 
pseudoscalar components of the Higgs doublets Hu^d do not mix. 
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Nevertheless the impact of exphcit CP violation on low energy phenomenology is not 
negligible. Although CPV effects are small in limiting cases for the singlet vev (s <^ VuiVd 
OT s ^ Vu,Vd) or very large tan/3, large scalar-pseudoscalar mixings can occur in regions 
where s ~ 0{v) and tan/3 ~ C'(l) [400]. As a consequence, the lightest Higgs mass can 
decrease by 10-30 GeV compared to the CP conserving case. However, avoiding an excessive 
contribution to the neutron EDM would require comparatively heavy gauginos and squarks 
with masses ~ C(TeV). 

Explicit CPV in the 'L^-invariant NMSSM with radiative corrections 

In order to obtain scenarios that are still in agreement with the bounds on EDMs, while 
at the same time inducing a moderate amount of CP violation (as desirable for electroweak 
baryogenesis) , one can consider NMSSM scenarios where explicit CP violation in the Higgs 
sector is induced only through radiative corrections to the Higgs masses and couplings [404] . 
Explicit CP violation in the Higgs sector of the NMSSM at the one-loop level, including 
radiative corrections from third generation squarks, has been studied in [402,403,405,406]: 
due to the possibly complex soft SUSY breaking terms in the squark sector, new CPV 
terms will appear in the neutral and charged Higgs effective Lagrangian via these radiative 
corrections. (The RGEs for the parameters of the NMSSM with explicit CP violation are 
given in [109,407].) 

Considering only the effect of the (dominant) corrections from third generation quark- 
squark loops, the new phases in addition to the three tree level phases in (8.5) are 0^^^, 
</)Aj, arising from At, Ab, respectively [402]. The minimisation conditions (with respect 
to Hij and Sj) reduce the five phases to three. (Note that the complex phases of the 
third generation quark Yukawa couplings can be reabsorbed by redefinitions of the quark 
fields.) An important consequence of taking these higher order effects into account is that 
scalar-pseudoscalar mixings between the two Higgs doublets (absent at tree-level) can now 
occur. 

CP violating effects originating from the chargino sector can also play a role. Contribu- 
tions from loops involving charged particles (W^ bosons, charged Higgs and charginos) were 
taken into account in [405]. It was found that the relative phase between the soft breaking 
SU(2) gaugino mass M2 and A, (pc, yields corrections to every element of the 5x5 neutral 
Higgs mass matrix proportional to siiKpc- Depending on the other low energy NMSSM 
parameters, the chargino induced CP violating corrections can lead to contributions to the 
lightest Higgs boson mass in the range —24 GeV < S^pyMni < 16 GeV [405]. Including 
all radiative corrections, M/^^ max increases with tan/3 (in contrast to the case without CP 
violation, see Section 3.2), and the upper bound of M^^^ max ~ 150 GeV is saturated for 
tan /3 ~ 30 [408]. 

Explicit CPV can also have implications on the production and decay rates of Higgs 
bosons at colliders [403,405]: at a linear e"*" e~ collider {^/s = 500 GeV and 1 TeV), the 
production cross sections for the neutral Higgs bosons in the NMSSM with explicit CPV 
can be much smaller than in the CP conserving case. Explicit CPV in the NMSSM can 
possibly be tested at the ILC with -y/i = 500 GeV, where at least one of the five neutral Higgs 
bosons would be produced via Higgs-strahlung, WW- and ZZ-fusion with corresponding 
cross sections of ~ 12, 15 and 1.5 fb, respectively. The decays of the neutral Higgs bosons 
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can also differ from those of the NMSSM without CP violation in the form of modified 
branching ratios of decays into pairs of ss and cc quarks [408] . 

In the nMSSM (with tadpole terms and k = 0, see Section 7.3), explicit (as well as spon- 
taneous) CPV is impossible at tree level, but radiative corrections from sbottom and stop 
squarks can also trigger CP violation in the Higgs sector. The associated phenomenology 
(mass spectrum and collider prospects) is analysed in [409]. 

In all cases of explicit CP violation in the Higgs sector, EDMs are expected to receive 
potentially large contributions. However, the phase combination responsible for the EDMs 
can be distinct from the one responsible for Higgs mixing, which potentially offers additional 
sources of CP violation for electroweak baryogenesis. It has been suggested [401] that the 
phases can be arranged in such a way that the combination inducing the neutron EDM can 
be suppressed, while the ones affecting the Higgs sector are allowed to be sizable. However, 
this requires a substantial fine tuning between the phases of the Higgs vevs and the phases 
of parameters in the Lagrangian. In order to ensure that SUSY contributions to the EDMs 
are below the experimental bounds, the SUSY spectra (notably gaugino masses) should be 
of C(l TeV) [399,400,410]. 

Modifying the NMSSM by promoting the phases of gaugino masses and the trilinear 
couplings Ai to fields allows to address the strong and the SUSY CP problems simultane- 
ously [411]. Then, however, observable SUSY CPV requires non-minimal flavour structures 
beyond those associated with the Yukawa couplings. 

CP violation in multi-singlet extensions of the NMSSM 

In order to keep sizable CPV in the Higgs sector (contributing to baryogenesis at the 
electroweak scale) without excessively large induced EDMs, multi-singlet extensions of the 
NMSSM have been considered. In [404] the NMSSM is extended by two additional vev- 
less singlets (S" and 5"'), and explicit CPV in the form of complex couplings is allowed 
only between fields which have no tree level couplings to quarks or leptons. Nevertheless, 
through loop corrections involving the new singlets, small CP violating phases can appear 
for the EDMs (consistent with current bounds). 

Spontaneous CP violation in the NMSSM with two singlets S and S", taking into account 
radiative effects, has also been proposed [412]. Both Higgs doublets and singlets can develop 
complex vevs, and one now has three independent phase combinations {6 = y^u+fd+fs+'Ps', 
6 = 3(ps and 6' = 3ips')- Contrary to the NMSSM (with either spontaneous or explicit 
CPV), the addition of the second singlet implies that scalar-pseudoscalar mixings between 
the Higgs doublets can occur already at tree level. Although the Higgs spectrum still 
exhibits light states, it is possible to find regions in parameter space where the lightest 
Higgs mass is as large as 80 GeV, barely viable in view of LEP bounds. A general analysis 
of CP violation at the one-loop level in the Higgs sector of the MSSM extended by an 
arbitrary number singlet fields has been performed in [413]. 
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8.2 R-parity violation and neutrino masses in the NMSSM 



The phenomenological implications of R-parity violation (RpV) in the NMSSM are exten- 
sive (for a review of RpV in the MSSM, see [414]). Here we will briefly summarise some of 
the most relevant consequences, as neutrino mass generation. 

R-parity is a discrete symmetry corresponding to the remnant of a group of continuous 
?7(1)r, transformations, whose charges Rp are assigned as [415] 



Historically, R-parity was introduced in early supersymmetric extensions of the SM in order 
to maintain the baryon [B) and lepton (L) number conservation laws of the SM [13, 18, 
415-417]: conserved R-parity naturally leads to conserved B and L, while RpV requires a 
violation of at least B or L. If R-parity is conserved, R-odd particles are necessarily pair 
produced, the lightest R-odd particle (the LSP) is stable and, if neutral and colourless, a 
good dark matter candidate. 

In the MSSM, the most general renormalisable superpotential including RpV terms is 



where Wapc denotes the R-parity conserving part, and z, j, /c = 1, 2, 3 are generation indices. 
The terms ~ /ij, ~ Xijk and ~ X'^^^ violate L, and the term ~ violates B. (RpV can 
also occur spontaneously through the vev of a neutral R-odd scalar, necessarily a sneutrino 
in the MSSM or NMSSM, which also violates L.) Note that the four superfields Hd and Lj 
have the same gauge quantum numbers, and in the presence of the bilinear terms and 
it is no longer possible to distinguish them unambiguously. 

If B and L are violated simultaneously, this can lead to fast proton decay (and other 
nuclear conversions). Generally, RpV can induce single production of neutralinos and 
charginos, LSP decays into SM particles, flavour violating Z couplings to leptons, vio- 
lation of charged current universality in the lepton and quark sectors, additional sources of 
mixing (and CP violation, if the couplings are complex) in neutral mesons, new rare decays 
as lepton flavour violation (LFV) in the charged lepton sector. The baryon and/or lepton 
number violation induced by RpV can be an ingredient useful for baryogenesis. 

One of the most desirable effects of the L number violation stemming from both bilinear 
and trilinear R-parity violating couplings (if the vevs of the sneutrinos are non-vanishing) 
is the mixing between Higgs bosons and sleptons Lj (in the scalar sector) and between 
leptons and neutralinos/charginos. The latter implies enlarged chargino and neutralino 
mass matrices (respectively 5x5 and 7 x 7, in the MSSM). These new mixtures allow 
to obtain massive neutrinos without the inclusion of right-handed neutrinos, and open 
the possibility of explaining neutrino data (mass squared differences and mixing angles) 
with a minimal superfield content. Bilinear R-parity violating couplings can induce, at the 
tree level, a single massive neutrino via a "seesaw-like" mechanism, where much heavier 
gauginos and higgsinos play the role of right-handed neutrinos. In order to account for 
realistic neutrino masses and mixings, higher order contributions (1-loop), arising from 
trilinear Aj^fc, \[-f^ and/or bilinear R-parity violating couplings, must be taken into account. 




+1 for SM particles including all Higgs bosons 




W^MSSM 



WK^c + ^^^H^■L, + -\,,L,■L^ E%, + X[.,L,-Q^ D^^^ + -\'r,D,-D^Unk , (8.9) 
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In addition to the MSSM RpV terms, the NMSSM superpotential can contain a trihnear 
couphng 

^NMSSM — K S Hu ■ Li , (8.10) 

and a corresponding contribution ~ A\^SHuLi to the soft SUSY breaking Lagrangian. As 
in the MSSM with R-parity violating /ij terms, the couphngs Aj in (8.10) induce numerous 
mixings (see e.g. [418]): the charged Higgs bosons mix with charged sleptons, and charginos 
with charged sleptons, leading to enlarged mass matrices. The neutral CP-even Higgs 
bosons mix with (left-handed) sneutrinos, now leading to 6 x 6 mass matrices. In the 
NMSSM one finds 8 neutral fermions resulting from the mixing of (left-handed) neutrinos 
and neutralinos. 

Note that in the presence of bilinear RpV soft breaking terms, electroweak symmetry 
breaking generally leads to non-vanishing sneutrino vevs, and these have to be taken into 
account in the conditions for vacuum stability (UFB directions and CCB minima) in this 
class of models [419]. The RG equations of the R-parity violating couplings in the NMSSM 
have been given (and analysed for fixed points) in [420,421], and all possible baryon and 
lepton number violating operators are classified in [421]. 

A common phenomenological feature of all SUSY models violating R-parity is the oc- 
currence of an unstable LSP. The lightest neutralino (or the lightest stau) can decay via its 
R-parity violating couplings into SM particles. An unstable LSP can still remain a viable 
dark matter candidate provided its lifetime is sufficiently long (of the order of the age of 
the universe). This possibility has been investigated recently in the NMSSM - including 
constraints from neutrino data - in [422] with the result that a gravitino LSP can be suf- 
ficiently long-lived. If the LSP is short-lived (hence irrelevant to the solution of the dark 
matter problem), it must decay sufficiently fast as not to affect the successful predictions 
of Big-Bang nucleosynthesis. In both cases severe constraints on the R-parity violating 
couplings can be derived [414]. 

At colliders, one can expect displaced vertices from LSP decays [414]. Its branching 
fractions (depending on its composition) can shed light on the R-parity violating couplings: 
if the LSP is a singlino dominated neutralino, it would be possible to distinguish the R- 
parity violating NMSSM from the R-parity violating MSSM at future colliders [423]. 

Let us review some phenomenological implications of RpV, including the possibility of 
adding right-handed neutrino superfields to the NMSSM, and mechanisms of neutrino mass 
generation with and without RpV. 

8.2.1 Massive neutrinos in the NMSSM 

Neutrino masses from R-parity violation 

Dirac neutrino-higgsino masses, which do not require the introduction of right-handed 
neutrino fields, can be generated from the R-parity violating MSSM-like first term ~ /ij 
in (8.9). However, the resulting neutrino mass matrix is only of rank 1 at tree level, i.e. 
only one neutrino is massive, which must be improved by loop effects. Furthermore the 
smallness of the neutrino masses must be understood, which requires some fine tuning of 
parameters. 
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Adding only the terms ~ fii in (8.9) to the NMSSM, the required amount of fine tuning 
can be reduced due to possible compensations between gaugino and singlino exchange 
diagrams (at tree level) [424]. Such a compensation is equally possible if only the NMSSM 
specific terms (8.10) are present [418,424]. Once loop corrections are included, a hierarchical 
spectrum and mixing angles in agreement with observation can be obtained (albeit still at 
the expense of some fine tuning) [418]. Natural values for the heaviest neutrino mass lead 
to bounds on ~ 0(10"^ - lO'^) [418,424]. 

If both terms ~ fii in (8.9) and those of (8.10) are present, two massive neutrino states 
can be obtained at tree- level (in contrast to the MSSM) [424,425]. This offers the appealing 
possibility of reproducing neutrino data without having to resort to loop effects. However, 
again the observed mass hierarchy and mixing angles can be obtained only at the price of 
substantial fine tuning of the parameters: although neutrino data by itself would allow for 
ratios \fii/ fies\ as large as ~ 10~^, potential contributions to LFV processes (such as fi ^ ej 
or /i — )■ eee) constrain the values of |/ij//ieff| ^ 10~^. 

In [424] it has been pointed out that the NMSSM singlet can generate thermal lep- 
togenesis (which would later be converted into a baryon asymmetry of the universe) via 
decays involving the R-parity violating trilinear couplings. The neutrino refraction indices 
(relevant for the neutrino propagation in matter, as in supernovae) have been studied in 
the NMSSM, including radiative corrections, in [426]. 

Right-handed neutrino sector 

If one adds three generations of right-handed neutrino superfields A'^^j (with, by con- 
vention, negative R-parity of its scalar components in contrast to the singlet superfield), 
the superpotential can include the following additional terms: 

and corresponding contributions to the soft SUSY breaking Lagrangian. The first two 
R-parity conserving terms are identical to those present in the MSSM with right-handed 
neutrinos, the third term ~ /tj is R-parity violating, and only the last R-parity conserving 
term is specific for the NMSSM. The fermion and scalar components of Nm will mix in 
general with all neutral fermions and scalars. 

Models with right-handed neutrino superfields Nr and a scale and R-invariant superpo- 
tential (M^'^j = and /tj = in (8.11)) allow for the spontaneous generation of mass terms 
for Nr through the last term ~ A'^^j in (8.11): when the singlet gets a vev s, a Majorana 
mass = \^^^ s (preferably of 0(TeV)) is dynamically generated [427]. These models 

have been shown to possess regions in parameter space where left- and right-handed sneu- 
trinos acquire vevs, so that R-parity is spontaneously broken and effective /xf^-terms in (8.9) 
are generated as well. In this case, the model contains two distinct sources for neutrino 
mass generation (the ordinary seesaw with right-handed neutrinos and contributions from 
the /i^*^-terms) and both solar and atmospheric neutrino data can be accommodated. 

Solving the fi-problem with right-handed neutrino superfields and RpV 

From the point of view of particle content, an economic and interesting possibility 
is to replace the singlet S by the right-handed neutrino superfields Njn (which is not 
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compatible with an R- invariant superpotential). This class of models is also denoted as the 
/ii^SSM [428]. In addition to quark and lepton Yukawa couplings - including the first term 
in (8.11) - the superpotential of the /ii^SSM contains the following terms: 

W = ~X,N'j,,Hu-H, +|%^ (8.12) 

Then the vevs of the scalar components of Njn generate an effective /i-term {fi^s = 
Tlii^i^Ri)-! ^^"^ Higgs vevs as well as the vevs vri generate both Dirac and Majorana 
masses for the neutrinos. The observed neutrino masses originate from two different sources: 
from the Yukawa couplings h^ij in (8.11) roughly of (9(10~^) (and a seesaw mechanism with 
right-handed neutrino Majorana masses of the order of the electroweak scale from the terms 
~ Kj), and an R-parity violating mixing of neutrinos with neutralinos. (In an earlier study 
of the effects of a coupling ~ N'^^HuH^ [429], an explicit Majorana mass term for the right- 
handed neutrinos of the order of the GUT scale was added, and cosmological constraints 
on the lifetime of the neutralino LSP - unstable due to the R-parity violating couplings - 
were derived.) 

The phenomenology of this class of models has been extensively studied in [423,428,430- 
432]; in large regions in parameter space, viable solutions can be found (avoiding Landau 
poles, and in agreement with collider constraints on Higgs and sparticle masses, and current 
neutrino data). In spite of the additional Higgs-sneutrino mixings, the upper bound on the 
lightest Higgs mass turns out to be similar to the one obtained in the framework of the 
general NMSSM [Mh^^ ^ 140 GeV) [430]. Regarding the neutrino masses, both normal 
and inverted hierarchies can be obtained, and two large and one small mixing angle can 
be generated [423] even with a simple flavour diagonal structure for the neutrino Yukawa 
couplings. In [432] it was shown that bimaximal and tri-bimaximal neutrino mixing could 
be accommodated with diagonal and degenerated Yukawas for fi— and r-neutrinos, with 
degenerated vevs for the left-handed sneutrinos. Allowing for complex Higgs and sneutrino 
vevs, spontaneous CP violation can occur in the lepton sector, leading to non-vanishing 
phases in the lepton mixing matrix. 

R-parity conserving NMSSM with right-handed neutrinos 

"Standard" seesaw mechanisms (type I, II and HI) can be also incorporated in the 
R-parity conserving NMSSM, in which case additional higher dimensional operators (of 
dimension 6 or 7) can arise from integrating out heavy states [433] . Even with a TeV-scale 
seesaw mechanism, these can have sizable couplings to SM leptons, and may be detected at 
the LHC. Within minimal left-right extensions of the NMSSM one can aim at a simultaneous 
explanation of the smallness of neutrino masses via the seesaw mechanism, and a solution 
of the tachyon slepton problem arising in Anomaly Mediated SUSY breaking models [434] . 

Within the R-parity conserving NMSSM, the vev of the singlet can be the only source 
for non-vanishing Majorana masses for the right-handed neutrinos [435,436] via the last 
term in (8.11). Remarkably, a right-handed sneutrino is a viable dark matter candidate in 
such a model (see also Section 9.3). 
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8.3 £/(iy-extensions of the NMSSM 

Extensions of the SM gauge group by one (or several) f/(l)' gauge symmetries can arise 
naturally from GUTs (see below), string- inspired constructions [437-440], solutions of the 
yU problem in GMSB [357] etc.; for recent reviews and a more complete list of references 
see [183,441-443]. 

The simplest class of U{iy extended SUSY models involves one additional SM singlet 
S, charged under the additional f/(l)', whose vev is responsible for the U{iy breaking. 
The U{iy charges of the matter and SU{2) doublet Higgs supermultiplets are not unique, 
but typically bilinear Higgs couplings (as the /i-term) are forbidden whereas a XSHuH^ 
term is allowed, in which case f/(l)' models are similar to the NMSSM (sometimes they 
are denoted as the UMSSM). The term ~ kS^ in the superpotential of the NMSSM (which 
stabilises the potential for s — )■ oo) is no longer possible, but its role is now played by the 
U{iy D-terms. Note that the Peccei-Quinn symmetry, that appears in the NMSSM in the 
absence of the cubic singlet term, is embedded in the gauged f/(l)' so that the would-be 
Peccei-Quinn axion is eaten by the Z'. Also the NMSSM domain wall problem is avoided, 
since the discrete Za-symmetry is embedded in the continuous U{iy as well. 

The singlet vev generates simultaneously an effective /x-term (which should not be too 
large in order to avoid fine tuning, see Section 6.4) and the mass of the new Z' boson. Limits 
on Mz' (from direct searches at the Tevatron) are model-dependent, especially on the Z' 
couplings to quarks and leptons (see [89,443] and references therein), and are typically 
Mz' > 600 - 900 GeV for a Z' that decays into SM fermions. 

Loops from matter, charged under both the SM and the U{iy gauge groups, generate 
a, Z — Z' mixing mass term parametrized by the angle azz', which is constrained by elec- 
troweak precision data to be smaller than O{10~^), the exact value depending again on the 
chosen t/(l)' charges. Both constraints from Mz' and azz' require a quite large singlet vev 
s > 0{1 TeV), and hence a sufficiently negative (and large) value of m| (see, e. g., [444]). 

All possible U{iy charge assignments allowing for a XSHuH^ term imply at least mixed 
anomalies between the U{iy and the SM symmetry groups. Their cancellation requires the 
introduction of new exotic fermions (hence superfields) which are vector-like with respect 
to the SM, but chiral under U{iy (see, e. g., [440,445-448]). Then, unification of the gauge 
couplings implies the introduction of additional exotics which are charged, but non-chiral 
under both the SM and U{iy gauge symmetries. These new states do not contribute to 
the anomalies, but do affect the RGEs for the gauge couplings. 

In what follows we briefly review some aspects of NMSSM-like models with one SM 
singlet 5*. Then we refer to some features of more general f/(l)' extended models. 

The most relevant modifications of the Higgs potential with respect to the Za-invariant 
NMSSM are the absence of all terms ~ k, but the presence of a new D term in the potential 
of the form 

Vd^'^' = \9l' {Qs\S? + QH,m^ + QHAH^? + e^oticsY , (8.13) 

where gz' is the U{iy gauge coupling, and Qi are the U{iy charges which obey Qs + Qh^ + 
Qhu = in order to allow for the term XSHuHd in the superpotential. 



69 



After symmetry breaking, the Higgs spectrum consists of a pair of charged Higgs bosons, 
three neutral CP-even scalars and only one pseudoscalar. A detailed overview of the Higgs 
mass matrices for various models, including the dominant top-stop one-loop corrections, 
can be found in [441]. 

Compared to the NMSSM, theoretically and phenomenologically allowed Higgs mass 
ranges have to be reconsidered: due to the additional quartic Higgs self-couplings from 
the D terms in (8.13), Higgs masses are generally larger; in particular, LEP bounds on 
Higgs bosons with SM-like couplings to the Z boson are easier to satisfy. The theoretical 
upper bound Mu^^a.^ on the mass of the lightest neutral Higgs boson increases to Mh^-^s.^ ~ 
170 GeV [449,450]. Radiative corrections from top-stop loops have been included in [449, 
451,452], as well as constraints from gauge coupling unification in [453]. 

Recalling that constraints from Mz' and azz' favour large values of s, the Higgs spec- 
trum is generally similar to the MSSM with a large /x-term [440]: Hi is SM-like, the heavy 
pseudoscalar is approximately degenerate with H2 and H^, completing a full SU{2) dou- 
blet. The mass of the heaviest singlet-like neutral Higgs is typically of the order of Mz'- 
Then, the lower bound on the lightest CP even Higgs mass implied by LEP constraints 
turns out to be in general similar to the one of the MSSM {M^-^ > 90 GeV). 

However, in particular cases one can obtain different configurations for the Higgs spec- 
trum as large doublet-singlet mixing, light singlet dominated states, etc. [443,454,455]. 
In [456] a scenario with two light CP-even Higgs bosons has been proposed in order to 
explain the (light) excesses of events in Higgs searches at LEP (see Section 5). Note that 
in the absence of a light pseudoscalar one cannot relax the lower LEP bound on a doublet- 
like Higgs allowing for its decays into a pair of light pseudoscalars as in the NMSSM, see 
Section 5.1. ^ 

The neutralino sector of the NMSSM is extended by the Z' gaugino Z' with a soft SUSY 
breaking mass M^,, and possible B — Z' mixing terms [457]. A detailed discussion of the 
6x6 neutralino mixing matrix can be found in [378,458-460]. In the limit gz' s ^ ^z'l 
singlino and the Z' gaugino mix to form an approximate Dirac fermion with mass ~ Mz'-, 
which mixes very little with the other four neutralinos. Conversely, in the large M^, regime 
one has a very heavy Majorana fermion Z\ and a much lighter singlino with a "seesaw" 
mass ~ M'^i/M^,. In all cases the new states will likely affect sparticle decay chains. 

CPV has also been considered in f/(l)' extended models [183,408,461,462]. However, 
at the tree level the Higgs sector is CP conserving, and spontaneous CPV (after including 
radiative corrections) cannot be realised in the allowed parameter space [408]. Explicit 
CPV in the Higgs sector can be possible via radiative corrections, through the explicit 
CPV phases present in the stop mass matrix or the additional phases in the exotic quark 
sector of the model. These can significantly alter the mixing between the pseudoscalar and 
the heaviest scalar [462]. Due to the extended neutralino sector, phases in A and in the 
slepton trilinear couplings can contribute to the electron EDM [461]. EDM bounds strongly 
constrain the effective phase in the chargino mass matrix, but all other SUSY phases remain 
largely unconstrained, thus alleviating the SUSY CP problem [462]. 

f/(l)' extended models have also important cosmological implications, since the extra 
states can modify the nature of the LSP. Early studies [463] suggested that Z' mediated 
neutralino annihilation could provide important (or even dominant) contributions to the 
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correct relic density. More recent analyses [183, 378, 464-467] (for a general discussion, 
see [468] and references therein) have shown that compatibility with the WMAP bounds 
on the relic density can be achieved in the case of a singlino-like LSP, with a small higgsino 
component, via annihilation into a Z'\ enhanced couplings to the Z also allow annihilation 
via a Z resonance [464]. The prospects for observation of the LSP in dark matter direct 
detection experiments have been considered in [467]. These models offer also non-neutralino 
dark matter candidates: scalar right-handed sneutrinos can be viable dark matter candi- 
dates due to the possibility of annihilation through the Z' [469]; other possibilities include 
exotic LSPs [470], see also Section 9.3. 

Dirac neutrino masses in the experimentally preferred range can also be generated 
within [/(l)'-extended generalized SUSY models [471-474]. If the ?7(1)' symmetry for- 
bids the Yukawa couplings to right-handed neutrinos in the superpotential, effective Dirac 
masses can be generated from (hard) SUSY breaking Yukawa couplings or, at the loop level, 
from non-holomorphic trilinear couplings to the "wrong" Higgs boson, which are naturally 
small [473]. SUSY models of sterile neutrinos can be naturally realised in ?7(1)' extended 
models [472]: sterile neutrino masses of (!?(leV) and mixings among the active and sterile 
neutrinos can then be obtained via the high-dimensional operators, generated by integrat- 
ing out the heavy fields. In extended models where the vev s generates Dirac neutrino 
masses through a dimension-4 term in the superpotential, the decays of the LSP and NLSP 
(the two lightest sneutrinos) could explain the PAMELA anomaly [475]. 

We recall that, once R-parity is violated, proton decay can be excessively fast since 
both baryon number and lepton number are violated in general. In models with 

RpV, the proton can still be sufficiently stable in a natural way, since the f/(l)' symmetry 
can forbid the renormalisable (dimension 3 and 4) and the most dangerous dimension 5 
operators [445]. In [448] models are discussed which can still contain either lepton number 
violating or baryon number violating renormalisable interactions, but never simultaneously 
both. In [476] models have been put forward where bilinear lepton number violating terms 
in the superpotential allow to accommodate the observed pattern of neutrino masses and 
mixings, but the most dangerous dimension 3 and 4 baryon number violating operators 
(and, possibly, all = 1 operators) are absent. 

Production (and decays) of the additional Z' gauge boson and exotic states at colliders 
are reviewed in [443,477]. The exotic quarks and squarks could appear in decays of the Z', 
or may manifest themselves indirectly, since the production of the heaviest Higgs boson (via 
gluon fusion) can be dominantly mediated by loops of these exotics [452]. The production 
of neutral Higgs bosons in this class of models has been discussed in [478]. 

At e~ colliders of higher energy, the Higgs sectors of f/(l)' extended models and 
the NMSSM can be distinguished in the limit of large trilinear Higgs couplings [199]. At 
the ILC, measurements of the Z Z couplings may help to distinguish scenarios with a 
light, leptophobic Z' from the NMSSM, since the couplings are expected to be substantially 
smaller in the former case [479] . In the case of TeV Z' bosons, searches for Z' decays into 
higgsinos can also test if a gauge symmetry is indeed responsible for the absence of a bare 
/i-term in the superpotential [480]. 

There are innumerous possibilities for the construction of more general ?7(1) '-extensions 
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of the MSSM including, for instance, flavour symmetries [481]. A systematic analysis and 
classification can be found in [445]. In what follows we consider some particular cases, some 
of which are motivated by GUTs or string-inspired constructions. 

U{iy and a secluded singlet sector 

A possible way to avoid the tension arising from a small enough /igg and a sufficiently 
large Mz' is to consider the secluded f/(l)'-extended MSSM (sMSSM), where a separation 
between /ies and Mz' is achieved by the introduction of four SM singlets S, Si, S2 and 
5*3 [454]. The ordinary sector contains the NMSSM (or USSM) singlet S, while the three 
additional singlets Si belong to a secluded sector which couples to the ordinary sector only 
via U{iy gauge and possibly soft SUSY breaking terms. All singlets are charged under 
U{iy such that their vevs contribute to Mz', but only S generates a fics of the order of the 
electroweak scale. The SM singlet-dependent part of the superpotential is given by 

W^^'^ = XSHdHu+\sSiS2Ss. (8.14) 

The soft SUSY breaking Lagrangian contains soft SUSY breaking masses and trilinear 
couplings for the new singlets {m'^., rn^SiSji ^^s)- Unwanted global symmetries have to be 
broken explicitly to avoid the potential occurrence of two massless experimentally excluded 
Goldstone bosons in the spectrum. For appropriate parameters, 5'i,2,3 acquire large vevs 
along nearly F-flat and D-flat directions of the scalar potential so that one obtains Mz' S> 
/ieff in a natural way. A can be larger than in the NMSSM (without leading to a Landau 
singularity), since the additional gauge coupling in its /3- function slows down its increase 
towards high scales. As a consequence, the lightest CP-even Higgs mass could be as large 
as ~ 174 GeV already at tree level [454]. In general the Higgs and neutralino sectors 
of these models are quite complicated, and we refer to [455] for a detailed discussion. A 
singlino-higgsino LSP can account for an acceptable cold dark matter relic density [464]. 
Regions in parameter space exist, that can explain the experimental deviation of the muon 
anomalous magnetic moment from the SM and yield an acceptable cold dark matter relic 
density without conflict with collider experimental constraints [482]. 

In the sMSSM, explicit CP violation in the Higgs sector can be induced at the tree 
level by the non-zero phase of fn'g^g^ [483]. Large values of the latter phase can allow for 
a neutral Higgs boson lighter than ~ 90 GeV, compatible with LEP constraints due to 
suppressed couplings to the Z boson. Unlike the MSSM and NMSSM, SCPV can occur at 
the tree level. The possibility of electroweak baryogenesis in the sMSSM has been discussed 
in [182]. 

Additional f/(l)' in string- inspired Eq 

As mentioned before, gauge extensions of the SM gauge group by one (or several) non- 
anomalous U{iy gauge symmetries can naturally arise from a (string-inspired) GUT gauge 
group as Eq (see, e.g., [484-486]). 

In an Eq GUT, the matter sector includes three 27-plets containing the three ordinary 
families of quarks and leptons (including right-handed neutrinos), three families of candi- 
date Higgs doublets [Hdi, Hui) and singlets (Si), and three families of extra colour triplets 
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{D[^ D'-). Anomalies are cancelled generation by generation within each complete 27 repre- 
sentation. Gauge coupling unification requires the addition of a further pair of Higgs-like 
multiplets, H' and H', arising from incomplete 27' + 2/ representations. The breaking of 
Eq occurs via 

Ee 50(10) X f/(l)^ ^ SU{5) x U{1)^ x f/(l)^ . (8.15) 

The low energy gauge group is assumed to be SU{3) x SU{2) x f/(l)y x f/(l)', where 
U{iy = f/(l)x or a specific linear combination of U{1)^ and U{1)^ motivated by the 
breaking of Eq by Wilson lines. 

In [484,485], patterns of f/(l)' and electroweak symmetry breaking were studied in this 
class of Eq models, i^g-type relations between the Yukawa couplings are not imposed, 
and only those conserving baryon and lepton number are kept. Then, the most general 
superpotential allowed by gauge invariance is of the form [485] 

= hl.Ut Q, Huk + /^J, A' Q, H,, + A,,fc H,, H^, + /ig,,5, D'. D'^ , (8.16) 

where /c are family indices. When some of the Higgs fields acquire vevs, electroweak 
and U{iy symmetries are broken. It is always possible to work in a basis where only one 
family of Higgs doublets and singlets develops a vev: one can define 8 = 8^, = and 
Hu = Hu3, while the remaining families are treated as additional exotic doublets and singlets 
without vevs. In order to obtain an acceptable low energy phenomenology, universality at 
the GUT scale of the soft SUSY breaking masses has to be relaxed. (Universal boundary 
conditions for the soft breaking terms could be accommodated by identifying f/(l)' with a 
more general linear combination of U{1)^ and U{1)^ [485].) 

Another E^ inspired version exhibiting an NMSSM-like low energy spectrum and al- 
lowing for explicit CP violation in the Higgs sector has been considered in [487]. The 
electroweak gauge symmetry of the SM is enlarged to 8U{2) x f/(l) x f/(l)i x f/(l)2, where 
both extra U{1) symmetries originate from Eq. The Higgs sector contains two doublets and 
two singlets (only one being allowed to couple to the doublets due to the underlying Eq 
symmetry). Due to top-stop mediated radiative corrections the model allows for explicit CP 
violation in the Higgs sector at the one-loop level, leading to scalar-pseudoscalar mixings 
among the five neutral Higgs bosons. 

The Exceptional SSM 

Specific realisations of E'e-inspired U{iy extensions offer the possibility of incorporating 
a seesaw mechanism (to explain the smallness of light neutrino masses) in a natural way. To 
this end the extra f/(l)' surviving at low energy (denoted by U{1)n) must correspond to a 
particular combination of U{1)^ x U{1)^, in (8.15) under which the right-handed neutrinos 
(also contained in the 27 representation) transform trivially: 

U{1)n = U{iy = U{1)^cos9n + U{l)^smeN with = arctan . (8.17) 

This is the so-called Exceptional SSM (EgSSM, see [488-495] and references therein). 

As discussed in [488] , due to the vanishing charges under U{l)i\f, right-handed neutrinos 
can acquire very heavy Majorana masses suitable for the "standard" seesaw mechanism. 
This also opens the possibility of baryogenesis from leptogenesis, and avoids constraints on 
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the mass of the Z' arising from nucleosynthesis. Although B — L is automatically conserved 
in i?6-inspired SUSY models, some of the new Yukawa interactions could violate baryon 
number resulting in rapid proton decay. Proton stability, successful leptogenesis, non-zero 
neutrino masses and the absence of exotic coloured dark matter relics can be guaranteed 
with the help of a new conserved R-parity-like quantum number and corresponding baryon 
and/or lepton number assignments for the exotic particles. Furthermore, it is mandatory 
to impose a hierarchical structure for the additional Yukawa couplings as well as a discrete 
Z2 symmetry as to avoid excessive flavour changing neutral currents. 

Due to the modified (two-loop) beta functions in the EgSSM, unification of the gauge 
couplings can be achieved for values of as{Mz) closer to the measured central value than 
in the MSSM [490]. 

The EgSSM Higgs sector includes three CP-even, one CP-odd and two charged states. 
In the MSSM limit A <C 1, the new states become very heavy and decouple, rendering the 
model essentially indistinguishable from the MSSM. The analysis of RG flow of the gauge 
and Yukawa couplings shows that the absence of a Landau singularity allows for larger A 
and lower tan /3 at the electroweak scale than in the MSSM and NMSSM, which affects 
the Higgs, neutralino and chargino spectra. When A > (71, the lightest Higgs scalar can 
be somewhat heavier than in the NMSSM; including two-loop corrections, the maximally 
allowed value is around 150 GeV. Vacuum stability constraints push the masses of the 
heaviest CP-even, CP-odd and charged states above 1 TeV, so that only the lightest scalar 
would be within reach of the LHC and an ILC [489]. 

The superpartners of the Z' and the singlet contribute to the neutralino spectrum, while 
the chargino sector remains unchanged. For X ^ gi, the heaviest neutralino and chargino 
states are the neutral and charged Hu,d higgsinos; the lightest chargino is wino-like, while 
the LSP is a bino-like neutralino [488]. The dark matter prospects for this class models 
have been addressed in [491] (see Section 9.3). 

The many new additional exotics form vector-like multiplets which could possibly be 
produced and detected at the LHC; the new Z' could be visible if lighter than 4-5 TeV. 
Within the constrained EgSSM (assuming universal soft SUSY breaking parameters at 
high energy) a detailed study of the phenomenology and impact for LHC discovery has 
been carried out in [492,493]. Scenarios for early LHC discovery were investigated in the 
low (mo, M1/2) regime, and a set of benchmark points was proposed. 

Eg models based on a A27 family symmetry broken close to the GUT scale allow in 
addition to address the flavour problem: quark and lepton masses and mixing angles can 
be accounted for, with tri-bimaximal neutrino mixing [494,495]. 



74 



9 Dark Matter in the NMSSM 



As in the MSSM, the hghtest neutrahno in the NMSSM is a candidate for cold dark matter 
(DM) in the form of WIMPs (weakly interacting massive particles), whose relic density 
should be in agreement with the constraints from WMAP, 0.094 < Qh"^ < 0.136 at the 2 a 
level [343,344]. (For a review of DM in the MSSM see [496].) Due to the differences in the 
neutralino and Higgs sectors, the properties of the LSP in the NMSSM can be significantly 
distinct from those in the MSSM, both in its nature and/or in the processes relevant for the 
LSP relic density and its detection: the LSP can be dominantly singlino-like, and/or the 
additional Higgs bosons can contribute to the LSP annihilation and detection processes. 
A neutralino LSP in the form of a higgsino-singlino mixed state in the NMSSM, giving 
rise to the required dark matter density, was considered first in [335,497-500], where often 
universality constraints on the soft SUSY breaking terms were imposed. 

The large number of processes relevant for the computation of the LSP relic density 
(annihilation and coannihilation cross sections) and of LSP detection in the NMSSM will 
not be discussed here in detail, but all of them are included in the code MicrOMEGAS (see 
Appendix D), and described in [501-503]. 

Subsequently we discuss some aspects of the LSP relic density and its detection in the 
NMSSM which differ from the MSSM: in Section 9.1, we discuss the processes relevant for 
the dark matter relic density; in Section 9.2 we consider the detection of an NMSSM dark 
matter candidate; finally, we will address specific NMSSM scenarios for DM in Section 9.3. 

9.1 Neutralino relic density 

The basic mechanisms that lead to a reduction of the neutralino relic density are essentially 
the same as in the MSSM: neutralino annihilation via s- and/or t-channel exchanges into 
gauge boson, Higgs boson and fermion pairs, and coannihilation of the neutralino with a 
heavier state (typically the NLSP). The cross section for these mechanisms can be strongly 
enhanced by s-channel resonances, occurring when ttt-lsp + '^(n)lsp is close to the mass of 
the particle exchanged in the s-channel. 

Distinctive NMSSM scenarios can be manifest in several ways, depending on the regions 
in the NMSSM parameter space: the singlino component of the LSP can be sizable (e.g. 
mixed higgsino-singlino LSP), and one can even have (nearly) pure singlino LSPs. Secondly, 
the extra scalar and pseudoscalar Higgs bosons can lead to rapid LSP annihilation through 
s-channel Higgs resonances. Moreover, in the case where the additional scalar /pseudoscalars 
are light, new final states can be kinematically allowed: for example, annihilation into Z Hi, 
Hi Hi, Hi Ai and Ai Ai can significantly contribute to the reduction of the neutralino relic 
density, either via s-channel Z, Hi, Ai, or t-channel (heavier) neutralino exchange. Here 
the NMSSM specific couplings A and k can play an important role. 

A detailed discussion the viability of DM candidates in the (unconstrained) NMSSM, 
using recent cosmological constraints and numerical tools, is given in [501]. In what follows, 
we briefly discuss some illustrative examples of how the new features of the NMSSM are 
manifest, separately considering annihilation and coannihilation. 
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9.1.1 Annihilation 

Here we briefly consider neutralino annihilation processes that reflect the new features of 
the NMSSM. For a detailed discussion, albeit in the framework of the MSSM, we refer 
to [496]. Detailed formulae for (MSSM) annihilation cross sections can be found, for exam- 
ple, in [504]. 

• Annihilation via s-channel Higgs exchange (CP-even or CP-odd) 

The annihilation cross section is proportional to (neglecting interference terms) 

2 

/,(s,mJo,m|,m^,) , (9.1) 

where A, A' denote SM fermions, gauge or Higgs bosons, HaXiXi the Higgs-neutra- 
lino-neutralino couplings and HaXX' the Higgs couplings to the flnal states (see 
Appendix A). is the width of Ha, and we have introduced a generic function fs 
that depends on the s-channel momentum and on the masses of the initial and flnal 
states [504]. A similar expression is obtained for the exchange of a pseudoscalar (with 
the appropriate replacements in the couplings, masses, and different fs). 

• Annihilation via s-channel Z exchange 
In this case, one approximately has 

2 

5(^(s,mJo,m^,m^,,M|) , (9.2) 

where the SM ZXX' couplings are proportional to electroweak gauge couplings, and 
{ZxiXi) = 5-2/(2 cos 61^^) {N^^-Nf^}, with A13 and A14 (according to (A.7)) denoting 
the higgsino components of the lightest neutralino Xv Fo^' the function Qs see [504]. 

• Annihilation via t-channel neutralino exchange 

As in the MSSM, neutralino annihilation can also take place via t-channel sfermion, 
chargino and neutralino exchange. In the latter case, annihilation into a pair of singlet- 
like Higgs bosons {HiHi, AiAi, or HiAi) can play a signiflcant role for singlino- 
like LSPs [501]. As an example, the cross section for XiXi ~^ HiAi via t-channel 
neutralino exchange reads 

a{x'lx'AH^A,) oc (i/a?X°)' (^iX?X°)'/^s(s, mjo, m^^, m^J ; (9.3) 

the scalar /pseudoscalar Higgs-neutralino-neutralino couplings are given in (A. 14) and 
hs denotes an auxiliary function encoding the mass dependence [504]. In this case, hs 
decreases with the neutralino mass as m~o. 

Xi 

• Resonant s-channel annihilation 

LSP freeze-out occurs at temperatures T ~ mLSp/20, where the LSPs are non- 
relativistic and s ~ 4mLSP- neutralino annihilation proceeds via the s-channel 
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exchange of a particle (Higgs or gauge boson) with mass ~ 2m^o, the pole in the 
annihilation cross section can lead to a rapid decrease in the neutralino relic den- 
sity. This resonance phenomenon can allow for compatibility with WMAP bounds 
in cases where the relic density would be otherwise too large. It is important to 
stress that for the case of resonant annihilation, the relic density is very sensitive 
to the mass differences of the LSPs and the exchanged particles. Hence, regions in 
the NMSSM (or MSSM) parameter space where WMAP compatibility is achieved 
via resonant annihilation typically correspond to narrow bands (for the NMSSM see, 
e.g., [346,347,501]). 

• Coannihilation LSP-NLSP 

In the MSSM and in the NMSSM scenarios exist where no efficient annihilation mech- 
anism is available. However, the correct relic density can still be obtained by coan- 
nihilation processes. If the NLSP (typically the second lightest neutralino X2 ^ 
sfermion) is not much heavier than the LSP, so that it is still present in the thermal 
plasma at the time of LSP freeze-out, coannihilation processes can efficiently reduce 
the LSP relic density. 

In the NMSSM, coannihilations with heavier neutralinos or charginos can occur via 
s-channel Z, or Higgs exchange, for example Xi + X2 ""^ ~^ XX', or via 
t-channel x'i^xf exchange [501]. For a sfermion NLSP (usually a light, mostly right- 
handed stau), s-channel and t-channel (the latter via neutralino or chargino exchange) 
coannihilation lead to fermion final states. Just like in the case of s-channel LSP-LSP 
annihilation, resonances can also increase the LSP-NLSP annihilation cross section: 
in this case the mass of the exchanged particle should be ~ ^lsp + rn^i^sp- 

• NLSP-NLSP annihilation 

This process, usually also labeled coannihilation, occurs when the dominant anni- 
hilation process of R-odd particles is via NLSP annihilation: NLSP + NLSP — )■ X. 
Provided the NLSP-LSP mass difference is small, the correct LSP relic density is 
achieved by processes that maintain the LSP and the NLSP in thermal equilibrium 
(implying tzlsp ^ ''^nlsp for the corresponding abundances): LSP + X ^ NLSP + X', 
where X,X' are light SM particles (quarks and leptons). 

The reaction rate for the NLSP annihilation process NLSP + NLSP — )■ X is given by 
(nNLSp)^cr, where a is the thermally averaged cross section. On the other hand, the 
reaction rate for the processes that help to keep the LSP and the NLSP in thermal 
equilibrium (LSP + X — NLSP + X' and its inverse) depends on jt-lsp nx o' . Even if 
the latter cross section a' is possibly suppressed (as, for instance, for a singlino-like 
LSP in the NMSSM), the process LSP -FX NLSP -hX' is typically faster than the 
annihilation process NLSP -|- NLSP — ?■ X, since near the freeze-out temperature the 
abundance nx of quarks and leptons is ~ 10^ larger than the abundances of the LSP 
and NLSP (for mLSP ~ ''t^nlsp) [505]. This allows to dilute the LSP density as fast as 
the NLSP density, and such "assisted coannihilation" can play a very important role in 
reconciling the relic density of an NMSSM LSP with WMAP measurements [346,501]. 
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In the MSSM, the LSP is usually a bino- or a higgsino-like neutralino (or a bino-higgsino 
mixed state). These scenarios can also be found in regions of the NMSSM parameter space 
where, apart from the enlarged Higgs sector and heavy neutralinos, the spectrum is MSSM- 
like (see [346,501]). Before turning to NMSSM specific annihilation processes, we briefly 
discuss these MSSM-like scenarios. 

The relic density of a nearly pure bino-like LSP is typically difficult to reconcile with 
WMAP measurements: the only available annihilation channel for a pure bino is via t- 
channel sfermion exchange (as all couplings to gauge bosons require a higgsino component), 
or via coannihilation with a sfermion (typically stau or stop) NLSP. The former requires 
light sfermions (as the cross section is oc (m^)"^), while the latter requires the sfermion to 
have a mass close to the one of the bino-like LSP. 

If the LSP has a non-negligible higgsino component, it can efficiently annihilate into 
a pair of gauge bosons (WW, ZZ) and doublet-like Higgs bosons (W^H"^, ZH, HA) 
via s-channel Z or Higgs exchange, as well as through t-channel neutralino and chargino 
exchange. For sufficiently heavy LSPs [rriy^o > rrit), annihilation into a pair of top quarks 
via Higgs exchange can significantly contribute. In the large tan /3 regime, where the Higgs 
couplings to bb and t~^t~ are enhanced, bb and t~^t~ final states are possible via s-channel 
Higgs exchange (not necessarily resonant). 

A tiny higgsino fraction already opens the possibility of annihilation through resonant 
s-channel Higgs exchange. However, to overcome the smallness of the Higgs-XiXi coupling 
in this gnificant fine tuning is required so that the corresponding Higgs boson has 

the appropriate mass ~ 2m^o, which is not always possible in the MSSM. 

Higgsino-like LSPs can also coannihilate with neutralinos and charginos. In the MSSM, 
it can happen that coannihilation is excessively efficient in this case, so that the relic 
density is too small. In the NMSSM, a non-vanishing singlino component can reduce the 
LSP couplings (and hence the coannihilation cross sections), leading to viable values for 
the relic density [501]. 

9.1.2 NMSSM-specific annihilation processes 

In the NMSSM, the richer scalar/pseudoscalar Higgs sector allows to have the correct relic 
density in large regions of the parameter space, from low to high tan/3 [346,501]. The 
possibility of light scalars and pseudoscalars consistent with LEP constraints also implies 
that new final states are kinematically open [346,501]. For example, nearly pure binos can 
annihilate via i^i-resonances into Ai^i. 

A distinctive feature of the NMSSM is the possibility of a nearly pure singlino-like LSP. 
However, due to the small couplings to SM particles, singlino-like LSPs tend to have a 
too large relic density due to comparatively small annihilation cross sections. Still, there 
are several possible mechanisms that allow to render a singlino-like LSP compatible with 
WMAP bounds [346,501]: 

The couplings of a singlino-like LSP (with N^^ ~ 1) to CP-even Higgs bosons Ha are 
approximately given by ~ —^/2Sa3K, (see (A. 14)), with Sas replaced by Pas for CP-odd 
Higgs bosons. Accordingly any CP-even (or CP-odd) Higgs state, whose singlet component 
Sa3 (or Pa3) is Sufficiently large, can give an important contribution to LSP annihilation 
via s-channel exchange, if k is not too small. 
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If the singlino-like LSP has a small higgsino component, resonant s-channel annihilation 
via scalar /pseudoscalar Higgs or Z bosons allows for viable NMSSM scenarios if m^o < 
50 GeV (occurring in regimes of k ^ A, and not very large /ieff), or even m^o ~ a few GeV 
provided that 2m^o ~ Mu^^Ma^- Annihilation can then proceed into XiXi ~^ bb or, for 
the case of Hi exchange and a hght pseudoscalar, into XiXi ~^ ^i^i- 

For a heavier mixed singlino/higgsino-like LSP {m^o > 100 GeV) and large A, the 
annihilation can occur via t-channel Xi exchange into pairs of mostly singlet-like Hi and 
Ai due to the enhanced XiXiHi{Ai) couplings ~ A. 

An LSP with a large singlino component can also coannihilate efficiently with heavier 
neutralinos (provided their mass difference is < 10 GeV), especially with a higgsino-like 
X2, leading to Xi + X2 ^ heavy doublet-like Higgs exchange. Such coannihilation 

cross sections can be further enhanced by s-channel resonances: if the heavy states H3 and 
A2 belong to the heavy Higgs doublet and are nearly degenerate, one may even obtain a 
"double resonance" with two s-poles from both H^ and A2 exchange. 

For a nearly pure singlino {Nf^ > 99%) in regimes of very small A and (^^ ^ A < 0.01), 
no efficient annihilation mechanisms are available, and even LSP-NLSP coannihilation (e.g. 
with heavier neutralinos/charginos, or sfermions) fails to diminish the LSP relic density 
sufficiently. In this case, the LSP density can be reduced by NLSP-NLSP annihilation 
(assisted coannihilation) in cases of, for instance, a bino-like NLSP or a sfermion NLSP. 
Assisted coannihilation is relevant for the fully constrained NMSSM [216,217] as discussed 
in Section 7.1, in which case it requires nearly degenerate Xi ^i-nd fi masses. 

However, to this end the reaction rate for the process LSP + X — NLSP + X' (~ 
nhsp'^x o'') has to be at least as large as the NLSP annihilation rate (~ ''^nlsp'^)) other- 
wise the LSP will no longer be in thermal equilibrium with the NLSP near the freeze-out 
temperature, but can be considered as decoupled. This condition allows to derive a rough 
lower bound on A: the process LSP + X — NLSP + X' is made possible only through the 
non-singlet component of the LSP; for a dominantly singlet-like LSP, the analysis of the 
neutralino mass matrix shows that this non-singlet component is proportional to A. Hence, 
the cross section a' is proportional to A^, and a' /a ~ \^ / g^- As mentioned above, the 
density nx of light SM particles is about nx ~ IO^tinlsp near the freeze-out temperature. 
Then, since tt-nlsp ~ '^lsp, the reaction rate for LSP + X NLSP + X' is larger than the 
NLSP annihilation rate only if \^ / gl > 10'^ or A > 10"^ [216,217]. 

9.2 Dark matter detection in the NMSSM 

In the NMSSM, the prospects for direct and indirect detection of neutralino DM can also be 
substantially different from what is expected in the MSSM (for a review of WIMP detection 
see [496]); subsequently we discuss direct and indirect detection separately. 

9.2.1 Direct dark matter detection 

Direct detection of neutralino dark matter (and generically WIMP dark matter) proceeds 
via the measurement of the recoil energy deposited by the scattering of WIMPs on nuclei 
in a detector. The energy transfer arising from the elastic scatterings typically lies around 
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0(10 keV), with predicted event rates below 0.1 events/kg/day. WIMP-nucleus interac- 
tions can be spin- independent (SI) or spin-dependent (SD). The SD contribution (via the 
axial-vector interaction) is only non-zero if the target nucleus has a net non- vanishing spin 
(unpaired nucleons), which is assumed to be carried by "odd-group" nucleons: protons or 
neutrons, whichever is most unpaired. The SI (scalar) contribution is proportional to the 
mass squared of the nucleus, and almost always dominates for nuclei with A > 30, like 
those used in most modern detectors (e. g. Germanium and Xenon). The total elastic cross 
section is the sum of both contributions; below we will focus on the SI contribution. To 
obtain the direct detection rate, one needs to take into account several astroparticle physics 
ingredients such as the WIMP density and the WIMP velocity distribution near the earth. 

With one remarkable exception, all dark matter direct detection experiments have only 
set upper limits on the WIMP-nucleon scattering cross section {ctn, N = n,p). Xenon [506] 
has reported a§ < 4.5 x 10^^ pb for a WIMP mass of 30 GeV {a^^ < 8.8 x 10"*^ pb for 
mwiMP ~ 100 GeV) at 90% CL, and CDMS [507] has set a limit a§ < 4.6 x 10"^ pb 
for a WIMP mass of 60 GeV at 90% CL. The best bounds on SD cross sections are from 
Xenon [508] ((t^° < 5 x 10"^ pb) and KIMS [509] (a^° < 0.18 pb). 

On the other hand, a positive result has been reported by the DAMA experiment (which 
looks for the so-called annual modulation signature): ctjv ~ 0.2 — 1.0 x 10~^ pb for WIMP 
masses between 30-100 GeV [510]. It has been noted in [321] that, in the NMSSM, very 
light bino- or singlino-like neutralinos (with m^o ~ 6 — 9 GeV) could allow to reconcile the 
DAMA results with the negative searches by other collaborations. 

In the absence of CP violation, and in the zero momentum transfer limit, the elastic 
scattering of a neutralino LSP with a quark in a nucleon can be described by the effective 
low energy four-fermion Lagrangian (in the Dirac fermion notation): 

jCcS = af XiXi m + «f X?757mX? Ql^l^q , (9.4) 

where ^^'^^ denotes the spin-independent (spin-dependent) interaction with up- and 
down- type quarks (described by spinors q = u,d). The contribution of the dominant SI 
term to the Xi-nucleon cross section is given by (see, e.g., [503]) 

4' = -I fl (9.5) 

where Jn/^n '^gC^q^/i^q, and each term in the sum over quarks must be weighted by 
the appropriate hadronic matrix elements. 

The scalar neutralino-nucleon interaction a^^ arises from s-channel squark and t-channel 
Higgs (or Z) exchange at the tree level, while neutralino-gluon interactions contribute at 
the one-loop level (where all (heavy) quark loops should be taken into account). The squark 
exchange term is identical to the MSSM, the only difference in the NMSSM arising if the 
LSP has an important singlino component implying a significant reduction of a^^'^. 

Significant differences in the NMSSM can arise from the t-channel Higgs contribu- 
tions [347,467,511-515], which can dominate the SI interaction provided the LSP higgsino 
component is large enough. One has (see, e.g., [514]) 

«f oc -j^ (i/aX?X?) (Haqq) (9.6) 

Ha 
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where the quark Yukawa couphngs (Haqq) and the Higgs-neutrahno-neutrahno Yukawa 
couphngs (HaXiXi) are given in Appendix A. The hght Higgs doublet contributions are often 
dominant, ahhough in the large tan/3 regime the contributions of the heavier doublet-like 
Higgs can also be important, as their couplings to neutralinos and quarks are enhanced [347]. 

For a singlino-like LSP, the cross section for the scattering off, e.g., a strange-quark is 
proportional to 

a§ ^{af'^f ^^S!,Sl. (9.7) 

Hence the exchange of a light Hi can lead to large direct detection cross sections within reach 
of the present generation of detectors [514,515]. (For a singlino-like LSP, the MSSM-like 
s-channel squark exchange diagram is in general too small to lead to cross sections close to 
the experimental limits. NMSSM specific contributions to the spin-dependent cross section 
can occur from a higgsino-like LSP, scattering off nuclei via Z boson exchange [347].) 

As an example, we display in Fig. 9 (from [515]) the SI cross sections in a large region of 
the unconstrained NMSSM parameter space, which is compatible with collider, cosmological 
and B physics constraints. (A heavy sfermion sector ~ 1 TeV is assumed, except for slepton 
masses of ~ 150 GeV to comply with {g — 2)^). Large detection cross sections can originate 
from the presence of light singlet-like Higgs (5^3 > 0.9) with a mass Mh^ ~ 50 GeV, see 
the right-hand panel of Fig. 9. Neutralinos within the reach of dark matter experiments 
are typically mixed singlino-higgsino states, with a mass between 50 — 130 GeV. The upper 
bound on the neutralino mass is due to the lightest stau becoming the LSP. If the slepton 
mass is increased, heavier neutralinos can be the LSP, but the resulting SUSY contribution 
to {g — 2)^ is generally too small. 

9.2.2 Indirect dark matter detection 

It is also possible to detect dark matter indirectly, looking for distinctive signals of WIMP 
annihilation products: gamma-ray fluxes, peculiar gamma-ray spectral features (like a sharp 
monochromatic peak at ~ m^o arising from the loop induced annihilation XiXi ~^ ll)^ 
sizable amounts of antimatter from pair annihilation in the galactic halo, and also energetic 
neutrino fluxes from annihilations in the center of the sun or earth. For a review in the 
MSSM we again refer to [496]. 

Indirect detection of NMSSM neutralino dark matter has been considered in [516]. Since 
the indirect detection rates generically scale as inverse powers of m^o (the pair annihilation 
rate being fixed by the requirement of a correct relic abundance), light neutralinos - as 
possible in the NMSSM - could give significantly enhanced rates. Additionally, new con- 
tributions (annihilations through a very light pseudoscalar) can enhance monochromatic 
Xi Xi 7 7 gamma-ray lines. The neutrino production rates from neutralino annihilation 
inside the earth and sun could also be distinct from MSSM expectations (see, e.g., [516] 
and references therein). 

NMSSM scenarios exist where the annihilation of sub-GeV neutralinos could explain the 
511 keV gamma-ray emission reported by INTEGRAL/SPI [517]. These gamma-rays can 
originate from low energy positrons possibly arising from the annihilation of an exception- 
ally light DM candidate. In the NMSSM, a very light neutralino ('tt.^o ~ 100 MeV) com- 
bined with a light pseudoscalar, could potentially produce the 511 keV gamma-rays [321]. 
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Figure 9: Scatter plots of the spin-independent neutralino-nucleon cross section in the general 
NMSSM for tan/3 = 5, with the remaining parameters in the ranges 0.01 < A, k < 0.7, 110 GeV < 
Ma < 430 GeV, 110 GeV < fi^s < 300 GeV, -800 GeV < Ax < 800 GeV, and -300 GeV < 
< 300 GeV. All the points fulfil all the experimental constraints and have a relic density in 
agreement with 0.1 < Qh'^ < 0.3 (astrophysical bound, grey dots) or the WMAP constraint (black 
dots). On the left as a function of the neutralino mass (GeV), and on the right as a function of 
the lightest scalar Higgs mass (GeV). From [515]. 



However, efficient LSP annihilation via resonant s-channel pseudoscalar exchange would 
require Ma^ ~ 2m^o ± 10 MeV. 

Recently, the PAMELA experiment [518,519] reported a cosmic ray positron excess, 
with a positron over electron fraction that appears to rise at energies from 10 to 100 GeV. 
However, the same detector does not report any obvious antiproton excess in the same 
energy range. In order to interpret such a positron excess (and the absence of a p excess) as 
dark matter annihilation products, it would be required that the DM particles dominantly 
annihilate into leptons (either due to dominant couplings to SM leptons, or annihilation into 
intermediate particles, sufficiently light so that they cannot decay into hadrons). Neutralino 
annihilation in the NMSSM (for m^o ~ 160 GeV, and a very light pseudoscalar) [329, 
520,521] has been put forward to explain the PAMELA excess: relatively light neutralinos 
(typically bino-higgsino mixtures) can annihilate as XiXi ~^ HiAi via A2 exchange, with Ai 
so light that only decays into a pair of leptons (muons, or even electrons, for M^^ < 2m^) are 
kinematically allowed. In this case, the lightest scalar Hi also dominantly decays into AiAi, 
so that the final state of the LSP annihilation indeed consists of very energetic leptons. 
Extensions of the MSSM by a singlet and an additional heavy lepton [522] (exhibiting two 
DM components) have also been proposed in order to accommodate the PAMELA data. 

9.3 Specific scenarios 

The (semi)-constrained NMSSM 
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The relic density in the cNMSSM was first discussed in [335,339,342]. In the meantime 
constraints on Higgs and sparticle masses as well on the dark matter relic density [343,344] 
have become tighter (and RGEs/radiative corrections are known to a higher accuracy), with 
the result that some of the regions in the parameter space of the cNMSSM considered in 
these early studies are no longer phenomenologically viable. On the other hand, a singlino- 
like LSP was believed to be excluded [339,342], which turned out to be incorrect after the 
inclusion of additional (co-)annihilation channels. 

In more recent studies, the parameter space of the semi-constrained NMSSM (relaxing 
the universality requirements for m| and at the GUT scale, see Section 7.1.3) has been 
thoroughly investigated. Regions allowed by constraints from collider and flavour physics 
and a WMAP compatible relic density have been identified in [346,347]. In the semi- 
constrained NMSSM, less options for LSP annihilation than in the general NMSSM are 
available: 

(i) MSSM-like mechanisms: annihilation of a bino- or bino-higgsino-like LSP through 
sfermion or Higgs boson exchange; coannihilation of a bino-like LSP with sfermions; 

(ii) annihilation near a pseudoscalar singlet-like Higgs resonance in the s-channel (A ~ 0.1, 

any tan /3); 

(iii) (assisted) coannihilation of a singlino-like LSP with a higgsino-like NLSP (A ^ 1, large 
uiQ and tan (3); 

(iv) (assisted) coannihilation of singlino-like LSPs with bino-like NLSPs, with the bino 
rapidly annihilating through a Higgs resonance (A <^ 1, large tan/3); 

(v) (assisted) coannihilation of a singlino-like LSP with a stau or stop NLSP (A ^ 1, small 
mo). 

In Fig. 10 (from [347]), we display some regions in the (mo, M1/2) parameter space, where 
a DM relic density within (or below) the WMAP bound is obtained through different anni- 
hilation mechanisms. (Also indicated are other relevant theoretical and phenomenological 
constraints.) The dominant annihilation mechanisms depend on (mo,Mi/2), and on the 
values of the remaining parameters A, Aq, A^^ and tan/3 of the semi-constrained NMSSM. 
On the left, these are chosen as A = 0.1, Aq = —900 GeV, = —60 GeV and tan/3 = 5, 
and the LSP is bino-like. Within the green horizontal band around M1/2 ~ 350 GeV, the 
correct relic density is achieved through a pseudoscalar resonance (ii). In the green region 
next to the blue region excluded by a stau LSP (at small mo), the dominant mechanism 
is bino-stau coannihilation (i), whereas for small mo and M1/2 ~ 250 GeV the dominant 
mechanism is bino-stop coannihilation. On the right we show an example with small A and 
a singlino-like LSP. In the small mo regime efficient singlino annihilation is achieved via 
coannihilation with a nearly degenerate lightest stau (v). 

The prospects for direct DM detection have also been considered in [347]. Present 
bounds already constrain such regions at large tan /3 where the heavy Higgs doublet is 
relatively light and/or its couplings to an LSP with a non- vanishing higgsino component 
and quarks are enhanced. Large cross sections near or above present limits can also be 
obtained for A ~ 0.1, where the LSP annihilates via a pseudoscalar singlet resonance 
(ii), and for bino-stau coannihilation regions (i). Hence, the semi-constrained NMSSM 
parameter space with a dominantly bino- or higgsino-like LSP has good prospects of being 
probed by the future direct detection experiments. On the other hand, for small tan /3 or 
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Figure 10: Constraints in the (mo,Mi/2) plane in the semi-constrained NMSSM: regions 
excluded by theoretical constraints or by LEP/Tevatron searches on sparticles are indicated 
in black, by charged sfermion LSP in blue, from violation of LEP limits on Higgs searches 
in red, by constraints from B-physics in pink and disfavoured by {g — 2)^ in violet. The 
regions allowed by the upper bound on the DM relic density are displayed in green. On 
the left: A = 0.1, Aq = -900 GeV, = -60 GeV and tan/3 = 5. On the right: A = 0.01, 
Ao = 300 GeV, A^ = 50 GeV and tan/3 = 10. Here, in addition, the regions with a 
singlino LSP are shown in grey/cyan (with a relic density above/below the WMAP limit). 
From [347]. 

TeV-range heavy Higgs and squarks, large scale detectors will be required. The detection 
of a singlino LSP is beyond the reach of (even far future) large scale detectors. 

The cNMSSM where = Aq is imposed at the GUT scale, but ms is allowed to differ 
from mo, has been investigated in [348] (see Section 7.1.3). It is still possible to have an LSP 
with a large bino or higgsino component which satisfies WMAP constraints via neutralino- 
stop (and stau) coannihilation or via Higgs resonances. For most points in this parameter 
space, direct detection of DM would appear possible for upgraded large scale detectors. For 
larger A, a singlino-like LSP can annihilate via a CP-odd Higgs resonance in the s-channel. 

In the fully constrained NMSSM [216,217], discussed in Section 7.1, the LSP is always 
singlino-like, and A < 10^^. The only available mechanism to reduce the relic density is 
"assisted coannihilation" (v) with the stau NLSP (provided A > 10"^), which requires a 
near degeneracy of the Xi ^i-nd fi masses. Since direct (or indirect) detection of a singlino 
LSP relies on its non-singlet components of 0{X), the prospects to detect the LSP in the 
cNMSSM are quite dim: WIMP-nucleon cross sections are extremely small, and indirect 
detection of the products of the annihilation process XiXi ~^ (with a pseudoscalar 
sufficiently light to be produced on-shell) also appears impossible, as the AiXiXi coupling 
is proportional to k which is also tiny in the cNMSSM (k < 10"'^). Hence, the fully 
constrained NMSSM can be excluded by the direct or indirect detection of a WIMP-like 
dark matter candidate. 
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The nMSSM 

The neutrahno rehc density has also been studied in the framework of the nMSSM (see 
Section 7.3). In [176,379] it has been shown that, for large A, one can simultaneously 
achieve the correct VLh"^ and a strong first order phase transition leading to successful EW 
baryogenesis. This occurs for a very light LSP with a mass m^o ~ 30 — 40 GeV, mostly 
singlino-like but with a sizable higgsino component. The dominant annihilation modes are 
via (possibly resonant) s-channel Z, Hi and Ai exchange [176,378-380], the latter being 
favoured by the constraint from [g — 2)^ [380]. 

Regarding DM detection [379], the dominant contributions to SI and SD interactions 
arise typically from Higgs and Z t-channel exchange, respectively. Due to the smallness of 
the XiXi^ coupling, the SD cross section lies beyond the reach of next generation of DM 
experiments. On the other hand, the SI cross section could be within reach of the current 
and next generation of DM experiments for sizable values of A (leading to an enhancement 
of the XiXi^a coupling) and a for comparatively light Higgs spectrum. 

U{1)' extensions 

In U{iy extensions (see Section 8.2.1), the inclusion of an additional f/(l)' gaugino 
(bino') and the absence of the singlet self-coupling term in the superpotential can signifi- 
cantly modify the nature and properties of neutralino dark matter relative to the MSSM 
and NMSSM [378,463-468]. Just like for the nMSSM, the absence of the term in the 
superpotential restricts the annihilation modes of the singlino. Compared to the NMSSM, 
new annihilation channels include t-channel x^ mediated annihilation into a pair of Z's, 
and Z' (resonant) s-channel exchange. The analysis of [468] shows that, depending on 
the nature of the LSP (higgsino, mixed singlino- higgsino-bino', singlino) and on the re- 
maining spectrum, the correct DM relic density can be obtained via s-channel resonances 
{H2/ A/ Z' / H^). Pure singlino neutralinos would not directly couple to the singlet Higgs 
boson, but they couple strongly to the Z' so that even very heavy LSPs can efficiently 
annihilate. The prospects for direct detection of neutralino DM in t/(l)' extensions have 
been addressed in [463,468]. 

Dark matter in supersymmetric U{\y models with a secluded f/(l)' breaking sector 
(sMSSM) has also been studied [464]. For a light singlino-higgsino LSP with a mass 
< 100 GeV and an enhanced coupling to the Z, annihilation through the Z-resonance 
leads to a viable relic density over a large part of the parameter space without violating 
LEP constraints. 

A comparative study of unconstrained versions of the MSSM, NMSSM, nMSSM, t/(l)' 
extended models (UMSSM) and secluded models (sMSSM) has been carried out in [378,467], 
identifying allowed ranges for the mass of the LSP in each case as displayed in Fig. 11. 
Prospects for direct detection of neutralino DM were also compared in [467]: the cross 
sections for the U{1)' extended models are similar to those of the NMSSM, and can be 
within reach of SuperCDMS and WARP. In view of the (comparatively) lighter LSPs in 
secluded models as well as in the nMSSM (see Fig. 11), spin-independent proton scattering 
cross sections could be detectable at CDMS. 

In the exceptional SSM (EgSSM), which in addition to the Z' comprises three families 
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Figure 11: Neutralino relic density versus the lightest neutralino mass for the MSSM and 
NMSSM-like extensions, with 0.099 < Qh"^ < 0.123 within the blue lines. Too large annihi- 
lation cross sections (leading to a too small Qh"^) can originate from the Higgs boson pole 
in the s-channel in the MSSM, NMSSM and UMSSM for m^o ~ MhJ2 ~ 60 GeV, and 
from the Z boson pole at m^o ~ Mz/2 in the n/sMSSM. (From [467].) 



of Higgs doublets and singlets, it has been shown in [491] that the two families of "inert" 
doublet higgsinos and singlinos (i.e. whose scalar partners are vev-less) lead to a decoupled 
neutralino sector with a naturally light LSP, that successfully accounts for Qh"^ provided 
that tan/3 < 2. 

As referred to in Section 8.2.1, non-neutralino DM candidates have also been put for- 
ward in the framework of 11(1)' extensions: exotic LSPs [470], and right-handed sneutri- 
nos [469,475]. The correct relic density of right-handed sneutrinos can be obtained via 
t-channel Z' exchange in large regions of the parameter space, and provided the Z' mass is 
not excessively large, s-channel resonant Z' exchange also provides an efficient annihilation 
mechanism. The prospects for the detection of these DM candidates have been investi- 
gated [469]. U{1)' extended models with sneutrinos as the LSP as well as the NLSP [475] 
have been proposed to explain the PAMELA observations, while accounting for a relic 
density in agreement with the WMAP bounds. 

Other variants 
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In large-A models [128] (see Section 3.2), the LSP is higgsino-like in most of the phe- 
nomenologically viable regions of parameter space. Then, the relic density is in general 
too small; only for tan/3 < 1.7 regions with a substantial DM abundance can be found: 
annihilation via s-channel Z exchange can be reduced by increasing the singlino component 
of the LSP, or via a partial cancellation of the higgsino-Z coupling. 

In the NMSSM including right-handed neutrino superfields with couplings to the sin- 
glet Higgs (thus providing a dynamically generated non-vanishing Majorana neutrino mass 
term), the properties of sneutrino dark matter have been discussed in [436]. For natu- 
ral values of the input parameters, right-handed sneutrinos with a mass in the range of 
5 — 200 GeV can reproduce the observed dark matter relic density without being excluded 
by direct dark matter searches. In this case, the predicted direct detection cross sections 
are generally within the reach of future experiments. 



87 



10 Conclusions and outlook 



In this paper we reviewed the present status of knowledge of the theoretical and phe- 
nomenological aspects of the NMSSM, including physics at colliders, low energy precision 
observables and the properties of the dark matter candidate, focusing on the possible dif- 
ferences with respect to the MSSM. 

We hope that the direct production of sparticles and/or Higgs bosons at the Teva- 
tron/LHC will provide us information about possible supersymmetric extensions of the SM 
in the near future. Subsequently it will be essential to study the precise nature of the super- 
symmetric extension of the SM as thoroughly as possible. Depending on the experimental 
results, it may not be easy to distinguish the MSSM from non-minimal extensions as the 
NMSSM or variants thereof, but possibly the properties of the Higgs sector and/or the 
neutralino sector (in particular the nature of the LSP) will signal the presence of an addi- 
tional singlet superfield. Clearly, these signals must be well understood in order to be able 
to interpret the data correctly. In view of the various scenarios that are possible within 
the NMSSM (as unconventional Higgs decays, and/or a neutralino with a large singlino 
component), additional efforts in the form of simulations of events corresponding to such 
scenarios are still required. 

On the other hand, the properties of the Higgs and/or the neutralino sectors could 
correspond to those expected within the MSSM - which would not prove, however, that 
the MSSM is the correct supersymmetric extension of the SM at the TeV scale: at least 
in the decoupling limit (see Section 2.2.1) - and without a singlino-like LSP - the physics 
of the NMSSM becomes indistinguishable from the one of the MSSM. The question, for 
which values of parameters (like A) the two models can be distinguished experimentally, is 
complicated due to the large number of unknown parameters, and the answer depends on 
the nature and the precision of future experimental data. Also here, additional studies are 
desirable (at least once data is available). 

In any case it will be essential to combine the various pieces of information from low 
energy precision observables, physics at colliders and the properties of dark matter; ideally, 
these will allow us to pin down the fundamental theory at the TeV scale without too many 
ambiguities. 
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Appendix A. Mixing matrices and tree level Higgs cou- 
plings 

The Feynman rules of the NMSSM have been described first in [50] (including the quartic 
couplings not given here), and subsequently in [52]. 



A.l Mixing matrices 

First we define the mixing matrices which diagonalize the Higgs, neutralino and chargino 
mass matrices (after the addition of radiative corrections). 

The CP-even 3x3 mass matrix (2.22) in the basis i/f = {HdR, H^r, Sr) (in the 
SLHA2 conventions [49]) is diagonalized by an orthogonal 3x3 matrix Sij rotating the 
basis H^'"''^, 

such that the mass eigenstates Hf^^^^ are ordered in increasing mass. 

The CP-odd 2x2 mass matrix Mp (2.24) in the basis {A = cos /3 + sin (3 Hai, Sj) 
was obtained after the rotation by the angle /3 (2.25), which allowed to omit the Nambu- 
Goldstone boson. It can be diagonalized by an orthogonal 2x2 matrix such that the 
physical CP-odd states Af^^ (ordered in mass) are 

^mass _ P[^A + P[^Sl, 

(A.2) 

In terms of the weak eigenstates Af^^^ = {Hdi, Hui, Sj), the mass eigenstates {A"'^^^, A^^^) 
are given by 

where 

Pii = sin/3i^i , P,2 = cos{3Pi, P^3 = Pl2■ (A.4) 
The inverse relation reads explicitly (omitting the Nambu-Goldstone boson) 

ZJ TD A mass 1 d a mass 

ndl = -fll/i^ + 121^2 ) 

ZJ TD A mass 1 d A mass 

n-ul = -fl2^l + -f22^2 5 

Si = Pi3A5"^^^ + P23Ar'' ■ (A.5) 

Omitting again the Goldstone boson, the charged weak eigenstates contain a com- 
ponent of the physical charged Higgs boson given by 

= cos (3H^ , = sin /3 . (A.6) 

The (symmetric) 5x5 neutralino mass matrix Ai^ (2.32) in the basis 

= {—iXi,—iXl,ip^,xl:'l,xl:s) is diagonalized by an orthogonal real matrix Nij, such that 
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the physical masses m^o ordered in |m^o| are real, but not necessarily positive. Denoting 
the 5 eigenstates by we have 

X° = iV.,^° . (A.7) 

The diagonalization of the (not symmetric) 2x2 chargino mass matrix X (2.35) in the 
basis ip~ , if)^ (see (2.33), (2.34)) requires different rotations of il)' and into the mass 
eigenstates x"? X"*" 

X- = Ur, X^ = V^^ (A.8) 

with 

cos6';7 sin6'{/ \ / con 6y sin6'y 



U=[ , V=[ . (A.9) 

\^ — sm Uu cos ojj J y — sm Uy cos uy J 

Finally the top squark mass matrix (2.36) is diagonalized in terms of the mass eigenstates 
ti, t2 (with < m^^), the bottom squark mass matrix (2.37) in terms of the mass 
eigenstates bi, 62, and the tau slepton mass matrix (2.38) in terms of the mass eigenstates 
^i! ^2 by the rotations 

ii = cos 6't + sin 6*7- tfl , ^2 = cos 6'^ t/? — sin 6't , 

bi = cos 9b b^ + sin 9b bji , ^2 = cos 9b bji — sin 9b b^ , (A. 10) 

Ti = cos 9t tl + sin 9^- tr, = cos 9^ r/j — sin 9^ . 

Now we proceed to give the Higgs couplings (Feynman rules) in terms of physical mass 
eigenstates. 



A. 2 Higgs couplings 

Higgs-quarks/leptons 

We consider Higgs couplings to quarks and leptons of the third generation only, from which 
the remaining couplings can be deduced easily. First we recall the relations between the 
quark and lepton masses rrit, mb, rrir and the Yukawa couplings ht, hb and h^-, where 



V 



2 



vj (related to by (2.13)): 



rrit nib nir 



V sin f3 V cos f3 v cos f3 



(A.ll) 



(Corrections to hb are discussed in Section 3.3.) Below, the left and right-handed top 
(bottom) quark and lepton bi-spinors are denoted by ti, {pi, bu) and i/^-^, tl, tr, and 
the quark couplings are diagonal in the colour indices. (For simplicity, the index '^'^^^ of the 
CP-even mass eigenstates H^^^^ and the CP-odd mass eigenstates A^^'^'^ will be omitted 
subsequently.) The denominators \/2 in the couplings originate from the rescaling of the 
kinetic terms of real Higgs bosons with respect to complex Higgs bosons. 
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HiTLT'k 

H-tLb% 



ht ^ 

"7! " 

hb ^ 

71 " 



5,; 



il 



■ hb „ 

■ hr p 

ht cos /3 
— /ife sin /3 
— /it- sin /3 



(A.12) 



Higgs-gauge bosons 



H,{p)H+{p')W- 
A{p)H^{p')W- 

H,{p)A,{p')Z^ 
H^{p)H-{p')Z, 



q'^ g2 

: gf,^^—^{vdSii + VuSi2) 



V2 



9l 



g^y—={vdSii + VuSi2) 
v2 
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COS /3S'i2 - sin (5Sii){p - p') 



i — {cos 13 Pi2 + sin/3Pii)(p - p% 



9l - gl 



{p-p)^ 



(A. 13) 



Higgs-neutralinos / charginos 

Here, the couplings to two neutralinos involve the combinations n^j* = NiaNj^ + Nn,Nja of 
the neutralino mixing matrices: 
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(^-^P^,U,2V,2 - ^{Pa2UilV,2 + PalU^aV.l)^ 



Sill [j 

H^Xix] ■■ A cos pUi2Nj^ - —j=-U,2{9iNji + g2N^2) + 92 sin ^UiiN,^ 

cos/3 

A 



COS 

H'xtx] ■■ A sin PV,2N,5 + —fv^2{9iNji + 92N,2) + 92 cos (5V,iNji 

V2 



j5 

-^{PaJi^ - Po^) + |(i^in^ - i^^nj)) (A.14) 



Triple Higgs couplings 



A^A-p.123 I tt333 



where 



2 2 

+^ (-.(n^L^ - n-) - vMl - n^L^)) 

+^ i'^^'Il - >^Kll) (A.15) 



■^^oftc — SaiSbjSck + SaiScjShk + SfjiSajSck 

''rSbiScjSak + S^iSajShk + SciShjSak ■ (A. 16) 
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where 



where 



313\ 



333 
abc 



-Llr2±|'nl23 , TT213 , TT312N _ '^^'^ -rT333 
2 2 

+^ (-.(n^L^ - n-) - tan- - n-)) 

+^ (A {liTc - ^Tc - ^T) + ^Kf^) (A.17) 



^abc ~ Sai{PbjPck + PcjPbk) ■ (A. 18) 



+^nf^ (A.19) 



= 2SaiCjCk with Ci = cos /3 , C2 = sin /3 . (A.20) 



In contrast to the previous couphngs, the triple Higgs couphngs in the general NMSSM 
differ slightly (by the terms ~ fi') from the ones in the scale invariant NMSSM. The quartic 
Higgs (and Higgs-gauge) couplings, which are less relevant for the Higgs phenomenology 
(i.e. Higgs-to-Higgs decays) can be found in [50]. 



Higgs-squarks / sleptons 

Here we assume diagonal Yukawa couplings in family space which we denote as hua, hda, 
hea, a = 1,2, 3, with hu3 = ht, hd3 = hb, hes = hr. The soft trilinear couplings are Aua with 
= At, Ada with = Ab, and A^a with Aez = Ar. 

The squarks are denoted as Uia, Dia, Ujia, Dua, the charged sleptons as ELa, Ejia, and 
the sneutrinos as uia- 
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g{H+,ULa,DLa 
g{H+,ULa,DRa 
g{H+,URa,DLa 

g{H+,URa.DRa 



9l 


9l 


12 


4 






9l 




12 





In the case of the first two famihes (a = 1,2), these states can be considered as mass 
eigenstates. For the third family (apart from the tau sneutrino), the couphngs of the mass 
eigenstates defined in (A. 10) depend on the corresponding squark and slepton mixing angles. 
We denote the couplings of the squark/slepton weak eigenstates (= mass eigenstates for 
a = 1, 2) by g{- ■ ■) for all three families; the couplings of the mass eigenstates of the third 
family are subsequently expressed in terms oi g{. . .) with a = 3. 

Higgs-squarks: 

g{Hi, ULa, Ula 
g{Hi,DLa,DLa 
g{Hi, URa, URa 
9iHi, DRa, DRa 

g{Hi, Ulo, URa 

g{Hi, Dia-, DRa 
g{Ai, Ula, URa 
g{Ai,DLa,DRa 



V2(} 



h'i„v„,Sio - ^ (u„S, 



VdS; 



( hlaVdSii + Y i^uSi2 - VdSil) 



-h. 



V2 



{AuaSi2 — ficsSil — XVdSis) 



(AdaSil — ^esSi2 



V2 

hda 

VuVd, 2 



{A^aPi2 + fJ'csPil + >^VdPi3) 



{AdaPil + fic{iPi2 + Al^n-Pis) 
"-ua "'da) 1 



-hdaicos /3 /i + sin /3 Ada) 
-/i„a(sin/3 /i + cos/3 A^a, 

hua^da 



(A.21) 
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Third family squarks, CP -GVGIl SCctlcirS Si'. 

Hiiiii : cos 9^ g{Hi, Uls, Uls) + sin 6^ g{Hi, Urs, Um) 

+2 cos 6t sin Or g{Hi, Ul3, Um) 
Hii2i2 : sm9^g{Hi,UL3,UL3) + cos9^g{Hi,UR3,Uii3) 

-2 cos 9t sin Or g{Hi, Ul3, Urs) 
Hiiii2 : cos Or sin Ot {g{Hi, Ur3, Ur3) - g{Hi, Ul3, f^Ls)) 

+ (C0S^2 _ gjj^^2) g(^^,^ u^^^ jj^^^ 

Hibibi : cos9%g{Hi,DL3,DL3) + sm9%g{Hi,DR3,DR3) 

+2 cos 9b sin 9Bg{Hi, Dl3, Dr3) 
Hib2b2 : sin 9l g{Hi, Dl^, Dl^) + cos 9l g{Hi, Dr^, D^) 

-2 cos 9b sin 9Bg{Hi, Dl3, Dr3) 
Hibib2 : cos 9b sin 9b {g{Hi, Dr^, Drs) - g{Hi, Dl3, Dl3)) 

+ (cos 9l - sin 9l) g{H,, Dl3, Dr,) (A.22) 

Third family squarks, CP-odd scalars Af. 

AMi ■■ 

AM2 : 

Aiiii2 : g{Ai,UL3.UR3) 

AAbi ■ 

Aib2b2 : 

AAb2 : (?(v4i,Di3,^iJ3) (A.23) 
Third family squarks, charged scalar H^: 

H^hbi : cos9TCos9Bg{H+,UL3,DL3) + sm9Tsm9Bg{H+,UR3,Dn3) 

+ cos 9t sin 6^5 Uls, Drs) + sin cos 9b g{H^ , U^s, ^ls) , 

H+iih : -cos9Tsm9Bg{H+,UL3,DL3) +sm9TCos9Bg{H+,UR3,Dn3) 
+ cos 9t cos 9b g{H^,UL3,DR3) - sm9Tsm9B giH'^,UB3, Dls) , 

i7+t2^i : - sin 9t cos 9b g{H+, Ul3, Dls) + cos 9t sin 9b g{H+, Ur^^, Dr^) 

- sin 9t sin 9Bg{H^, Ul3, Drs) + cos 9t cos 9b g{H^, Ur3, Dls) , 
H+kh : sm9Tsm9BgiH+,UL3,DL3) + cos9TCos9BgiH+,UR3,DR3) 

- sin 9t cos 9b g{H+ , Ul3, Drs) - cos 9t sin 9b g{H+ , Ur3, Z^ls) • (A.24) 
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Higgs-sleptons: 

g{Hi, Ena, Ebo) = V2 (^hl^VdSn + y {vuSi2 - VdSn)^ , 



g{Hi,iyLa,l'La) 
gi^i, Eia, Ejia) 
g{H'^,ELa,l^La) 



7! 



2VuVd 

92 — 



(A.25) 



Third family sleptons (staus), CP-even scalars Hi (the couphng HiUis^Ls does not differ 
from (A.25)): 

HiTiTi : cos 9l g{Hi, El3, Els) + sin 6l g{Hi, Er^, Ers) 

+2 cos Ol sin Ol g{Hi, Els, Ers) 
HiT2T2 : sin 9lg{Hi, Els , ^L3 ) + cos 9lg{Hi, Er^ , Er^ ) 

-2 cos 9l sin 9l g{Hi, Els, Ers) 
HiTiT2 : cos6'i sin6'L {g{Hi, Er^, Er-^) - g{Hi, Els, Els)) 

+ {cos9l-sm9l)giH,,ELs,ER,) 



(A.26) 



Third family sleptons (staus), CP-odd scalars Af. 



AiTiTi 

AiT2T2 
AiTiT2 








g{Ai, Els, Ers) 
Third family sleptons (staus), charged scalar H~^: 



(A.27) 



H'^tiUls ■ cos9L^^^{gl - hj^) - sm9Lhes{cos (3 jj, + sin/3 A^s) 
H'^r2VLs ■ - sin 9 L^^^ {gl - hl^) - cos 9 lKs{cos (3 iJ, + sin P Aes) 



(A.28) 



Yukawa induced radiative corrections to triple Higgs couplings 

As in the MSSM [523-528], the dominant radiative corrections to triple Higgs couplings 
originate from top/bottom-quark loops, which are regularized in the ultraviolet by the 
corresponding contributions from squark loops. Here we give these radiative corrections in 
the NMSSM, which just require to generalize the Higgs mixing matrices of the MSSM to 
the NMSSM, see Appendix A.l. 
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However, we confine ourselves to the contributions involving a potentially large loga- 
rithm: we assume approximately degenerate squark masses given by Msusy- The infrared 
cutoff inside the dominant logarithm is given by the mass rrig of the quark inside the loop 
or the mass of the decaying Higgs boson through this vertex (whichever is larger); we will 
denote it by Mh here. Then the potentially large logarithm is 

t ^ In (MIusy/ max(mj, M^)) . (A.29) 

The dominant corrections to the triple CP-even Higgs interactions (A. 15) are then given 

by 

AH^H,H^ : {htv^n'f^ + htv.Ulll) (A.30) 

with H^^j^ as in (A. 16). The dominant corrections to the CP-even/CP-odd Higgs interactions 
(A. 17) read the same, 

AH^A,A, : {htvuUlll + htv.Ulll) (A.31) 

where now H^^j^^ is given in (A. 18). The dominant corrections to the CP-even/charged Higgs 
couplings (A. 19) read 



3t 



8V2 7r-L 

+hlv4hlUl^^ + h^,{Uf^ + Ul''))] (A.32) 



with n^'' as in (A. 20). 
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Appendix B. Renormalisation group equations 



In this Appendix we give the renormahsation group equations in the DR scheme, vahd above 
the SUSY breaking scale, for the parameters of the general NMSSM. The conventions are 
t = InQ^, and the U{1)y gauge coupling gf is defined in the Standard Model normalisation 
(related to the GUT or SU{5) normalisation hj gf = ^ (gf^^) )• 

B.l Gauge and Yukawa couplings 

The RGBs are known to two-loop order [529-531] (for the NMSSM see [72,73]); it suffices, 
however, to include the Yukawa couplings of the third family only. (The CKM mixing 
matrix is not considered here.) 



167r 



2^ 

dt 



dt 



'.dgi 
dt 
,dh^ 
Itt 



167r 



167r' 



> dhl 

HI 



dt 



g\ /199 2 , ' 
= h't{Qh^t+hl + X'-^-^gl-3gl-^-lgl 



2 26 2 

■9s - -jK - y-. 



^-^hl - Qhl - 2A^) . 



2Agl - Qhl - Qhl - 2hl - 2\^ 



9l 



Ugl - Ahl 



' ^hl 



5ht - 3A' 



hl\'-2\\' + 2glhl+ ^, 

2743 4 _ 
79i + y^2 



5hlhl - 3hl\^ - hlhl - AhlX^ 

'^gfhl + Qglhl + IQglhl 
5 „ „ 13fi „ „ 



+ — ^ 
167r2 

- hlX 

X5 A 16 A 5 9 9 136 9 9 9 9 \ 

y^2 - y^3 + 3^i^2 + y^^i^s + §^2^3 1 > 



162 



7 2 _ o 2 _ 1^ 2 ' 



= hl(^6hl + hl + hl + X^-^g! 

+ f - 22ht - 5ht - 3ht - 3A^ - 5hlhl - Shlhl - 3hl\^ 

167r^ \ 

2A2«:2 + Iglhl + ^ 2(^2^2 ^ g^2^2 ^ ^g^2^2 



AhlX'' 

5 22, ^0 22io22l 
3^1^2 + ^^1^3 + 8^2^3 ) ' 



1435 4 
V9i 



162 

= /i2(^4/i2 + 3/i; 

167r2 V 



15 4 16 4 5 2 2 40 
y^2 - y^3 + 3 ^1^2 + ^5 



,2 , - -igl - 3g: 



3hl\'' 



- 2\^k'' 



|L _ 10/^4 _ - 3A^ - 'dhlhl - 3hlX^ 

75 ,4 , ^ 



- 2^?/.^ - Iglhl + 6^2^/^^ + 16glhl + y^7^ + y 
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IGtt^^ = X^(3h',+3hl + hl + 4X' + 2K'-gl-3gl 

+ ^ ( - lOA^ - 9ht - 9ht - 3ht - 8/€^ - 9\^h^ - 9\^hl 



167r2 V 



- ?,\'hl - 12X^k' - Qhlhl + 2glX^ + ^gfh', - \glhl + 2glhl 
+ QglX^ + IQglhl + IQglhl + '^gf + ^gl + Sgjgi 



dt \ J 167r2 

- 18/i?A^ - IShlX' - 6hlX^ + QglX^ + IS^^^A^J . (B.l) 

B.2 Gaugino masses 

The RGEs to two-loop order can be found in [532,533]. 

26 14 \ 

+ -hl{At - Ml) + -hl{A, - Ml) + Qhl{Ar - Ml) + 2X\A^ - Mi) J , 

167r2^ = (72'^2 + yII^ {^gl (Ml + M2) + 50(?2M2 + 24(72 (M2 + M3) 

+ 6/1? (Aj - M2) + Qhl{A^ - M2) + 2hl{A, - M2) + 2A2(Aa - M2)) , 

IGtt^^ = -3^7|M3 + ^ (y^7?(Mi + M3) + 9gl{M^ + M3) + 28^?32M3 

+ Ahl{A-M^)+Ahl{A,-M^)^ . (B.2) 

B.3 Trilinear couplings 

The two-loop /3-functions are known for the parameters of a general softly broken SUSY 
theory [94] . We concentrate on the trilinear couplings involving the squarks and sleptons of 
the third family; those of the first two families do not play an important phenomenological 
role, except for the muon trilinear coupling A^ (= Ae2) which enters the muon anomalous 
magnetic momentum {g — 2)^. 



99 



levr^ 



.dAt 
dt 



,dAb 



167r' 



,dA^ 
'dt 



167r' 



) dAi^ 

IF 



+ 



QhlAt + hlA, + A^A, + ^glM, + + ^^^^Ma 

^ - 44/1^ A - - 6AM, - hhlhl (A + A) 

3/i2A2(Ai + A,) - hlhl{A, + A,) - 4/i2a2(A + A^) 
hl\^ [Ar + Ax) - IX^h" (Aa + A«) + 2^72/,2 (^^ _ Ml) 

^(72/,2(^^ _ Ml) + ^glh\{At - M2) + 16^?|/i,2(A - Ms) 
2743 32 5 



136 



glgl{Mi + M3) - 8glgl{M2 + M3) 



7 16 

6hlAt + hlAt + /i^A, + A^Aa + -(??Mi + ?,glM2 + y (^s'Mg 

^ - 44/1^, - IQhtAt - Gh^Ar - 6AMa - 5hlh^, {Af, + At) 

- 3hlhl{A, + A.) - 3hlX^{At + Ax) - ^hl\\At + Ax) 

- 2A\2 {Ax + A«) + (A - Ml) + (A - Ml) 
+ 2(72/^2 (^^ - Ml) + 6(?2'/i' (^6 - M2) + 16(73^/1^ (A, - M3) 

- ^-^g\M, - 15glM2 + fglM, - ^g^gl {M, + M,) 



40 
27 



(7?^?32(Mi + M3) - 8^72^73(^2 + M3) 



4hlAr + 3/igA + X^Ax + 3^?Mi + 3^2'M2 



167r2 



- 20/i^ A - 18/i^ A - 6AMa - 9hlhl {Ar + Ah) - 3hlX^ (A + Ax) 

- 3hlhl {At + A) - 3/i?A2 (A + Ax) - 2X^K^ {Ax + A) 

+ 2^7i2/i2 (A - Ml) - yihl{A, - Ml) + Gglhl{A^ - M^) 

+ IGglhl {Ah - M3) - 75^7^ Ml - l^glM^ - 3glgl (Mi + M2) 

1 



3hiAh + KA, + AMa + 3^fMi + 3^2M2 + 



167r2 



- \8hlAh - Gh\Ar - 6AMa - 3h\hl {At + Ah) - 3h\X^ {At + A^ 

- 2A\2 (Aa + A) - ^ (7?/^' (A - Ml) + 2g\hl {A, - Mi) 

+ 16glhl{Ah - Ms) - 75gtM^ - IdglM^ - 3glgl{M^ + M2) V 
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167r^ 



dt 



16r 



<dA^ 
dt 



4AMa + ShfAt + ShiAb + KAr + 2k' A^ + Mi + ?,giM2 
1 



- QA^/l? ^ _ 9^2;^2 ^ _ 3^2^2 ^ 

- 12A2k2(Aa + A,) - 6/i?/i2(Ai + ^) + 2glX\Ax - Mi) 
+ \glh]{A, - Ml) - ^^?/.^(^ - Ml) + 2glhl{A^ - M,) 
+ QglX' {A^ - M2) + IQglhl (A - M3) + le^^^/z^ (A,. - M3) 

- 23^^Mi - IS^^Ms - ?,glgl{Mi + M2) 

1 



Qk'A^ + 6AMa + 



167r2 



- 24k2 A^ [A^ + Aa) - 18/i?A2 {At + Aa) - 18/i^A2 {A, + Aa) 

- 6/iiA2(A. + Aa) + 6(7^A2(Aa - Ml) + l^gl\'{A^ - Ivh)^ . 



(B.3) 



B.4 Squark and slepton masses 

Let us define the following quantities: 



M/ 
M^ 
M2 
M^ 



+ 



m: 



m: 



2 



2 I 2 

,2 , ^2 



+ A]^ 
+ Ah 



A- 



A ) 



3m^ + Af. 



Q 



2m^ + m|, - mi + m|] + 



Tr[: 

^? ( - "^Qs + ^'^c^s - 3mi^J +hl{-ml^- 2ml^ + 3mj^J 



2 2 



/ii(mi3 - 2m|3 + m^J + X'{ml^ 



m 



3 



16 



-m^ 



-m 



(7 



-m 



D 



-mi + 2m, 



-^2(^Tr[mJ - mi] + mj,^ - m\ 



-m^ 

3 ^ 



3 ^ 



2 2 

3 ^ 



2mi 



m^ 



1 . 
+ -(m 

2m 



— m 



(7 



m 



D\ ' 



gl [Ti [3m| 
^732Tr[2m^ + 



m 



mr 



m 



(B.4) 



where m denote matrices in family space. The two-loop RGEs then read 
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levr^^ = 5a,hlM^ + 5a,hlMl--glMl-^lMl--glMl + -gli 

106asht{M^ + A') - lOSa3ht{M^ + Al) 



- + m2 + 2AA^) - 6a3hlhl [Ml + Ml + 2A,A.) 

- (5,3/i^'(M2 + Ml + + ^5azglhl{Ml - 2MMt - Mi)) 

+ l5,,glhl{Ml - 2Mi(A - Ml)) + ^^^M,^ + - ^-^gtMl 

1 16 
+ -glgl{Ml + m2 + M1M2) + -glgl{Ml + M| + M1M3) 

+ 16(72(^32(M2 + M| + M2M3) + + + ^gla2 + ^^?3'^3 



1 



167r2 



mazht{Ml + A?) - 25aXhl{Ml + + 2 A A) 



- 25aihl\\Ml + M2 + 2AiA;,) - haiglhl{Ml - 2MMt - M,)) 

1712 64 
+ 65,3^?2'^'(M,^ - 2M2(A - M2)) + -^g{Ml - ^-^glMl 

+ ^f^^^^M^ + Ml + M1M3) - ^gli' + + ^(73^3 

+ Teb ( ~ ^^^'^^/^^(M^ + Al) - 25a,hlhl{Ml + Ml + 2AA) 

- 25aihlhl{Ml + M2 + 2AA,) - 25a:ihlX\Ml + + 2^^,) 

+ '^-5a:,glhl{Ml - 2MMh - Mi)) + Q5a^glhl{Ml - 2M2{A, - M^)) 

404 64 64 

+ ^^7^M,^ - y^?3X^ + ^^7?^?3W + M3^ + M1M3) 

drn^ 1 1 / 

levr^^ = 5a,hlMl - glMl - S^^^'M^^ - -gl^ + " e^as/^' (M^ + A^ 

- SSashlhl {Ml + Ml + 2A,A,) - d^^hlX^ {Ml + Ml + 2A,A^) 
+ 25.3^7?/^^(M,2 - 2Mi(A. - Ml)) + ^-^glMl + ^glMj 

+ 2,glgl{Ml + M^ + M1M2) - (7^^' + ]^glai + ^(72^2 
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- Q6a3hlhl{M^ + + 2A^A,) - 26a3hlX^{M^ + + 2A^Ax) 

- 25a39lhl{M^ - 2MMr - Ml)) + Q5asglhl{M^ - 2M2{A, - M,)) 

+ 156^f M2 + 2gie + 2gla^] . (B.5) 



B.5 Higgs masses 

Similarly, the two-loop RGEs for the Higgs soft masses are 

16vr2^^ = 3h',M^ + X'Ml-glM^,-3glMl + ^gl^ 

+ ^,(^-lSht{M^ + A^,)-6X^{Ml + Al) 

- 3h^X (M2 + + 2 At A,) - 3hl\^ {Ml + Ml + 2 A^a) 

- hlX^ {Ml + Ml + 2ArAy) - 2X^K^ {Ml + Ml + 2A^A^) 

+ ^^7?/i?(M2 - 2Mi(A - Ml)) + lQglhl{Ml - 2M^{At - M3)) 
+ |^^?fMi2 + y^?2'M| + 3(7?(72'(Mi^ + M| + MiM2) 

+ ^,[-l<{Ml + Al)-QK{Ml + Al) 

- 6A^ (M^^ + Al) - 3hlhl {Ml + Ml + 2^^ A) 

- 3h]X^{Ml + Ml + 2AtAx) - 2X^k^ {mI + Ml + 2AaA,) 

- yihl{Ml - 2MMb - Ml)) + 2glhl{Ml - 2M,{A, - M,)) 

fiQ 33 

+ lQ9lhl {Ml - 2Ms{A, - M3)) + ^-^gtMl + ^g^M^ 
+ 3glgl{Ml + Ml + M1M2) - gli' + ^(^iVi + ^gla2 

16n^^ = 2X'Ml + 2K^Ml + -l-,(^-8X\Ml + Al)-lQ^\Ml + Al) 

- eX^hl {Ml + Ml + 2AxAt) - eX^hl {Ml + Ml + 2AxAb) 

- 2X^hl{Ml + Ml + 2AxA^) - 8X\^{mI + Ml + 2AxA,) 

+ 2glX^{Ml - 2MMx ' Mi)) + QglX^Ml - 2M2(Aa - M^))] . (B.6) 
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B.6 Additional parameters of the general NMSSM 

The two-loop RGEs for the SUSY conserving fi and fi' terms are 



32n'^ = f,{3h^,+3hl + hl + 2\'-gl-3gl) + ^^[-9ht-9ht 

4 



- 3/1^ - 6A^ - 6hX - 3hl\^ - 3hlX^ - hlX' - AX't,^ + ^glh] 

- '^glhl + 2glhl + IQglhl + IQglhl + + ^-gl + 3glgi) , 
16vr2^ = j_,'(^2X' + 2K'] + J^(-4X'-8K'-Qh^X' 



dt \ J 167r2 

- 6hlX' - 2hlX' - 8X^n' + 2glX' + Qg^X^^ . (B.7) 

For the corresponding soft SUSY breaking terms m| and m'g (see (2.9)) we have 

3271^^ = 3hl {ml + 2^1 At) + 3hl [ml + 2/iAfe) + [ml + 2/iA^) 

+ 2X^ [3ml + 2;Uy4A) + 2AKm'| - gl [ml - 2^Mi) - 3gl [ml - 212M2) 
+ ( - 9/i^ [ml + ifiAt) - 9ht [ml + ifiAb) - 3ht [mj + 4/iv4,) 

- 2A' [7ml + IS/xAa) - 6/1? /i^^ [ml + 2/i( A + A)) 

- 3h^tX^ [5ml + 2/i(3A + A^)) - 3hlX^ (Sm^ + 2/i(3 A + Ax)) 

- hlX^ [5ml + 2fi{3Ar + Ax)) - AX^k^ [ml + 2fi{Ax + A^)) 

- 8A^K (m'l + /xMa) - 8Xk^ (m'| + 

+ ^glh^timl + 2/i(A - Ml)) - l9lhl[ml + 2/i(A - Mi)) 

+ 2gfhl[ml + 2ii{Ar - Mi)) + AglX^[ml - ^iMi) 
+ l2glX^ [ml - /iMa) + 16^^/^t2 j|^2 ^ 2^(At - M3)) 

23 

+ l&glhl[ml + 2/i(A - M3)) + y (/^(m^ - 4/iMi) 
+ Y^?2K - 4/iM2) + 3glgl[ml - 2/i(Mi + M2))) , 
IGtt^^^ = 2A2 (m'l + 2/i'AA) + 4^2 (m'| + /i'A^) + 4AKm2 



+ ( - (m'l + 4/iMA) - 8k' (2m'| + 5/iM.) 

- 8X^t^ (2m'| + /i'(3AA + 2A^)) - QX^h^^ (m'| + 2fi'{Ax + A)) 

- 6X^hl (m'l + 2/i'(AA + Ab)) - 2X^hl (m'| + 2fj.'{Ax + A,)) 

- 8A^K(m^ + fiAx) - 12XKhf[ml + jjAt) - 12XKhl[ml + fiA^) 

- AXnhl [ml + jj^A^) + 4:XKgl [ml - fiMi) + UXngj [ml - /iMa) 

+ 2AVKI + V(^A - Ml)) + 6A2^?2(m'| + 2^'[Ax - M^))) . (B.8) 
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Finally for the singlet tadpole terms, the two-loop RGEs read 



- 3hlX^ - hlX' - AX^k' + gfX^ + 3glX^'^ , 

levr^^ = X''{^s + 2Ax^F)+f^'{^s + 2A^^p) 

+ 2X{ml{Ax + fi') + fi{mjj^ + m^J) + K(m'|(A, + fi') + 2/i'm|) 

- ?>X^h] l^s + 2(A, + A^)^f) - SX'^hl {^s + 2iAx + A,)^p) 

- X''hl{^s + 2{Ax + - 4A\2(e5 + 2(Aa + A^)^f) 

- QXhl (ml (Aa + At + fi') + fi{M^ + At{Ax + /i') + m^^ + m^j) 

- 6A/i^(^m2(AA + A + fi') + ;u(Mfe2 + A(Aa + fi') + m^^ + m^j) 

- 2A/i2 (^m2(AA + A, + + ii{M^ + v4,(Aa + fi') + mj,^ + rn^j) 

- AX^(ml{2Ax + fi') +fi{Ml + Ax{Ax + fi') + + m^j) 

- AX^K (w!l [Ax + A« + /i') + /i' [Ml + (A, + /i') + 2m|) ) 

- 4«:3(^m'|(2A« + /i') + /i'(M2 + + /i') + 2m2)) 

+ Xg\ (Zml [Ax + /i' - Mi) + 2/i (m^^ + - AxM^ - /x'Mi + 2Ml) 

^X[2ip{Ax-M^)^is)) 

+ "iXgl (?>m\ [Ax + /i' - M2) + 2/i (m^^ + m^^ - Ax^h - /i'Ms + 2Ml) 

^X[2ip{Ax-M^)^is)) \ ■ (B.9) 
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Appendix C. Radiative corrections to the Higgs masses 

As explained in Section 3.2, the following steps are required for a systematic calculation 
of the radiative corrections to the Higgs masses in the present approach: i) the gauge and 
Yukawa couplings at the scale Msusy must be determined; ii) the effective potential and the 
Higgs wave function normalisation constants must be computed with an ultraviolet cutoff 
MgusY; and iii) the pole masses have to be evaluated from the effective action. Here we give 
the necessary formulae for the calculation of the radiative corrections to the Higgs masses 
to the order described in Section 3.2. 

C.l Yukawa and gauge couplings at the scale Msusy 

First the Yukawa couplings (we confine ourselves to the quarks of the third family) have to 
be deduced from the pole masses. The running top quark Yukawa coupling at the scale rrit 
is given in terms of the top quark pole mass as^ 



For the computation of /ij(MgusY), we consider only contributions involving potentially 
large logarithms (the leading logarithmic approximation or LLA) as ln(M|ugY/'m^) and 
more, which depend on the masses of the various particles appearing in the loops. For 
instance, a complete SU (2) multiplet of Higgs bosons can have (nearly degenerate) masses 
given by Ma (the mass of the MSSM-like CP-odd scalar, cf. Section 2), if M| = Mp^u ~ 
Ai"^ ^ M|. These states will only contribute to the radiative corrections to ht from 
scales > M\. In addition, charginos and neutralinos can have masses below Msusy 
and induce threshold effects; within the LLA it is sufficient, however, to approximate their 
masses by the diagonal elements of the corresponding mass matrices (2.35) and (2.32). 
Then, pure gaugino and higgsino states appear simultaneously in the Higgs wave function 
renormalisation diagrams, and only the mass of the heavier state appears in the potentially 
large logarithm. (On the other hand, effects from squark/gluino loops can be neglected in 
the LLA, since the squark masses define the ultraviolet cutoff MgusY-) One obtains (with 

^Here we use the MS relation between the pole mass and the running mass. In the DR scheme the 
second factor A/Ztt would read S/Stt, which would decrease ht{mt), and hence the lightest Higgs mass, 
by ~ 1%. Since we have not included subdominant (single) logarithms in the two-loop corrections to the 
effective potential we are not sensitive to the scheme in which the running top quark mass is defined, which 
leads to a theoretical error of > 1% on the mass of the lightest Higgs. 




(C.l) 
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V 



Ks + /i'/2, and the scales of couplings without arguments are arbitrary in the LLA) 



SUSY J 



ht{mt) ( 1 + -^a,(mt)ln(M|usYMt) 



47r 



-4/7 



647r2 



9/1? 



SUSY 



(-9 cos^ + (3 cos^ /3 - + 2 cos^ ln(M|/m; 



-2A^ ln(max(/i,^ff, 4z/'^)/Ms^usy) - ^?^ In (max (/i.^^, Mi^)/Ms^usy) 
-3^2' ln(max(/i,\, M22)/M|usy)] ) • (C.2) 

In the case of /i;,, we assume a given value of the running h quark mass m^yiMz) (typically 
~ 2.9 GeV), from which h\,iMz) can be obtained as in (A. 11). Integrating the one-loop 
(QCD) RGEs, hb{mt) is given by 



and subsequently /i;,(Msusy) by 
hbiMsvsy) = 

x(l + 



23 
12^ 



as{mt) In 



Ml 
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'^hl - ^gf - ^gl ] ln(M|usYK) 



+ (-9 sin^ phi + (3 sin^ (3 - l)hf - 2 sin^ {3hl) \n{Ml/ 
\^ ln(max(/i^g, Au'' 

-3gl ln(max(/i,\, M2)/Msusy) 



2A2 ln(max(/i,\, 4i/'2)/M|usy) - ln(max(/i,\, M2)/MsV) 



(C.4) 

The electroweak gauge couplings are assumed as given at the scale Mz e. g. in the 
on-shell scheme in terms of the Fermi coupling Gp, Mz and Mw as g^ = 4a/2G'^M^, 
g\ = 4\/2Gi?(M| — M^) (the renormalisation scheme is irrelevant in the LLA). For their 
computation at the scale Msusy, we include threshold effects of the potentially heavy 
SU{2) Higgs multiplet with a mass Ma, higgsinos with mass /Xefr, SU{2) gauginos with 
mass M2, and sleptons of a common mass (in the LLA). Squarks with masses ~ Msusy 
do not induce threshold effects at scales below Msusy, and the top quark threshold is not 
considered since ~ Mz inside logarithms in the LLA. 
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The strong gauge coupling at Msusy is computed somewhat more precisely; the one-loop 
RGE is solved exactly below and above the top quark threshold: 



,2/ 



^3'(^susy) = gl{m,){l + ^f^[71n(MsVK) -21n(M3V/M|)])" • (C.8) 

C.2 Higgs wave function renormalisation constants 

Within the present approximation it suffices to consider wave function renormalisation 
constants for the weak eigenstates Hu-, and S. In the case of contributions from quarks 
and leptons of the third family we use an infrared cutoff m^, which simplifies the remaining 
top- and bottom-quark induced corrections to the pole masses in section C.4 below. In the 
case of contributions from gauge bosons, charginos and neutralinos we confine ourselves 
to potentially large logarithms, i.e. we use Mw ~ Mz and the diagonal elements of the 
chargino/neutralino mass matrices inside the logarithms. 

If M\ ^ m^, we include the corresponding threshold effects: the wave function renor- 
malisation constants of heavy Higgs states on their mass shell ~ M\ do not receive con- 
tributions from momenta < M\ in the LLA, and subsequently this effect is translated into 
the wave function renormalisation constants for the weak eigenstates via a rotation 

by the angle (3. Then we obtain for Zh^, Z^^ and Zs in the Landau gauge (we recall 
z/' = + /i72): 

Zu^ = 1 + (3/^? ln(M|usY/'^?) -Ii9l + ^gl) HMksy/Ml) 



,2 



+ cos^ P{3hl + K- ^h'i) HMl/mi) + | ln(M|usY/ max(/i,\, M^)) 
+ M ln(M|usY/ max(/.,\, M^)) + ln(MsV/ max(/i,\, Au"))) , (C.9) 
Zh, = l + j^{i3hl + hl)\niMi^sy/m',)-^igf + 3gl)HMl^,y/M^ 



2 



sin^ f3{3h', -hi- 3hl) HM\/ml) + f HMIvsyI max(/x,\. Ml)) 

^gl 



+ ^ln(MsW/max(/i,\,M|)) + A2ln(MsW/max(/i,\,4z.'2))) , (C.IO) 
(a^ ln(M|usY//^cff) + ln(M|usY/4i^")) • (Cll) 



^7r2 



These wave function renormalisation constants will multiply the corresponding kinetic 
terms in the effective action, which results from the addition of quantum effects with 

< M|ugY to the "tree level" Lagrangian, wherein all couplings are defined at = 
Mlugy. In the next subsection, the contributions of quantum effects with Q"^ < M^^^y to 
the Higgs mass matrices are considered, in which - to start with - the Higgs fields are not 
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(yet) properly normalised. We will denote the Higgs vevs before the rescaling by square 
roots of the wave function renormalisation constants by Vu{Q), V(i{Q) and s{Q)^ which are 
related to the vevs Vd and s of the properly normalised Higgs fields by 



Vu{Q) = vj , Vd{Q) = vj ^Z^^ , s{Q) = s/^s- (C.12) 

In the LLA, and Vd are related to Mz and the to quark and lepton masses as in (2.13) 
and (A. 11). 

C.3 Corrections to the Higgs mass matrices 

According to the philosophy employed here, we denote subsequently by "tree level" mass 
matrices those of Section 2 where i) all couplings are defined at the scale = M|usy, and 
ii) all Higgs vevs are given by Vu{Q), Vd{Q) and s{Q) as defined in (C.12). 

The one-loop corrections to the effective potential are given by the Coleman- Weinberg 
formula (with an ultraviolet cutoff M|usy, and in the DR scheme in agreement with the 
two- loop RGEs for the soft SUSY breaking terms in Appendix B): 



In (mVmIusy) - \ 



(C.13) 



where the couplings and Higgs vevs in the mass matrices M of the particles in the loops 
are still those at the scale Mgugy. As a result of this procedure, the two-loop contributions 
to AV^ff and hence to the Higgs mass matrices - at least in the approximation considered 
here - become quite simple. 

The radiative corrections to the effective potential affect also the three minimisation 
equations with respect to Vu{Q), Vd{Q) and s{Q), which serve to replace the three para- 
meters Tn^j^, m'jj and m| {fJ^cS, Bcs and mg were defined in (2.7) and (2.14)): 



2 I „2 



VuiQ) {ml^ + l^ls + X'vjiQ) + ^-^{vliQ) - vim 
-VdiQ) (As(g)fieff + fhl) + -AKfr,^„(Q) = , 

V^Q) + /^efT + vliQ) + ^-^^iv'diQ) - vim 

-VuiQ) (As(g)fieff + fhl) + ^AV^sMQ) = ' 
s{Q) (ml + m!g + /i'^ + 2k^f + k.A^s{Q) + 2k^s^{Q) 
+\\vl{Q) + vim - 2XKv^m'd{Q) 



Hs + ^F/i' - \v^{Q)vdmAx + /i') + ^AKff,s(Q) = . (C.14) 

The contributions of AV^fj to the Higgs mass matrices are obtained by appropriate second 
derivatives of A\4fr (with respect to Vu{Q), Vd{Q) and s{Q) in the CP-even case), taking 
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care of the modified replaced mass terms: 



AMI = ^AKfr,.,., for 



v^ = VuiQ), VdiQ), s{Q) . (C.15) 

Below we give the corrections to the Higgs mass matrices for the CP-even, CP-odd and 
charged Higgs masses, separately for each particle species in the loops. 

CP-even scalars 

Top and bottom (s) quarks 

The top and bottom squark mass matrices have been given in (2.36) and (2.37) in Sec- 
tion 2. They are diagonalized as in (A. 10) with the help of angles 6t and 6b such that the 
eigenstates have masses M^^ and M^^ with M^^ < M^^ (and corresponding eigenstates 
and Mi^). It is useful to define the following quantities in terms of the mass eigenvalues 
and mixing angles: 

ft = .,0 ^ .,0 \MHn \ I - In 

'Mi+M? f M? 

gt = sin^ -j^ In — | 

^ ^ M? - M? \M? 




et = —rrit sin 29t In 



M? 

ti 




'M? +M? [Ml 
^ Ml - Ml \Ml 



' M? 

-mbSin2^Bln [t7t\ ■ (C-16) 



Ml 



Some of the radiative corrections due to top/bottom squark loops to CP-even, CP-odd 
and charged Higgs masses can be described by a shift of the trilinear soft SUSY breaking 
coupling A\ (at the scale Ml^gy) everywhere in the tree level mass matrices: 

3/72 3/)2 

A, ^ A', = A, + ^A/* + ^,A,h . (C.17) 

The remaining corrections h"^ = /ij(M|ugY) and = hKM^y^y) the CP-even 
mass matrix (2.20) in the weak basis {HdR, Hur, Sr) are given by 
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3/^2 / (Mi Ml 

AMlu = ^l-Alg, + AA,e, + Aml\ni^^ 

(HQ)? 9t , 



327r 



2 



/ 3/?2 3/)2 

3/)2 3/72 
A>l|,i3 = ^At;„(g)(A,(7,-2e,) + ^A2.(g)t;,(g)(4/,-<70, 

A>i|,23 = ^At;,(g)(A,(7,-2e,) + ^A2s(gXg)(4/,-<7,). 

(C.18) 

Additional contributions proportional to the electroweak gauge couplings gf and do not 
generate large logarithms and would be the same as in the MSSM (see [534] and refs. 
therein) . 



Dominant two loop corrections 



As mentioned above, the two-loop corrections to AV^cfr and hence to the Higgs mass matrices 
are relatively simple within the present approach and within the approximation where we 
neglect all terms without two powers of large logarithms; the ones given below can be 
obtained by integrating the relevant RGEs. The only complication arises, as above, from 
diagrams where Higgs bosons are connected to two top or bottom-quark lines which generate 
a (logarithmic) dependence on Ma- As infrared cutoff (in the case of diagrams with b quark 
lines) we use Mz (in the LLA, inside the logarithms), since the running Higgs masses 
at this scale are closer to the Higgs pole masses. Then one obtains 
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(C.19) 



(Contributions ~ g^ cancel.) 
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Dominant slepton contributions 



Contributions proportional to the electroweak gauge couplings gf and from squark loops 
do not give large logarithms, since we assume the squarks to have masses close to the 
ultraviolet cutoff Msusy- Sleptons could be considerably lighter, and here we assume a 
common slepton mass Mj. The corresponding contributions are best expressed in terms of 
the weak mixing angle 6^ and the auxiliary quantity 
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Then we have 



-4:±|Im|(9 sin^ ew + 3 cos^ 6^) In ( % 
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(C.20) 



(C.21) 



Chargino / neutralino contributions 

Here we assume that the gaugino mass terms Mi and M2 are of similar order of magnitude. 
Ml ~ M2 = Ml 2 (inside logarithms). We define the following potentially large logarithms: 

= ln(/i2g/M|usY) , 
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fii/ 



ln(4i.'VM|usY) , 

ln(max(Mi%,/^,2ff)/M|usY) 

ln(max(4i/'2,/i,\)/M|usY) 



(C.22) 



Again, some of the radiative corrections due to chargino/neutralino loops to all CP-even, 
CP-odd and charged Higgs masses can be described by an additional shift of the trilinear 
soft SUSY breaking coupling Ax on top of the corrections in (C.17): 

^ ;{glMi + 3glM2)LM,,. (C.23) 
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The remaining contributions are given by 
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(C.24) 



Higgs contributions 



The contributions from Higgs loops are quite involved even in the LLA: the eigenvalues of 
the tree level mass matrices (2.22), (2.26) and (2.28) can be anywhere between M| and 
M|ugY, and consequently Higgs loops can generate many different possibly large logarithms. 
Here we confine ourselves to the loop corrections to the diagonal element of the CP-even 
mass matrix after the rotation by the angle /3, which corresponds at tree level to (2.23), 
the upper bound on the lightest eigenvalue of A1|. These Higgs induced loop corrections 
will be denoted as Aniggs- Subsequently the rotation by the angle (5 can be reversed, and 
the contributions can be written in terms of corrections to the upper left 2x2 CP-even 
submatrix. (This is the procedure employed equally in [124]; the present expression for 
Aniggs differs however from the one in [124], since the conventions for the tree level mass 
matrices are different.) Hence, first we define 
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+ |Mi^(l _ sin4 2/3) + M±M)!A2(sin4 2/3 + sin' 2/3 - 2) 
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+ X\~ sin^ 2/3 - ^ sin^ 2/3 + 1) + X\X + k sin 2/3f^ In {{M%u + -^p,22)/^l) 
- K.^'^^'^'^T + /.Sin 2/3) f ^^(1 - sin^2/3) - A^(2 - sin^ 2/3) 
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(C.25) 



^•^1,11 = ^Higgs cos^ /3 , 

AM%22 = ^Higgs sin^ /3 , 

AA4|^i2 = Amggs sin /3 cos /3 . (C.26) 
Gauge boson contributions 

These are relatively simple in the Landau gauge; first we define 

A/fS/ 2 I 2\ 

Acauge = ^r!; 2 (-9 + 12sin2^H^ - 6sin^^vy)ln(M|usY/M|) ; (C.27) 
then we obtain 
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(C.28) 



Rescaling 



After adding up all previous corrections to AA^| (and performing the shifts (C.17) and 
(C.23) of Ax in the tree level expressions for A^|) we obtain the Higgs mass matrix in 
terms of Higgs fields which are not yet properly normalised, since the kinetic terms in the 
effective action are multiplied by the wave function renormalisation constants. Hence, the 
Higgs mass matrix elements have to be rescaled by appropriate powers of wave function 
renormalisation constants Zi, i = H^, H^, S: 

M'l,, = M%^I^Za, . (C.29) 

Herewith also dominant double logarithms (two-loop terms corresponding to an RG im- 
proved one- loop calculation) are generated on top of the ones already present in A^c. 

I„ the basi. - ln,„,n^n,S.), M'l, . diagona,i.ed by a„ orthogo„al^3 X 3 

matrix Sij in analogy to the tree level procedure (A.l). The eigenvalues of ■M.'s,ij (ordered 
in mass) will be denoted as M''s,i- This concludes the evaluation of the radiative corrections 
to the CP-even Higgs mass matrix; the pole masses will be computed in the subsection C.4. 
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CP-odd scalars 



Here we confine ourselves to radiative corrections of tlie order hff^ and h^^ to tlie mass 
matrix, except for tfie cfiargino/neutralino contributions ~ (7^, wfiicli are included as well. 
These radiative corrections are quite simple: all effects of the order h^^ are included by the 
shift (C.17) of A\, all effects induced by chargino/neutralino loops by the additional shift 
(C.23) of Ax, and all contributions of the order g"^ h^ f^ - as before - by the rescaling by the 
wave function renormalisation constants where it is sufficient to consider the terms ^ h^j^ 
in Zh^^h^ (which we will not indicate explicitly). 

The only complication originates from the fact that the rescaling by the wave function 
renormalisation constants is performed for the 3x3 mass matrix before the identification 
of the Goldstone mode, and has to be translated into the 2x2 mass matrix after dropping 
the Goldstone mode. The result for the rescaled CP-odd mass matrix elements A^'p^jj in 
the basis specified before (2.26) is 



M%, = A(B„j + 3Ks(Q))HiiM^iM_3K^,s(Q)_2m|-K,^'s(Q) 



VuiQ) , VdiQ) 



s{Q)J s{Q) ' 

>tpM2 = - '^^s{Q) - ii')^vl{Q)/ZH, + vI{Q)/Zh^ , (C.30) 

where /Xefr, -Bcs and m\ are defined as in (2.7, 2.14), but in terms of parameters (and s) at 
the scale Q = Msusy and with A'l instead of A^. 

The conventions for the diagonalization of the radiatively corrected CP-odd mass matrix 
are the same as those at tree level corresponding to (A. 2 - A. 5) above, leading to two 
eigenvalues A^'pj (ordered in mass). 

Charged scalar 

The precision of the radiative corrections is the same as in the CP-odd case. However, both 
the corrections of the order hf^^ and those induced by chargino/neutralino, gauge boson and 
slepton loops give rise to some additional terms on top of the shifts of Ax- Altogether one 
obtains 



(Q) , v,{Q) 



^Q^2 1 ■"b "n--SUSY/ ■■"f / ^^2 "n--SUSY/ 

^ T^l^l^ ln(M|usY/M|) + 2igy, - gt)LM,,^ , (C.31) 
where Lm^^ was defined in (C.22). 
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C.4 Pole masses 



In order to obtain the Higgs pole masses with the precision defined in Section 3.2 (exact to 
the order hf, hf, g'^h^, etc., but only in the LLA in the case of four powers of electroweak 
gauge couplings and Yukawa couplings A, k) it suffices to consider top- and bottom-quark 
corrections to the pole masses, which can be deduced from [534]. There, the corrections 
An(p^, Q"^) to the inverse tree level propagators n(p^) = — M^,.^^ were computed. In 
order to make contact with our approach, we first have to write 

An(p2, ^ ^n^p2^ ^ = o, q') . (c.32) 

The quantity AIl{p^ = 0, Q"^) corresponds to (the negative of) our radiative corrections 
to the masses squared before the rescaling by the wave function renormalisation constants. 
The quantity An(p^, Q"^) proportional to still contains potentially large logarithms. We 
have already included these large logarithms in our wave function renormalisation constants 
in (C.9) - (C.ll), but there remains a missing part 6Il{p^) which is defined as 

^(/) + A^(p^g2) = p' + Auip',Q') - mI^ + Auip' = o,q') 

= Z{Q^)p^-M^ + 6U{p^) (C.33) 

where = M^^.^^ — An(p^ = 0, Q^) are the radiatively corrected masses squared before 
the rescaling by the wave function renormalisation constants. Hence we obtain 

5Uip') =p' + Af[{p\ _ z{Q^) p" = An(p2, m?) (C.34) 

where we have used the fact that the relevant terms ~ h^j^in our wave function renormalisa- 
tion constants Z{Q'^) in (C.9) - (C.ll) were computed with an infrared cutoff m^. Hence, 
no large logarithms appear in the remaining pole corrections 5H(p^), which justifies the 
present method. Finally, one obtains for the pole mass from the last expression in (C.33) 
within the present approximation 

M^,, = M'/ZiQ') - 5U{MyZ{Q')) (C.35) 

where the first term corresponds to the results in Section C.3 (where M'^/Z{Q'^) = M'^), 
and 611 can be deduced from [534] according to the formulae above. In the case of mass 
matrices, M^^j^ is diagonalized by the same rotation which diagonalizes M'^ within the 
present approximation. 

In the CP-even case, one obtains in terms of the eigenvalues A^'l j and of the elements 
of the diagonalization matrix Sij 

- ^,Sf,[M%M^V^l)+{M'%-Aml)B{M%,,ml)) (C.36) 

where the logarithm ln(m^/m^) in the second line originates from the infrared cutoff in 
Zh^ in (C.IO), and leads to a finite pole mass correction for irib — )■ 0. B{M'^,m?) is defined 
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by 
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A similar expression holds for the CP-odd Higgs pole masses in terms of the eigenvalues 
and the elements of the diagonalization matrix Pjj defined in (A. 3): 
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For the charged Higgs pole mass one obtains (neglecting mf inside logarithms) 

3 



M'l + 



167r2 



{ht cos' ^ + hi sin' /3) 



(C.39) 



M 



/2 



m^ 



In 



M 



/2 



+ (mj + ml) 



mi 
Ml 



mt 



In 



mr 



/2 



+4:htmthhmh sin /3 cos /3 



mr 



M' 



mf 



In 



mr 



± 



mr 



+ 1 



117 



Appendix D. Public computer tools 



NMSSMTools [52, 88, 345] allows to compute the Higgs and sparticle spectrum, mixing 
angles and Higgs branching fractions. NMSSMTools includes NMHDECAY (where the 
parameters can be specified at the weak/SUSY scale), NMSPEC (where the parameters 
can be specified at the GUT scale) and a version where the soft terms satisfy GMSB-like 
boundary conditions. A connection to MicrOMEGAS is provided. 
Link: http://www.th.u-psud.fr/NMHDECAY/nmssmtools.html 

CalcHEP allows to compute cross sections in the NMSSM (used, e. g., by MicrOMEGAS). 
Link: http: //theory.sin p. msu.ru/~pukhov/ca lchep.htm I 

MicrOMEGAS [502,503] allows to compute all relevant cross sections for LSP annihilation 
and coannihilation, integrates the density evolution equations, determines the LSP relic 
density and detection cross sections. It can be called from NMSSMTools. 
Link: http:/ /wwwiapp. in2p3.fr/lapth/micromegas 

Spheno [535] allows to compute the sparticle two-body decay amplitudes in the NMSSM; 
a computation of the Higgs and sparticle spectrum is in progress. 
Link: http://ific.uv.es/~porod/SPheno.html 

WHIZARD is a multi-purpose Monte Carlo event generator, into which the NMSSM has 

been implemented [536]. 

Link: http://projects.hepforge.org/whizard 

Superlso [537] allows to compute fiavour physics observables in the NMSSM such as i? — )■ 
Xs7, the isospin asymmetry of 5 — )■ A'*7, and the branching ratios Bg — )■ B ^ t u^-, 

B ^ D T K ^ ^u^, Dg ^ T Vj. and ^ n v^. 
Link: http://superiso.in2p3.fr 
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